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ABSTRACT 


The  primary  objective  of  this  thesis  was  to 
develop  a  computer  analysis  method  to  facilitate  the 
analysis  of  the  tensile  splitting  test  for  low  temperature 
tensile  properties  of  asphaltic  concrete.  The  computer 
programs  developed  within  this  thesis  were  used  to  analyze 
the  results  of  the  tensile  splitting  test  for  cored  pave¬ 
ment  specimens  obtained  during  the  early  service  life  of  a 
four-lane  divided  freeway  in  central  Alberta.  In  addition, 
an  analysis  of  alternate  laboratory  preparation  methods , 
and  a  preliminary  investigation  into  the  effect  of  density 
upon  the  low  temperature  tensile  properties  of  asphaltic 
concrete  were  undertaken  within  the  scope  of  this  thesis. 

The  tensile  splitting  test  method  consists  of 
loading  a  4  -  inch  diameter  asphalt  concrete  cylinder  via 
loading  strips  across  a  diameter,  in  a  compression  testing 
frame  and  within  a  controlled  temperature  chamber  maintained 
at  a  constant  low  temperature.  Output  signals  from  a  load 
cell  and  two  series  connected  linear  variable  differential 
transformers  (attached  to  opposite  ends  of  the  specimen) 
are  monitored  on  a  two  channel  recorder. 

Five  computer  programs  were  established  to  analyze 


in 


the  tensile  splitting  test  results.  The  programs  vary  from 
the  basic  analysis  of  computed  parameters  measured  through¬ 
out  the  duration  of  the  test ,  to  the  Calcomp  plots  of  stress 
versus  strain  and  stiffness  versus  strain  for  all  the  test 
specimens  within  a  test  series. 

The  results  of  the  analysis  of  the  highway  test 
project  indicated  a  definite  correlation  between  the  low 
temperature  transverse  cracking  of  asphalt  pavements,  and 
the  failure  stress,  failure  strain,  and  failure  stiffness 
of  the  asphaltic  concrete  as  measured  by  the  tensile  split¬ 
ting  test  at  0°F.  The  analysis  of  the  alternate  laboratory 
specimen  preparation  methods  indicated  that  the  75  blow 
impact  compaction  yields  a  density  equivalent  to  that 
obtained  by  kneading  compaction.  Furthermore,  a  decrease 
in  the  density  reduces  the  failure  stress  and  failure 
stiffness,  but  increases  the  magnitude  of  the  failure  strain 

The  principal  conclusion  from  this  research  is 
that  the  established  computer  programs  permit  a  thorough 
study  of  the  low  temperature  tensile  properties  of  asphaltic 
concrete.  The  major  recommendation  from  this  study,  is  that 
in  the  design  and  evaluation  of  asphalt  pavements,  subjected 
to  low  temperatures,  the  tensile  splitting  test  and  analysis 
methods  should  be  adopted  to  determine  the  low  temperature 
tensile  characteristics  of  the  asphaltic  concrete. 
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CHAPTER  I 


INTRODUCTION 

Low  temperature  transverse  cracking  of  asphaltic 
concrete  pavement  surfaces  has  been  under  intensive  study 
during  recent  years.  Agencies  involved  in  the  design, 
testing  and  evaluation  of  pavement  components  for  use  in 
areas  having  a  significant  amount  of  sub-freezing  weather, 
have  recognized  that  the  development  of  transverse  cracks 
may  lead  to  a  rapid  loss  in  pavement  serviceability.  In 
order  to  determine  the  relative  susceptibility  of  a  pave¬ 
ment  to  low  temperature  transverse  cracking,  the  majority 
of  the  research  has  involved  tests  on  the  asphalt  cement 
and/or  the  asphaltic  concrete  under  various  service  and 
laboratory  conditions. 

STATE  OF  INVESTIGATIONS 

This  thesis  is  one  of  a  series  of  continuing 
research  programs  at  The  University  of  Alberta  concerned 
with  the  low  temperature  tensile  properties  of  asphalt 
cement  and  asphalt  concrete  cylinders.  The  previous 
research  conducted  by  Gillespie  (1966),  Christison  (1966) 
and  Hahn  (1967)  involved  the  development  of  a  laboratory 
test,  to  be  used  at  the  mix  design  stage,  for  predicting 
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the  cracking  behavior  of  a  given  asphalt  concrete  pavement. 
This  test  method  has  been  described  in  detail  by  Anderson 
and  Hahn  (1968).  The  research  work  conducted  within  the 
scope  of  this  thesis  involved  refinements  of  the  testing 
technique  as  well  as  the  establishment  of  computer  programs 
for  the  analysis  of  the  tensile  splitting  test  results.  In 
addition,  the  thesis  includes  the  further  application  of 
the  tensile  splitting  test  to  field  cores  at  various  temp¬ 
eratures  and  ages  for  different  asphalt  supplies ,  as  well 
as  laboratory  specimens  prepared  by  both  kneading  and  impact 
compaction  methods  and  tested  at  various  temperatures  and 
eompactive  efforts. 

OBJECTIVES  OF  THE  THESIS 

From  the  results  of  the  initial  Investigations  of 
the  use  of  the  tensile  splitting  test  for  asphaltic  concrete, 
further  investigation  was  warranted  and  the  primary  objective 
of  this  research  was  to  develop  a  computer  analysis  program 
or  series  of  programs  to  facilitate  the  analysis  of  test 
results  to  establish  parameters  for  the  low  temperature 
tensile  properties  of  asphaltic  concrete. 

Secondary  objectives  of  this  research  were  as 

follows : 

1.  To  use  the  established  program  or  series  of 

programs  to  demonstrate  the  trends  in  the  low 
temperature  tensile  properties  and  cracking 

frequency  during  the  early  service  life  of  a 
particular  highway  pavement. 


3. 


2.  To  conduct  an  analytical  comparison  of  alter¬ 
nate  laboratory  specimen  preparation  methods 
(i.e.,  kneading  and  impact  compaction)  in 
order  to  determine  characteristic  properties 
as  well  as  low  temperature  tensile  properties 
of  the  asphaltic  concrete. 

3.  To  conduct  a  preliminary  analysis  to  determine 
the  effect  of  compaction  upon  the  low  temper¬ 
ature  tensile  properties  of  asphaltic  concrete. 

LIMITATIONS  OF  THE  THESIS 

Some  of  the  more  significant  limitations  of  this 
thesis  are: 

1.  Inherent  difficulties  of  using  the  tensile 
splitting  test  for  asphalt  concrete  cylinders 
which  respond,  at  the  rate  of  loading  used,  as 
a  visco-elastic  material. 

2.  The  use  of  only  one  rate  of  loading  in  the 
testing  program. 

3.  The  variations  in  material  placed  under 
construction  or  field  conditions. 

4.  The  investigation  does  not  attempt  to  explain 
the  behavior  of  the  asphalt  concrete  cylinders 
in  terms  of  the  properties  of  the  binder. 
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ORGANIZATION  OF  THE  THESIS 

Chapter  II  provides  a  description  of  the  test 
method  and  analysis  of  the  tensile  splitting  test  as  applied 
to  asphalt  concrete  cylinders.  It  also  includes  a  brief 
review  of  the  development  and  previous  work  with  the  tensile 
splitting  test  at  The  University  of  Alberta,  as  well  as  a 
brief  description  of  the  computer  programs  developed  to 
facilitate  the  analysis  of  the  tensile  splitting  test  to 
establish  parameters  for  the  low  temperature  tensile 
properties  of  asphaltic  concrete. 

Chapter  III  is  devoted  to  the  application  of  test 
results.  It  provides  an  example  of  how  the  analysis  pro¬ 
grams  may  be  used  to  establish  trends  of  low  temperature 
tensile  properties  and  cracking  frequency  on  a  four-lane 
divided  freeway  facility  in  central  Alberta.  In  also 
includes  a  brief  discussion  of  the  methods  used  to  compare 
the  characteristic  properties ,  as  well  as  the  low  temperature 
tensile  properites ,  of  laboratory  specimens  prepared  by 
kneading  and  impact  compaction.  In  addition,  there  is  a 
brief  summary  of  the  preliminary  investigation  into  deter¬ 
mining  the  effect  of  compaction  upon  the  low  temperature 
tensile  properties  of  asphaltic  concrete. 

Chapter  IV  contains  the  conclusions  drawn  from 
this  investigation,  the  recommendations  presented  as  a 
result  of  this  investigation,  and  the  recommendations  for 
further  study. 
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CHAPTER  II 


TEST  METHODS  AND  ANALYSIS 

This  thesis  is  one  of  a  series  of  continuing 
research  programs  at  The  University  of  Alberta,  concerned 
with  the  low  temperature  tensile  properties  of  asphalt 
cement  and  asphalt  concrete  cylinders.  Previous  authors, 
Gillespie  (  1966  ),  Christison  (1966  )  and  Hahn  (.1967)  have 
presented  extensive  reviews  of  work  done  on  the  tensile 
behavior  of  asphalt  cement  and  asphaltic  concrete  ,  •  and 
therefore  a  literature  review  will  not  be  repeated  in  this 
report.  The  information  provided  within  this  chapter  only 
provides  a  brief  review  of  the  main  points  of  the  tensile 
splitting  test,  noting  any  revisions,  along  with  a  brief 
description  of  the  computer  programs  established  for  the 
analysis  of  the  tensile  splitting  test  results. 

TENSILE  SPLITTING  TEST 

The  theory  involved  in  the  tensile  splitting  test 
is  described  in  reports  by  Gillespie  (1966),  Christison 
(1966)  and  Hahn  (1967).  The  method  of  test,  materials,  and 
apparatus  are  described  in  detail  in  Appendix  A,  Sections 
1.1  to  6.4.4. 

Briefly,  the  tensile  splitting  test  consists  of 
loading  a  4  -  inch  diameter  asphalt  concrete  cylinder  via 
loading  strips  across  a  diameter,  in  a  compression  testing 
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frame  and  within  a  controlled  temperature  chamber  maintained 
at  a  constant  low  temperature.  Output  signals  from  a  load 
cell  and  two  series  connected  linear  variable  differential 
transformers  (attached  to  opposite  ends  of  the  specimen) 
are  monitored  on  a  two  channel  recorder.  The  analysis  of 
the  results  from  the  tensile  splitting  test  may  vary  from 
the  simplest  analysis  as  outlined  in  Appendix  A,  Section  7.1., 
to  one  of  the  computer  analysis  methods  described  in  this 
chapter . 


DEVELOPMENT  OF  COMPUTER  PROGRAMS 

In  previous  work  by  Hahn  (1967)  the  output  signals 
were  monitored  on  an  X-Y  plotter.  The  load  and  strain 
coordinates  for  each  specimen,  within  a  test  series,  were 
then  supplied  as  input  for  a  computer  program  to  compute 
and  tabulate  the  coordinates  of  the  average  stress-strain 
curve  for  that  test  series. 

The  emphasis  placed  upon  the  use  of  the  computer 
for  the  analysis  of  the  tensile  splitting  test  is  based 
upon  the  improvements  that  may  be  realized  in  the  data 
interpreation .  In  the  initial  stages  of  this  research, 
the  computer  program  was  established  for  use  on  The  University 
of  Alberta  IBM  360-67,  to  output  sufficient  data  to  identify 
the  test  specimen,  and  to  list  computed  values  of  the  various 
parameters  at  several  stages  throughout  the  tensile  splitting 
test.  The  next  major  revision  in  the  computer  analysis 
involved  changes  so  that  a  statistical  summary  of  all 
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properties  and  computed  parameters  would  be  available  at 
the  completion  of  the  analysis  of  the  last  specimen  within 
a  test  series.  The  final  stage  of  the  computer  analysis 
development  involved  the  establishment  of  programs  to  permit 
computer  plotting  of  various  parameters  with  the  use  of  The 
University  of  Alberta  Calcomp  Plotter,  model  770/663. 

Five  computer  programs  are  included  within  the 
report  and  are  described  in  detail  in  Appendix  A.  The 
following  section  is  included  to  serve  only  as  a  guide  to 
the  use  of  the  various  programs  and  is  therefore  brief  in 
nature . 

BRIEF  DESCRIPTION  OF  COMPUTER  PROGRAMS 

All  programs  are  established  so  that  the  printed 
output  comprises  one  page  per  specimen,  followed  by  a 
statistical  summary  for  the  complete  test  series.  The 
single  page  for  each  specimen  includes  a  summary  of  the 
properties  and  identification  of  the  test  specimen,  as  well 
as  various  computed  parameters  at  several  stages  throughout 
the  tensile  splitting  test.  The  statistical  summary  includes 
the  mean,  standard  deviation  and  coefficient  of  variation  of 
various  parameters  within  the  test  series. 

BASIC  COMPUTER  ANALYSIS 

In  order  to  use  the  "Basic  Computer  Analysis"  for 
the  tensile  splitting  test  the  only  details  required  are 
the  geometry  of  the  test  specimen  and  the  load-deformation 
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trace  from  the  two  channel  recorder.  The  parameters 
considered  within  the  output  of  this  program  are  the 
stress,  strain,  stiffness  of  mix  and  time. 

STRESS  is  the  tensile  stress  in  pounds  per  square 
inch  computed  using  a  formula  based  on  the  theory  of 
elasticity.  STRAIN  is  the  horizontal  strain,  in  inches 
per  inch,  due  to  the  biaxial  stress  condition.  STIFFNESS 
OF  MIX  is  the  stiffness,  in  pounds  per  square  inch, 
computed  by  considering  the  ratio  of  the  tensile  stress 
to  the  tensile  strain  computed  with  an  assumed  Poisson’s 
ratio  of  0.33.  TIME  is  the  elapsed  time,  in  seconds, 
from  "zero  time"  to  any  stage  within  the  execution  of 
the  tensile  splitting  test.  "Zero  time"  is  arbitrarily 
established  as  that  point  where  the  plywood  loading 
strips  have  compressed  sufficiently  to  allow  a  uniform 
load  application  upon  the  specimen  (See  Appendix  A, 

Section  7.2.7.).  The  "Basic  Computer  Analysis"  program 
is  described  in  detail  in  Appendix  A,  Section  7.3. 

A  sample  output  for  the  results  of  one  specimen 
analyzed  by  the  "Basic  Computer  Analysis"  program  is 
provided  in  FIGURE  2.1.  In  addition,  a  sample  statistical 
summary  from  the  "Basic  Computer  Analysis"  program  is 
provided  in  FIGURE  2.2. 

DETAILED  COMPUTER  ANALYSIS 

If  data  are  available  to  determine  the  volumetric 


■ 


9  . 

relationships  of  the  constituents,  the  "Detailed  Computer 
Analysis"  program  may  be  used.  The  program  is  intended 
for  analysis  of  test  results  in  a  manner  similar  to 
Heukolom  (1966)  involving  the  stiffness  of  bitumen  and 
the  stiffness  of  mix.  Furthermore,  the  analysis  involves 
the  determination  of  the  toughness,  i.e.  the  area  under 
the  stress-strain  curve,  as  well  as  the  determination  of 
the  work  input  per  unit  volume  using  the  work-energy 
concept  (Breen  and  Stephens,  1966).  The  "Detailed 
Computer  Analysis"  program  is  described  in  detail  in 
Appendix  A,  Section  7.4. 

A  sample  output  for  the  results  of  one  specimen 
analyzed  by  the  "Detailed  Computer  Analysis"  program 
is  provided  in  FIGURE  2.3.  In  addition  a  sample 
statistical  summary  from  the  "Detailed  Computer  Analysis" 
program  is  provided  in  FIGURE  2.4. 

COMPUTER  ANALYSIS  WITH  STRESS -STRAIN  PLOT 

The  printed  output  for  this  program  is  the  same 
as  that  for  the  "Basic  Computer  Analysis"  ;  but  the 
additional  information  is  provided  in  the  form  of  an 
8-1/2  by  11  -  inch  plot  of  stress  versus  strain  for  all 
the  specimens  within  a  test  series  plotted  on  a  single 
graph.  This  program  is  described  in  detail  in  Appendix 
A ,  Section  7.5. 

A  sample  plot  from  the  "Computer  Analysis  with 
Stress-Strain  Plot"  program  is  provided  in  FIGURE  2.5. 
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COMPUTER  ANALYSIS  WITH  STIFFNESS-STRAIN  PLOT 

The  output  for  this  program  is  similar  to  the 
"Computer  Analysis  with  Stress-Strain  Plot"  except  for 
the  fact  that  the  plot  will  be  stiffness  of  mix  versus 
strain  for  the  complete  test  series.  This  program  is 
described  in  detail  in  Appendix  A,  Section  7.6. 

A  sample  plot  from  the  "Computer  Analysis  with 
Stiffness-Strain  Plot"  program  is  provided  in  FIGURE  2.6. 

COMPUTER  ANALYSIS  WITH  STRESS -STRAIN  AND  STIFFNESS-STRAIN 
PLOTS 

This  program  allows  a  printed  output  as  previously 
described,  as  well  as  a  plot  of  stress  versus  strain 
adjacent  to  a  plot  of  stiffness  of  mix  versus  strain  for 
the  same  test  series.  For  convenience  of  plotting,  and 
also  to  allow  valuable  comparisons,  two  such  combinations 
are  plotted  on  a  17  by  22  -  inch  plot.  Samples  of  this 
form  of  output,  to  a  reduced  scale,  are  given  in 
Appendices  B,  C  8  D.  This  program  is  described  in 
detail  in  Appendix  A,  Section  7.7. 

A  sample  plot  from  the  "Computer  Analysis  with 
Stress-Strain  and  Stiffness-Strain  Plots"  program  is 
provided  in  FIGURE  2.7. 
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FIGURE  2.1  Sample  Output  for  The  Results'  of  One  Specimen  Analyzed  by  the  "Basic 

Computer  Analysis"  Program. 


SUPPLY  NUMBER  5268  L  WATSONVILLE  AGG.  COMPACTION!  :  75  BLOWS 
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IGURE  2.2  Sample  Statistical  Summary  From  the  "Basic  Computer  Analysis"  Program 


SUPPLY  NUMBER  5268  8  WATSONVILLE  AGG.  COMPACTION  :  75  8L0WS 
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FIGURE  2.3  Sample  Output  for  the  Results  of  One  Specimen  Analyzed  bt  the  "Detailed 
Computer  Analysis"  Program 
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FIGURE  2.5 


Sample  Plot  from  the  'Computer 
Analysis  with  Stress-Strain  Plot" 
Program 
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FIGURE  2.6 


Sample  Plot  from  the  "Computer 
Analysis  with  Stiffness-Strain 
Plot"  Program 
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Sample  Plot  from  the  "Computer  Analysis  with 
Stress-Strain  and  Stiffness-Strain  Plots"  Program 


FIGURE  2.7 


CHAPTER  III 


APPLICATION  OF  TEST  RESULTS 

This  chapter  is  structured  as  two  separate  ex¬ 
amples  of  the  application  of  test  results  obtained  from 
the  tensile  splitting  test. 

The  first  part  is  the  analysis  of  test  results 
of  the  tensile  splitting  test  conducted  on  asphaltic  con¬ 
crete  cores  from  a  four-lane  divided  freeway  facility  in 
central  Alberta.  Cores  were  obtained  at  construction  as 
well  as  at  the  following  service  lives:  18  months,  32  months 
and  34  months.  The  cores  were  subjected  to  the  tensile 
splitting  test  at  test  temperatures  of  20°,  10°,  0°  and 
-10°F. 

The  second  part  is  the  analysis  of  alternate 
laboratory  specimen  preparation  methods  to  determine  the 
impact  compaction  required  to  yield  specimens  equivalent  to 
that  obtained  by  kneading  compaction.  Comparisons  were  made 
on  the  basis  of  air  voids  and  unit  weight,  as  well  as  the 
low  temperature  tensile  properties  of  the  asphaltic  con¬ 
crete.  Other  information  included  in  the  second  part  of 
this  chapter,  is  a  summary  of  the  preliminary  investigation 
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into  the  effect  of  compaction  upon  the  low  temperature 
tensile  properties  of  asphaltic  concrete. 


PART  A:  HIGHWAY  2-D-2/1  &  2/2  TEST  PROJECT 

During  the  summer  of  1966  the  Alberta  Department 
of  Highways  in  cooperation  with  the  Alberta  Research  Counci 
undertook  an  extensive  highway  testing  program  as  part  of 
the  final  asphaltic  concrete  paving  of  a  four-lane  divided 
freeway  facility  in  central  Alberta  (Shields  et  al.,  1969). 
In  order  to  obtain  additional  information  on  the  possible 
environmental  conditions  contributing  to  low  temperature 
transverse  cracking,  and  the  influence  of  asphalt  source, 
the  Department  of  Highways  decided  to  incorporate  within 
one  paving  contract  three  different  sources  of  200-300 
penetration  grade  asphalt  cement.  The  three  sources  re¬ 
presented  major  suppliers  of  asphalt  in  the  province,  and 
will  be  referred  to  in  the  subsequent  discussions  as 
Asphalt  Supply  Numbers  1,  2  and  3. 

Normal  construction  quality  control  procedures 
of  the  Department  of  Highways  were  maintained  throughout 
the  project.  In  addition,  an  extensive  evaluation  and 
testing  program  was  implemented,  which  was  primarily 
directed  towards  determining  changes  in  the  asphalt  cement. 
As  well  as  tests  run  upon  the  asphalt  cements  and  the 
mixtures ,  field  cores  were  obtained  from  each  of  the  three 
test  sections  for  further  laboratory  testing.  The  majority 
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of  the  laboratory  tests  conducted  on  the  field  cores  were 
of  a  conventional  nature,  but  as  indicated  by  Shields  et  al 
(1969)  certain  of  the  tests  were  specialized  and  would  be 
reported  in  due  course . 

The  tensile  splitting  test  was  one  of  the  special 
tests  referred  to  specifically,  and  the  subsequent  dis¬ 
cussion  outlines  the  relationships  between  asphalt  source, 
cracking  frequency,  and  the  low  temperature  tensile  pro¬ 
perties  of  the  asphalt  concrete  cored  pavement  specimens. 

In  addition  to  this  correlation,  Hahn  (1967)  established 
a  similar  correlation  between  asphalt  source  and  the  low 
temperature  tensile  properties  of  asphalt  concrete  labor¬ 
atory  specimens,  prepared  using  Asphalt  Supply  Numbers  1, 

2  and  3  and  the  design  values  for  each  mix.  A  report  on 
this  work  is  given  by  Anderson  and  Hahn  (1968). 

In  addition  to  the  field  cores  obtained  at  the 
time  of  construction,  cores  were  also  obtained  from  the 
three  test  sections  at  service  lives  of  18,  32  and  34 
months.  Of  the  specimens  obtained,  214  were  subjected  to 
laboratory  tests  involving  density  determinations ,  followed 
by  the  tensile  splitting  test  at  temperatures  of  20°,  10°, 

0°  and  -10°F . 

The  computer  analysis  incorporated  the  application 
of  each  of  the  five  programs  mentioned  in  Chapter  II.  How¬ 
ever,  the  two  programs  of  greatest  significance  for  this 
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study,  were  the  "Detailed  Computer  Analysis"  and  the 
"Computer  Analysis  with  Stress-Strain  and  Stiffness- 
Strain  Plots".  The  output  from  the  Detailed  Computer 
Analysis  for  the  214  specimens  total  298  pages,  therefore 
reference  will  only  be  made  to  the  statisitcal  summaries 
included  at  the  end  of  each  series  rather  than  to  individual 
test  results.  A  sample  of  the  output  from  the  "Detailed 
Computer  Analysis"  program  is  referred  to  in  Appendix  A, 
Section  8.2.2.  The  Calcomp  plots  of  stress  versus  strain 
and  stiffness  versus  strain  for  each  of  the  test  series 
represented  by  the  214  specimens  are  included  in  Appendix 
B. 

ANALYSIS  OF  RESULTS  FROM  TENSILE  SPLITTING  TEST 

Upon  analysis ,  it  was  found  that  the  results  for 
the  18  month  series  exhibited  higher  failure  stresses, 
lower  failure  strains  and  higher  failure  stiffnesses  than 
would  be  expected  by  observing  the  established  trends  from 
the  "At  Construction"  series  to  the  32  and  34  month  series. 
In  attempting  to  explain  this  apparent  anomaly  it  was  noted 
that  the  coring  of  all  other  specimens  was  completed  when 
the  air  and  pavement  temperatures  were  above  freezing.  In 
the  case  of  the  18  month  cores  the  coring  was  conducted 
at  air  temperatures  at  or  below  -30°F.  The  overall  effect 
produced  by  coring  at  such  a  low  temperature  is  not  fully 
understood,  but  the  fact  that  parameters  such  as  failure 
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stress,  failure  strain,  and  failure  stiffness  were  sign¬ 
ificantly  different,  as  compared  to  the  cores  at  service 
lives  both  preceding  and  following  the  18  month  test  series  , 
has  been  considered  as  sufficient  cause  to  warrant  further 
investigation  into  the  effects  of  coring  at  low  temperatures . 
Because  of  the  time  element  involved  that  type  of  reserach 
is  beyond  the  scope  of  this  thesis. 

In  the  plotting  of  trends  given  in  FIGURES  3.1, 

3.2,  and  3.3  points  are  plotted  to  represent  the  various 
parameters  at  construction  as  well  as  at  the  three  service 
lives  of  18  months,  32  months  and  34  months.  The  discussion 
of  trends  will  be  limited  to  the  change  in  parameters  from 
the  "At  Construction"  value  to  the  centroid  of  the  32  and 
34  month  value.  In  addition,  for  the  purpose  of  this  study, 
the  reporting  and  analysis  of  test  results  is  confined  to 
the  specimens  tested  at  0°F.  This  test  temperature  was 
chosen,  due  to  the  fact  that  the  previous  investigation 
of  the  same  materials  by  Hahn  (1967),  involved  the  tensile 
splitting  test  at  only  the  one  temperature  of  0°F. 

The  following  discussions  compare  the  change  in 
parameters  with  time  for  each  of  the  three  asphalt  supplies. 
The  discussion  and  explanation  of  results  are  made  largely 
with  reference  to  FIGURES  3.1,  3.2  and  3.3  augmented  as 
required  by  information  available  from  the  computer  output 
for  the  specimens  or  series  of  specimens  tested. 
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CHANGE  IN  DENSITY  WITH  TIME 


In  addition  to  the  plots  of  density  versus  age 
in  FIGURES  3.1,  3.2  and  3.3  the  information  in  TABLE  I 
provides  a  means  of  comparing  the  effect  of  service  life 
upon  the  unit  weight  and  air  voids  as  compared  to  the 
Marshall  Mix  Design  values.  (Note:  The  densities  and 
air  voids  reported  as  design  values  are  as  reported  by 
Hahn  (1967).  The  densities  and  air  voids  for  the  two 
ages  are  based  upon  the  volumetric  computations  completed 
in  the  "Detailed  Computer  Analysis"  for  the  specimens 
subjected  to  the  tensile  splitting  test.  The  "N" 
values  given  indicate  the  total  number  of  specimens 
tested  at  that  age). 

The  three  asphalt  cements  considered  had  the 
following  "as  supplied"  penetrations  and  viscosities 


,  19  67)  ; 

SUPPLIER 

PENETRATION 
at  7 7°F 

VISCOSITY 
at  140°F 

1 

265 

237 

2 

215 

610 

3 

217 

544 

Employing  current  terminology  (McLeod,  1967) 
Asphalt  Supply  Number  1  is  a  "low  viscosity"  asphalt 
cement,  Asphalt  Supply  Number  2  is  a  "high  viscosity" 
asphalt  cement,  and  Asphalt  Supply  Number  3  is 


TABLE  I 


SUMMARY  OF  DENSITIES  AND  AIR  VOIDS 
:  HIGHWAY  2-D-2/1  8  2/2 


SUPPLIER 

CONDITION 

UNIT  WEIGHT 
( PCF ) 

AIR  VOIDS 

(%) 

1 

Design 

14  4.5 

3.9 

(N=20) 

At 

Construction 

138.6 

8.1 

(  N  =  3  4  ) 

32  8 

34  Month  Ave . 

147.0 

2  .  5 

2 

Design 

144 . 6 

4.2 

(N=l 8 ) 

At 

Construction 

137.0 

9.3 

(N=31) 

32  8 

34  Month  Ave. 

142.0 

6 . 0 

3 

Design 

145.3 

4.1 

(N=  21) 

At 

Construction 

134 . 0 

12.1 

( — 1 
OO 

II 

S 

32  8 

34  Month  Ave. 

141.4 

7.6 

an  ''intermediate  viscosity"  asphalt  cement.  The  increase 
in  density  with  service  life  correlates  well  with  the 
viscosities  at  140°F.  Asphalt  Supply  Number  1,  the 
low  viscosity  asphalt  cement,  demonstrates  both  the 
greatest  magnitude  of  change  in  density  as  well  as  the 
greatest  percent  change  in  density  with  service  life. 

As  can  be  seen  from  a  comparison  of  FIGURES  3.1,  3.2 
and  3.3,  Asphalt  Supply  Number  1  also  exhibits  the 
highest  crack  frequencies .  This  observation  of  the 
rapid  densification  of  a  low  viscosity  asphalt  pavement 
verifies  the  affirmation  that  "The  use  of  softer  grades 
of  asphalt  to  minimize  transverse  cracking  is  questioned 
in  view  of  observed  rapid  densification  of  the  pavement 
surfaces  in  service"  (Shields  et  al,  1969). 

CHANGE  IN  FAILURE  STRESS  WITH  TIME 

Asphalt  Supply  Number  1  exhibits  the  greatest 
change  in  failure  stress  at  the  test  temperature  of  0°F. 
There  is  a  60  per  cent  increase  in  failure  stress  for 
the  length  of  time  considered  for  Asphalt  Supply  Number  1. 
The  Increase  in  failure  stress  over  the  same  period  for 
Asphalt  Supply  Number  2  is  only  6  per  cent;  this  is  con¬ 
sidered  to  be  insignificant  since  the  coefficient  of 
variation  for  the  average  failure  stress  at  each  age 
generally  ranged  from  6  to  15  per  cent.  Asphalt  Supply 
Number  3  exhibits  an  increase  in  failure  stress  of  30 
per  cent  over  the  same  service  life. 
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In  summary,  for  the  data  available  within  this 
study,  the  low  viscosity  asphalt  cement  mix  exhibits 
the  greatest  increase  in  failure  stress,  the  inter¬ 
mediate  viscosity  asphalt  mixture  experiences  a  lesser 
increase  in  failure  stress ,  but  for  the  high  viscosity 
asphalt  cement  mix  there  is  no  appreciable  change  in 
failure  stress  for  the  time  interval  considered. 

CHANGE  IN  FAILURE  STRAIN  WITH  TIME 

The  important  aspects  of  the  failure  strain 
appear  to  be  its  actual  magnitude  as  well  as  its  rate 
of  change  with  time.  The  low  viscosity  asphalt  cement, 
Asphalt  Supply  Number  1,  has  the  lowest  failure  strain 
at  the  time  of  construction  and  at  the  32  to  34  month 
age.  The  intermediate  viscosity  asphalt  cement,  Asphalt 
Supply  Number  3,  has  the  highest  failure  strain  at  the 
time  of  construction  and  at  the  32  to  34  month  age.  It 
is  interesting  to  note  that  the  crack  frequencies  corre¬ 
late  with  this  parameter  in  that  Asphalt  Supply  Number  1 
exhibits  187  cracks  per  mile  on  the  average,  and  Asphalt 
Supply  Number  3  exhibits  an  average  of  88  cracks  per 
mile  . 

CHANGE  IN  STIFFNESS  OF  MIX  AT  FAILURE  WITH  TIME 

It  is  to  be  expected  that  the  trends  in  stiffness 
of  mix  will  directly  reflect  the  changes  that  have 
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already  occurred  in  failure  stress  and  failure  strain. 
However,  the  stiffness  of  mix  consideration  does  allow 
a  single  approach  to  the  trends  that  occur  jointly  in 
the  failure  stress  and  failure  strain. 

The  low  viscosity  asphalt  mix,  Asphalt  Supply 
Number  1,  exhibits  a  185  per  cent  increase  in  stiffness 
over  the  time  interval  considered.  Asphalt  Supply 
Number  2  shows  a  33  per  cent  increase  in  stiffness  over 
the  same  period,  while  Asphalt  Supply  Number  3  demon¬ 
strates  a  95  per  cent  increase  in  stiffness  of  mix  at 
failure,  from  the  time  of  construction  to  the  32  to  34 
month  test  period. 

CHANGE  IN  CRACK  FREQUENCY  WITH  TIME 

The  crack  frequencies  indicated  in  FIGURES  3.1, 
3.2  and  3.3  are  the  average  number  of  cracks  per  mile 
based  on  the  total  number  of  miles  reported  (FIGURE  13, 
Shields  et  al ,  1969  ). 

The  marked  increase  in  cracking  in  all  three 
asphalt  supply  sections  over  the  winter  of  1968-69  may 
be  explained  by  the  fact  that  the  1966-67  and  1967-68 
winters  were  considered  equivalent  to  the  long-term 
average  for  central  Alberta,  but  the  winter  of  1968-69 
was  considered  to  be  one  of  the  most  severe  in  central 
and  northern  Alberta  in  the  past  75  years. 
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SUMMARY  OF  RESULTS  FOR  HIGHWAY  2  TEST  PROJECT 


TABLE  II  provides  a  summary  of  tensile  splitting 
test  information  for  the  Highway  2  Test  Project.  The  design 
values  of  failure  stress  and  failure  strain  were  obtained 
by  the  tensile  splitting  test  at  0°F  on  field  cores  at 
various  service  lives. 

In  conjunction  with  TABLES  I  and  II,  the  analysis 
of  the  tensile  splitting  test  data  for  the  Highway  2  Test 
Project  yields  the  following  results: 

1)  The  increase  in  cracking  frequency  was 
accompanied  by  a  decrease  in  failure  strain, 
and  an  increase  in  failure  stress  and  failure 
stiffness . 

2)  The  increase  in  density  with  service  life  was 
accompanied  by  a  decrease  in  failure  strain, 
and  an  increase  in  failure  stress,  failure 
stiffness,  and  cracking  frequency. 

3)  In  terms  of  a  relative  comparison,  the  pave¬ 
ment  section  with  the  highest  crack  frequency, 
also  had  the  lowest  failure  strain  and  the 
highest  failure  stress  and  failure  stiffness. 

4)  In  terms  of  service  life  Asphalt  Supply  Number 
1,  the  low  viscosity  asphalt  cement  mix,  ex¬ 
hibited  the  greatest  change  in  density,  failure 


TABLE  II 


SUMMARY  OF  TENSILE  SPLITTING  TEST  INFORMATION 
HIGHWAY  2-D-2/1  S  2/2  TEST  TEMPERATURE  0°F 


DESIGN  VALUES  (Hahn,  1967) 

SUPPLY 
NUMBER  1 

SUPPLY 
NUMBER  2 

SUPPLY 
NUMBER  3 

Failure 

Stress 

(psi) 

490 

480 

49  0 

Failure 

Strain 

(in/in) 

.  00074 

. 00151 

.00165 

AT  CONSTRUCTION 

Failure 

Stress 

(psi) 

302 

289 

194 

Failure 

Strain 

(in/in) 

.  00082 

.00105 

.00138 

Failure 

Stiffness  (psi) 

699  ,000 

599  ,000 

336  ,000 

32  MONTHS 

SERVICE 

Failure 

Stress  (psi) 

460 

305 

270 

Failure 

Strain  (in/in) 

. 00056 

.  00102 

.00083 

Failure 

Stiffness  (psi) 

1,679  ,  000 

579  ,000 

623,000 

34  MONTHS 

SERVICE 

Failure 

Stress  (psi) 

505 

310 

263 

Failure 

Strain  (in/in) 

.00047 

. 00061 

.  00082 

Failure 

Stiffness  (psi) 

2 ,296  ,000 

1,071,000 

704,000 

PAVEMENT  PERFORMANCE 

(Number  of  Cracks/Mile) 

After  1st  Winter 

4 

Nil 

Nil 

After  2nd  Winter 

87 

Nil 

4 

After  3rd  Winter 

18  7 

126 

88 

. 


- 
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stress,  failure  strain,  and  failure 
stiffness ,  and  the  highest  cracking 
frequency  per  mile. 

PART  B:  ANALYSIS  OF  ALTERNATE  LABORATORY  SPECIMEN 

PREPARATION  METHODS 

During  the  spring  of  1968,  Anderson  (1968)  pre¬ 
pared  91  laboratory  specimens  at  the  Institute  of  Traffic 
and  Transportation  Engineering,  University  of  California. 

He  used  one  aggregate,  that  being  the  Watsonville  aggregate 
(State  of  California  Grading  1/2"  Max.  Med,)  and  the  follow 
ing  asphalt  supplies  and  grades:  Chevron  60/70,  85/100, 
and  120/150;  Santa  Maria  85/100;  and  Golden  Bear  85/100. 
These  specimens,  prepared  with  strict  control  on  gradation, 
density,  and  asphalt  content,  were  compacted  by  the  knead¬ 
ing  compaction  method,  and  then  subjected  to  the  tensile 
splitting  test  at  temperatures  of  40°,  30°,  20°  and  0°F. 

As  a  further  example  of  the  use  of  the  computer 
analysis  methods,  the  test  results  which  were  available 
were  processed  using  the  "Detailed  Computer  Analysis" 
and  the  "Computer  Analysis  with  Stress-Strain  and  Stiffness 
Strain  Plots".  The  Calcomp  plots  of  stress  versus  strain 
and  stiffness  versus  _strain  for  each  of  the  test  series 
are  included  In  Appendix  C. 


COMPARISON  OF  KNEADING  AND  IMPACT  COMPACTION 


A  correlation  between  kneading  and  impact  com¬ 
paction  must  be  established  if  comparisons  are  to  be  made 
between  the  test  results  of  the  California  specimens  and 
those  tested  at  The  University  of  Alberta.  In  attempting 
to  establish  this  correlation,  laboratory  specimens  were 
prepared  by  means  of  the  impact  compaction  apparatus 
(Soil  test  Model  AP-195)  using  Watsonville  aggregate  and 
Santa  Maria  85/100  asphalt  cement.  Five  specimens,  sub¬ 
sequently  tested  at  0°F,  were  prepared  at  each  of  10,  50 
and  75  blows  at  the  required  mixing  and  compacting  tempera¬ 
tures  as  determined  by  the  temperature-vis cosity  relation¬ 
ship  for  the  Santa  Maria  asphalt  cement. 

The  purpose  of  the  laboratory  tests  was  to  deter¬ 
mine  the  compactive  effort  required,  using  the  impact 
apparatus,  to  achieve  the  same  density  as  the  kneading 
compaction  method.  It  was  found  that  the  density  and  air 
voids  at  75  blows  agreed  favourably  with  the  density  and 
air  voids  obtained  with  the  kneading  compaction  as  shown 
by  the  following  data: 


KNEADING 

IMPACT  COMPACTION 

COMPACTION 

75 

Blows 

50 

Blows 

10 

Blows 

Unit  Weight 
( PCF) 

153.56 

153 . 34 

152.11 

143 . 18 

Air  Voids 
(%) 

4 .46 

4.60 

5.36 

10 .92 

« 


. 


A  second  method  of  determining  which  impact  com- 
pactive  effort  corresponded  to  the  kneading  compaction  in¬ 
volved  a  plot  of  the  failure  stress  versus  the  failure 
strain  on  a  log  plot,  as  given  in  FIGURE  3.4.  It  should  be 
noted  that  the  91  points  plotted  and  referred  to  as  "Cali¬ 
fornia  Asphalt-Kneading  Compaction"  comprise  five  distinctly 
different  asphalt  cements  (including  three  different  asphalt 
suppliers)  and  tensile  splitting  test  temperatures  of  40°, 
30°,  20°  and  0°F.  The  plot  of  all  the  points  for  this  one 
aggregate  and  compactive  method  fall  within  a  unique  envelop 
as  indicated  by  Anderson  (1968). 

The  second  degree  equation  of  the  best  fit  curve 
(APL  Program  Library  Number  115 ,  University  of  Alberta) 
of  all  the  points  for  the  California  asphalts  and  knead¬ 
ing  compaction  is  as  follows : 

Y  =  4,858,000  X2  -  101,500  X  +  727, 

with  a  standard  deviation  of  57  psi,  considering  X  as  the 

independent  variable,  and  Y  as  the  dependent  variable. 

-4 

Where  X  =  failure  strain  X  10  (in/in),  and 

Y  =  failure  stress  (psi)  . 

By  substituting  the  average  failure  strain  for 
each  of  the  10,  50  and  75  blow  series  into  this  equation, 
a  deviation  from  the  best  fit  curve  may  be  determined.  On 
this  basis,  the  average  failure  stress  for  the  10  blow 


‘ 
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series  lies  5.5  standard  deviations  away  from  the  best 
fit  curve.  Similarly,  the  50  blow  series  is  2.9  standard 
deviations  below  the  best  fit  curve,  while  the  75  blow 
series  is  1.3  standard  deviations  below  the  best  fit 
curve. 

In  order  to  complete  the  cross-referencing  of 
asphalt  supplies  and  aggregates  between  Alberta  and  Calif¬ 
ornia  materials ,  three  sets  of  specimens  were  prepared  with 
Watsonville  aggregate  and  Alberta  Asphalt  Supply  Numbers 
1,  2  and  3,  as  described  in  Chapter  II.  The  specimens, 
prepared  using  75  blow  impact  compaction,  were  subjected 
to  the  tensile  splitting  test  at  a  temperature  of  0°F. 

The  failure  stress  and  failure  strain  values 
have  been  plotted  on  FIGURE  3.4.  The  failure  stress  of 
Asphalt  Supply  Number  1,  the  low  viscosity  asphalt  cement, 
lies  2.1  standard  deviations  below  the  best  fit  curve. 
Similarily,  Asphalt  Supply  Number  2,  the  high  viscosity 
material,  lies  0.6  standard  deviations  away  from  the  best 
fit  curve,  while  the  intermediate  viscosity  asphalt  cement, 
Asphalt  Supply  Number  3,  is  only  0.1  standard  deviations 
from  the  best  fit  curve  considered. 


, 
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PERTIES  OF  ASPHALTIC  CONCRETE 

The  results  of  the  previously  mentioned  specimens 
prepared  at  10,  50  and  75  blow  impact  compaction  have  been 
summarized  on  a  plot  in  FIGURE  3.5. 

The  trends  of  the  low  temperature  tensile  pro¬ 
perties  of  the  asphaltic  concrete  are  as  anticipated.  How¬ 
ever,  in  terms  of  magnitude  of  change,  the  7  per  cent  de¬ 
crease  in  density  from  75  blows  to  10  blows  results  in 
a  50  per  cent  loss  in  failure  stress,  a  36  per  cent  increase 
in  failure  strain,  and  a  64  per  cent  drop  in  stiffness  of 
mix  at  failure  when  the  specimens  were  subjected  to  the 
tensile  splitting  test  at  a  temperature  of  0°F. 

Summary  of  Results  for  Alternate  Laboratory  Specimen  Pre¬ 

paration  Metho d s 

The  analysis  of  the  alternate  laboratory  specimen 
preparation  methods  yields  the  following  results: 

1.  The  75  blow  impact  compaction  produced  a  density 
equivalent  to  that  obtained  by  kneading  compaction 
However,  further  research  should  be  undertaken  to 
establish  a  more  definite  trend  of  the  low  tempera 
ture  tensile  properties  to  allow  a  comparison  of 
the  effect  of  the  two  laboratory  specimen  pre¬ 
paration  methods . 
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CHAPTER  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  general  purpose  of  this  investigation  was  to 
develop  and  demonstrate  the  use  of  a  computerized  method  o 
analysis  for  the  tensile  splitting  test  used  to  determine 
the  low  temperature  tensile  properties  of  asphaltic  con¬ 
crete.  The  extent  of  the  achievement  of  this  purpose  will 
be  indicated  in  this  portion  of  the  thesis. 

CONCLUSIONS 

1.  The  computer  analysis  programs  developed  within 
this  thesis  permit  a  thorough  study  of  the  low 
temperature  tensile  properties  of  asphaltic  con¬ 
crete.  The  programs  are  written  with  a  built-in 
flexibility  so  that  they  may  be  used  by  others 
with  a  minimum  of  difficulty,  and  with  the  free¬ 
dom  of  revising  the  test  conditions  as  the  re¬ 
search  or  analysis  may  require. 

2.  Conclusions  from  the  Highway  2  test  project: 

(a)  An  increase  in  cracking  frequency  is  accom¬ 
panied  by  a  decrease  in  failure  strain,  and 
an  increase  in  failure  stress  and  failure 


stiffness . 


- 


(b)  An  increase  in  density  with  service  life 
is  accompanied  by  a  decrease  in  failure 
strain,  and  an  increase  in  failure  stress, 
failure  stiffness,  and  cracking  frequency. 

(c)  In  terms  of  relative  comparisons  the  pave¬ 
ment  section  with  the  highest  crack  fre¬ 
quency,  also  has  the  lowest  failure  strain 
and  the  highest  failure  stress  and  failure 
stiffness . 

(d)  In  terms  of  service  life,  the  low  viscosity 
asphalt  cement  exhibits  the  greatest  change 
in  density,  failure  stress,  failure  strain, 
and  failure  stiffness,  and  the  highest 
cracking  frequency  per  mile. 

Conclusions  from  the  analysis  of  alternate  lab¬ 
oratory  specimen  preparation  methods : 

(a)  The  75  blow  impact  compaction  yields  a  densi 
equivalent  to  that  obtained  by  kneading  com¬ 
paction  . 

(b)  For  a  given  mix,  and  given  test  conditions  a 
decrease  in  the  density  reduces  the  failure 
stress  and  failure  stiffness,  but  increases 
the  magnitude  of  the  failure  strain. 


RECOMMENDATIONS 


1.  In  the  design  and  evaluation  of  asphalt  pave¬ 
ments  subjected  to  low  temperatures,  the  tensile 
splitting  test  and  analysis  methods  should  be . 
adopted  to  determine  the  low  temperature  tensile 
characteristics  of  the  asphaltic  concrete. 

2.  In  future  research  programs,  involving  the  test¬ 
ing  of  field  cores  over  a  particular  service  life, 
the  following  steps  should  be  considered: 

(a)  establish  coring  locations  at  the  time  of 
construction  with  the  use  of  random  numbers , 

(b)  complete  all  future  coring  under  environmental 
conditions  as  close  to  the  original  as  isr 
feas ible , 

(c)  the  location  of  the  cores  to  be  obtained  for 
various  service  lives  should  be  in  the  same 
general  area  as  the  original  randomly  selected 
locations ,  to  allow  for  a  better  measure  of 
trends  due  to  time  alone ,  and 

(d)  keep  a  separate  series  for  the  upper  and  lower 
course  materials.  The  difference  is  minimal 
at  the  time  of  construction,  but  is  of  con¬ 
siderable  importance  as  the  surface  course 

is  subjected  to  the  weathering  and  aging 


elements . 
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RECOMMENDATIONS  FOR  FURTHER  STUDY 


1.  The  tensile  splitting  test  should  be  altered  to 
determine  the  effect  of  different  loading  rates, 
and  to  establish  values  for  Poisson's  ratio  for 
different  temperatures  and  different  loading  rates . 

2.  Investigate  the  relationship  between  the  low  temp¬ 
erature  tensile  properties,  as  measured  by  the 
tensile  splitting  test,  and  the  recovered  pro¬ 
perties  of  asphalt  cement  for  various  service 
lives . 

3.  Determine  the  value  of  the  additional  parameters 
computed  within  the  "Detailed  Computer  Analysis" 
program  (i.e.  ,  Stiffness  of  Bitumen,  Bitumen 
Strain,  Toughness,  and  Work  Input)  as  an  indicator 
of  cracking  susceptibility  for  asphaltic  concrete. 

4.  Investigate  the  possible  effects  of  temperature 
at  coring  on  the  low  temperature  tensile  proper¬ 
ties  of  asphalt  concrete  cylinders. 
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APPENDIX  A 


METHOD  OF  TEST  AND  ANALYSIS  FOR  THE 
LOW  TEMPERATURE  TENSILE  PROPERTIES  OF 
ASPHALT  CONCRETE  CYLINDERS  USING  THE 


TENSILE  SPLITTING  TEST 
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METHOD  OF  TEST  AND  ANALYSIS  FOR  THE 
LOW  TEMPERATURE  TENSILE  PROPERTIES  OF 
ASPHALT  CONCRETE  CYLINDERS  USING  THE 
TENSILE  SPLITTING  TEST 

1 .  Scope 

1.1  This  method  covers  the  procedure  developed  for 
determining  the  low  temperature  tensile  properties  of 
asphalt  concrete  cylinders  using  the  tensile  splitting  test. 
The  test  can  be  conducted  on  asphalt  concrete  laboratory 
specimens  and  cored  pavement  specimens  (Note  1) . 

NOTE  1  -  For  method  of  making  laboratory  speci¬ 
mens  see  The  Standard  Method  of  Test  for  Resistance  to 
Plastic  Flow  of  Bituminous  Mixtures  Using  Marshall  Apparatus 
(ASTM  Designation:  D  1559-65).  Alternate  methods  for  pre¬ 
paration  of  laboratory  specimens  may  be  used. 

2 .  Summary  of  Method 

2.1  The  tensile  splitting  test  method  consists  of 
loading  an  asphalt  concrete  cylinder  via  loading  strips 
across  a  diameter,  in  a  compression  testing  frame  and  within 
a  controlled  temperature  chamber  maintained  at  a  constant 
low  temperature.  Output  signals  from  a  load  cell  and  two 
series  connected  linear  variable  differential  transformers 
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(attached  to  opposite  ends  of  the  specimen)  are  monitored 
on  a  two  channel  recorder.  From  the  load  trace  and  the 
deformation  trace  a  limited  analysis  of  the  test  can  result 
in  failure  stress,  failure  strain,  and  a  manually  plotted 
stres s -strain  diagram.  The  limited  analysis  is  covered  in 
sections  up  to  and  including  Sections  7.1.5. 

2.2  More  detailed  analyses  involving  the  use  of  an 

IBM  360-67  Computer  and  a  Model  770/663  Calcomp  Plotter  at 
the  University  of  Alberta  Computing  Center  are  also  used 
and  are  described  in  detail  in  Sections  7.2.  to  8.5.1.  The 
computer  program  summarizes  the  physical  properties  of  the 
specimen  as  well  as  various  computed  parameters  at  several 
points  throughout  the  duration  of  the  test.  At  the  comple¬ 
tion  of  a  particular  test  series,  a  statistical  summary  giving 
the  mean,  standard  deviation,  and  coefficient  of  variation 
of  the  parameters  considered  is  provided  for  a  more  meaning¬ 
ful  analysis.  In  addition,  with  the  use  of  the  Calcomp 
Plotter  automatically  plotted  stress -strain  and  stiffness- 
strain  plots  are  developed  for  the  specimens  within  a  test 
series . 

3  .  Significance 

3.1  This  method  determines  the  tensile  stress -strain 

and  stiffness-strain  characteristics  of  asphalt  concrete  at 

low  temperatures  and  is  primarily  intended  to  assist  in  the 

design  and  evaluation  of  asphalt  concrete  with  respect  to 
th.erjmal  cracking. 
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4 .  Apparatus 
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4.1  Controlled  Temperature  Chamber  -  The  controlled 
temperature  chamber  shall  be  capable  of  maintaining  test 
specimens  at  a  constant  temperature  -2°  within  the  range  of 
+40°  to  -10°F  during  the  course  of  a  test.  A  temperature 
monitoring  device  shall  have  its  sensor  embedded  in  a  speci¬ 
men  of  similar  size  and  composition  to  the  specimen  which 

is  to  be  tested  and  shall  be  capable  of  measuring  temperature 
to  ±1°F. 

4.2  Loading  apparatus  shall  consist  of  the  following: 

4.2.1.  Compression  Testing  Frame'*'  -  The  Compression 
frame  shall  have  a  minimum  capacity  of  5  tons  and 
shall  be  capable  of  providing  the  rate  of  loading 
prescribed  in  Section  6.4.1. 

4.2.2.  Supplementary  Bearing  Bar  or  Plate  -  The 
supplementary  bearing  bar  or  plate  shall  conform 

to  the  specifications  for  this  item  in  the  Standard 
Method  of  Test  for  Splitting  Tensile  Strength  of 
Molded  Concrete  Cylinders  (ASTM  Designation:  C  496- 
66),  except  that  the  width  of  the  bearing  bar  or 
plate  shall  be  not  less  than  1.3  inches. 

4.2.3.  Bearing  Strips  -  The  bearing  strips  shall  con 

form  to  the  specifications  for  this  item  in  Method  C 

496-66,  except  that  nominal  3/16  inch  plywood, 

approximately  1/4  inches  wide  shall  be  used. 

1.  A  suitable  device  (Wykeham  Farrance  Mod.  57,  5  ton  com¬ 
pression  tester)  may  be  obtained  from  Wykeham  Farrance  Eng 
ingeering,  Ltd.,  127  Edinburgh  Avenue,  Slough,  Bucks,  U.K. 


' 
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4.3  Load  Measurement  Apparatus  shall  consist  of  the 

following : 

4.3.1.  Load  Cell^  -  The  load  cell  shall  have  a 
minimum  capacity  of  5  tons  and  shall  be  capable 
of  measuring  compressive  loading  to  -1  per  cent 
of  true  at  the  rate  of  loading  prescribed  in 
Section  6.4.1. 

Gauge  Points,  and  Marking  and  Mounting  Apparatus 
consist  of  the  following: 

4.4.1.  Gauge  Points  -  The ' gauge  points  shall  be 
3/8  x  3/8  x  3/16  inches  (-  0.001  inches  from  mean 
in  any  dimension)  brass  plates. 

4.4.2.  Gauge  Point  Jig  -  The  gauge  point  jig  shall 
provide  slots  for  marking  the  specimen  and  holes  for 
mounting  the  Gauge  Points  (Fig.  A1  of  Andepson  S 
Hahn,  1968). 

4.5  Deformation  Measurement  Apparatus  shall  consist  of 

the  following: 

3 

4.5.1.  Displacement  Gauges  -  The  displacement 
gauges  shall  be  two  series  connected  linear  variable 
differential  transformers  of  matched  sensitivity 
(within  5%)  and  be  capable  of  measuring  displace- 
ments  to  within  -0.00005  inches,  and  shall  have  a 

stroke  of  not  less  than  -0.010  inches. 

2.  A  suitable  device  (Kwoya  Musen  Load  Cell  Mod.  LC-5,5  ton) 
may  be  obtained  from  Kwoya  Musen  Kenkyujo  Co.,  Ltd.,  Tokyo, 
Japan . 

3.  Suitable  devices  Sanborn  Linear  Variable  Differential  Trans¬ 
formers  Mod.  595  DT  025)  may  be  obtained  from  the  Sanborn  Co. 
175  Wyman  Street,  Waltham  54,  Massachusetts,  U.S.A. 


4.4 

shall 


.  . 
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4.5.2  Displacement  Gauge  Core  and  Coil  Assemblies  - 
The  two  displacement  gauge  core  and  coil  assemblies 
which  hold  the  Displacement  Gauges  shall  be  made  of 
brass  (Fig.  A2  of  Anderson  S  Hahn,  1968). 

4.5.3  Displacement  Gauge  Calibration  Jig  -  The 
displacement  gauge  calibration  jig  shall  be  made  of 
brass  and  aluminum  (Fig.  A3  of  Anderson  6  Hahn,  1968). 
The  dial  gauge  which  comprises  a  portion  of  the 
displacement  gauge  calibration  jig  shall  be  a  0.00001 
inch  dial  gauge. 

4.6  Recorder14  -  The  Recorder  shall  be  of  a  type  having 

two  channels  capable  of  providing  a  load  trace  on  one  channel, 
from  the  output  signal  generated  by  the  Load  Cell,  and  a  de¬ 
formation  trace  on  the  second  channel,  from  the  output  signal 
generated  by  the  Displacement  Gauges. 

5  .  Test  Specimens 

5.1  Asphalt  Concrete  Laboratory  Specimens  -  If  Marshall 
specimens  are  to  be  tested  they  shall  conform  to  the  specifi¬ 
cations  set  forth  in  ASTM  Method  D  1559-65. 

5.2  Asphalt  Concrete  Cored  Pavement  Specimens  -  If 
cored  pavement  specimens  are  to  be  tested  they  shall  be  trimmed 
to  a  cylindrical  shape  (within  -0.01  inches  of  the  mean  length 
and  diameter)  having  a  diameter  of  4  inches,  -  0.10  inches 

and  a  length  of  less  than  4  inches. 

4i  A  suitable  device  (Sanborn  Recorder  Model  311)  may  be 
obtained  from  the  Sanborn  Company. 
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6 .  Procedure 

6.1  Calibration  consists  of  the  following: 

6.1.1.  Load  Cell  -  The  Load  Cell  shall  be 
calibrated  at  room  temperature  (if  temperature 
compensating)  or  at  the  test  temperature  (if  non¬ 
temperature  compensating),  on  a  Compression  Tester 
whose  load  accuracy  has  been  verified  to  -1  per  cent 
in  accordance  with  the  Standard  Methods  of  Verifi¬ 
cation  of  Testing  Machines,  ASTM  Designation: 

E4-64.  It  is  recommended  that  the  load  be  cali¬ 
brated  on  Channel  B  of  the  Recorder. 

6.1.2.  Dial  Gauge  -  The  dial  gauge  shall  be  cali¬ 
brated  while  on  the  Displacement  Gauge  Calibration 
Jig  described  in  Section  4.5.3,  using  machinist's 
gauge  blocks. 

6.1.3.  Displacement  Gagues  -  The  two  most  closely 
matched  Displacement  Gauges  shall  be  chosen  on  the 
basis  of  separate  calibration  of  each  available 
Displacement  Gauge  at  room  temperature  on  the  Dis¬ 
placement  Gauge  Calibration  Jig.  Then  the  two  Dis¬ 
placement  Gauges  shall  be  series  connected  and  cal¬ 
ibrated  on  the  Displacement  Gauge  Calibration  Jig 
(using  a  1.000  inch  gauge  length  at  null)  at  the 
test  temperature.  It  is  recommended  that  the  de¬ 


formation  be  calibrated  on  Channel  A  of  the  Recorder. 


- 
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6.2  Preparation  of  Specimen  for  Testing  is  as  follows: 

6.2.1.  Measurement  -  Determine  the  length  and  dia¬ 
meter  of  the  test  specimen  to  the  nearest  0.01  inch 
by  averaging  four  readings  at  each  dimension. 

6.2.2.  Marking  -  Mark  dimetral  loading  points  on 
each  end  of  the  specimen  in  the  same  axial  plane 
using  the  Gauge  Point  Jig  described  in  Section  4.4.2. 

6.2.3.  Gauge  Point  Attachment  -  Coat  one  side  of 

each  of  two  gauge  points  (described  in  Section  4.4.1.) 
with  warm  asphalt  cement.  Insert  the  two  coated 
Gauge  Points  through  the  holes  in  the  aligned  Gauge 
Point  Jig  and  press  firmly  onto  the  specimen. 

Store  the  specimen  horizontally  in  a  cold  atmosphere 
for  approximately  15  minutes  to  firmly  affix  the 
gauge  points  to  the  specimen.  Remove  the  specimen 
from  the  cold  atmosphere ,  invert  the  specimen 
(and  prop  in  a  manner  that  will  not  disturb  the  pre¬ 
viously  attached  Gauge  Points)  and  attach  the  other 
two  Gauge  Points  in  a  similar  manner. 

6.2.4.  Cooling  -  Immediately  place  the  specimen 
into  the  Controlled  Temperature  Chamber. 

6.3.  Preparation  for  Loading  of  the  Specimen  is  as 

follows : 

6.3.1.  Specimen  Inspection  -  After  the  specimen  to 
be  tested  has  reached  equilibrium  temperature,  in¬ 
spect  it  for  Gauge  Point  Slippage.  If  any  slippage 

is  evident  remove  the  Gauge  Points  and  repeat  steps 
6.2.3  and  6.2.4. 


. 
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6.3.2.  Positioning  -  Place  the  Load  Cell  on  the 
loading  ram  platen  of  the  Compression  Testing 
Frame.  Position  the  specimen  so  that  the  marked 
loading  points  are  in  a  vertical  plane  passing 
through  the  center  of  thrust  and  so  that  the  long¬ 
itudinal  axes  of  the  plywood  Bearing  Strips  are  in 
this  vertical  plane.  Raise  the  loading  ram  of  the 
Compression  Testing  Frame  just  enough  to  secure 
the  specimen  for  Displacement  Gauge  attachment. 

6.3.3.  Displacement  Gauge  Attachment  -  Tie  both 

of  the  Displacement  Gauge  Core  and  Coil  Assemblies 
to  some  point  on  the  Compression  Testing  Frame  to 
obviate  damage  after  specimen  failure.  Simultan¬ 
eously  place  the  rear  Displacement  Gauge  Core  and 
Coil  Assembly  onto  the  Gauge  Points  and  then  secure 
the  assemblies  by  tightening  the  alien  screws.  Re¬ 
peat  the  foregoing  attachment  procedure  for  the 
front  Displacement  Gauge  Core  and  Coil  Assembly. 

6.4.  Loading  and  Recording  procedure  is  as  follows: 

6.4.1.  Loading  Rate  -  Set  the  Compression  Testing 
Frame  to  a  nominal  loading  rate  of  0.06  in/min. 

The  actual  loading  rate  may  vary  from  the  nominal 
loading  rate  by  -10  per  cent  but  must  be  reproduc¬ 
ible  within  -1  per  cent. 
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6.4.2.  Loading  -  Engage  the  Compression  Testing 
Frame  and  return  to  the  recording  area  (Note: 
that  loading  will  not  begin  until  the  power 
supply  switch  for  the  servo  motor  is  closed.  This 
switch  should  be  located  in  the  recording  area, 
adjacent  to  the  recorder). 

6.4.3.  Recording  -  Balance  the  resistance  and  cap¬ 
acitance  of  the  Displacement  Gauge  Core  and  Coil 
Assemblies  using  the  procedure  indicated  in  the 
recorder  manual.  Engage  the  chart  movement  switch 
to  establish  a  chart  movement  rate  of  1  millimeter 
per  second  simultaneously  with  the  closing  of  the 
power  supply  switch  for  the  servo  motor.  Set  the 
attenuation  for  each  of  channels  A  and  B  at  the 
most  sensitive  calibrated  scales  that  will  be 

used.  Align  the  recording  pens  for  each  of  channels 
A  and  B  along  a  line  5  millimeters  in  from  the 
right  hand  markings  on  the  strip  chart.  As  the 
recorded  load  and  deformation  increase,  the  atten¬ 
uation  should  be  adjusted  so  that  all  recording 
takes  place  within  the  central  40  millimeters  of 
the  50  millimeter  tape  width.  All  attenuation 
changes  shall  be  indicated  on  the  strip  chart  to 
facilitate  establishment  of  proper  scales.  All 
other  pertinent  data  such  as  date  of  test,  test 
temperature,  rate  of  loading,  and  specimen  identi¬ 
fication  number  shall  be  recorded  on  the  strip 
chart  at  this  time. 
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6*^.4.  Termination  of  Test  -  Upon  failure  of  the 
specimen  disengage  the  chart  movement  switch,  turn 
the  attenuator  for  channel  A  to  the  off  position, 
and  disengage  the  power  supply  switch  for  the 
servo  motor.  Mark  the  failure  point  on  the  load 
trace  on  channel  B  (at  the  point  where  the  slope 
first  reaches  zero).  Project  the  corresponding  point 
onto  channel  A  to  establish  the  failure  strain. 
Disengage  the  Compression  Testing  Frame  and  ex¬ 
amine  the  fractured  specimen.  If  the  fracture 
surface  passes  under  a  Gauge  Point  the  test  shall 
be  rejected. 


7 .  Calculations 

7.1  Limited  Analysis  -  The  limited  analysis  may  be  con¬ 

ducted  with  simple  calculations  not  requiring  the  use  of  a 
computer.  The  results  obtainable  with  this  method  include 
failure  stress,  failure  strain,  and  the  coordinates  of  the 
stress-strain  diagram. 


7.1.1  Tensile  Stress  -  The  tensile  stress  at  any 
point  to  failure  shall  be  calculated  as  follows : 


T 


2P 
it  td 


where : 

T  =  tensile  stress,  in  pounds  per  square  inch 

P  =  applied  load  as  indicated  on  channel  B, 
in  pounds 


, 


* 


All 


t  =  specimen  thickness  in  inches  and 
d  =  specimen  diameter,  in  inches. 

7.1.2.  Strain  -  The  strain,  at  any  point  to  failure 
is  equal  to  the  deformation,  in  inches,  indicated 
on  channel  A  of  the  Recorder  (see  Note  2). 

NOTE  2  -  Due  to  the  biaxial  state  of  stress 
existing  within  the  cylindrical  specimen,  the 
displacement  measured  between  the  gauge  points  is 
a  result  of  both  compressive  stresses  in  the  vertical 
direction  and  tensile  stresses  in  the  horizontal 
direction.  The  term  strain  is  used  without  differ¬ 
entiation  as  to  its  cause.  If  tensile  strain  is 
desired,  as  for  calculation  of  a  stiffness  modulus, 
use  of  equations  applying  the  Generalized  Hooke's 
Law  is  necessary. 


7 . 1 .  3.  Failure  Strain^  -  The  failure  strain  shall  be 
considered  as  the  strain  corresponding  to  the  point 
on  channel  A  of  the  strip  chart  located  by  the  pro¬ 
jection  of  the  point  where  the  slope  of  the  load 
trace  on  channel  B  first  reaches  zero. 


7.1.4.  Tensile  Strength  -  The  tensile  strength  shall 
be  considered  as  the  tensile  stress  calculated  using 
the  load  corresponding  to  the  point  of  failure 
(see  7.1.3.). 


. 
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7.1.5.  Report  -  The  report  for  the  Limited 
Analysis  shall  include  the  following: 

(1)  Identification  number,  aggregate  identi¬ 
fication,  asphalt  cement  penetration  or 
viscosity,  and  asphalt  cement  supply, 

(2)  Test  temperature, 

(3)  Rate  of  loading, 

(4)  Specimen  diameter,  and  thickness, 

(5)  The  load  and  deformation  strip  chart, 

(6)  The  failure  strain, 

(7)  The  tensile  strength,  and 

(8)  Any  abnormalities  in  the  type  of  fracture. 


7.2.  Computer  Analyses 
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7.2.1.  The  computer  analyses  consists  of  five 
separate  programs  indicated  as  follows: 

1)  Basic  computer  analysis  (section  7.3.) 

2)  Detailed  computer  analysis  (section  7.4.) 

3)  Computer  analysis  with  stress-strain  plot 
(section  7.5.) 

4)  Computer  analysis  with  stiffness-strain 
plot  (section  7.6.) 

5)  Computer  analysis  with  stress-strain  and 
stiffness-strain  plots ( section  7.7.) 

The  programs  indicated  are  progressively  more 
complex  in  function  and  in  the  required  estab¬ 
lishment  of  the  data  decks.  The  relationship 
of  one  program  to  another,  with  respect  to  data 
cards ,  is  best  understood  by  making  reference 
to  the  flow  chart,  Figure  Al.  The  remaining 
topics  covered  in  Section  7.2.  indicate  the 
procedure  to  be  followed  for  obtaining  infor¬ 
mation  common  to  all  computer  programs. 

7.2.2.  Determine  the  number  of  sets  of  specimens 
to  be  processed  during  this  run. 

7.2.3.  Establish  a  descriptive  title  (80  columns 
or  less)  that  is  to  appear  at  the  top  of  each 
page  to  identify  each  particular  set.  Generally 
it  would  include  details  pertaining  to  asphalt 
cement  supply,  grade,  aggregate  source,  and  the 


‘ 


■ 
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FLOW  CHART  FOR  ARRANGING  DATA 
FOR  ALL 

COMPUTER  PROGRAMS 


FIGURE  A1 


Flow  Chart  For  Arranging  Data  For  All 
Computer  Programs 
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highway  or  test  section  description.  It  is  not 
necessary  to  include  specimen  numbers ,  test 
temperature ,  loading  rate ,  or  date  of  test  for 
a-11  of  this  is  included  in  a  different  format 
within  the  remainder  of  the  computer  analysis. 

1.2.4.  Note  the  test  temperature  in  degrees 
Fahrenheit  ,  the  rate  of  loading  in  inches  per 
minute,  and  the  test  date  for  each  of  the  sets 
to  be  processed. 

7.2.5.  Determine  the  number  of  specimens  tested 
within  each  series. 

7.2.6.  Summarize  previously  obtained  information 
as  it  applies  to  the  individual  specimens  tested. 
This  will  include  the  following: 

(1)  specimen  identification  number, 

(2)  weight  in  air,  in  grams, 

(3)  weight  in  water,  in  grams, 

(4)  diameter,  in  inches, 

(5)  thickness,  in  inches, 

(6)  the  number  of  load  and  strain  points 
picked  off  the  load  and  deformation  strip 
chart  (see  Section  7.2.7.),  and 

(7)  the  time  to  failure,  in  seconds. 

7.2.7.  Several  time,  load,  and  deformation 
values  are  picked  off  the  load  and  deformation 
strip  chart,  to  provide  the  information  required 
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"to  allow  the  computer  analysis  of  various  para¬ 
meters  throughout  the  loading  of  the  specimen, 
up  to  the  failure  point.  During  the  initial 
loading  period,  it  should  be  realized  that  a 
certain  amount  of  time  and  energy  are  required 
to  compress  the  plywood  loading  strips ,  and  that 
the  load  transferred  to  the  specimen  during  this 
period  would  be  of  a  very  small  order.  This  fact 
makes  it  necessary  to  establish  the  zero  time  by 
some  criteria  that  would  be  indicative  of  the 
time  at  which  the  plywood  strips  are  completely 
compressed,  from  this  point  forward  the  load  is 
being  carried  entirely  by  the  specimen.  This 
zero  time  point,  necessary  for  calculating  the 
work  input  (Section  7.4.),  is  taken  arbitrarily 
as  the  point  at  which  the  load  trace  shows  a 
marked  increase  in  slope  after  having  held  re¬ 
latively  constant  with  a  load  in  the  order  of 
300  to  600  pounds  (this  will  vary  with  the  type 
of  plywood) .  To  establish  a  sufficient  number 
of  points  throughout  the  testing  of  the  specimen, 
the  time  in  seconds,  load  in  pounds,  and  defor¬ 
mation  in  ten- thousandths  of  an  inch  should  be 
recorded  at  the  following  deformation  values  up 
to  the  failure  point  (See  Section  7.1.3.)  (for 
simplicity  all  values  of  deformation  are  express¬ 
ed  as  the  number  of  ten-thousandths  of  an  inch: 


. 
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1  ,  2,  3  ,  4  ,  5  ,  10  ,  15  ,  20  ,  30  ,  40  ,  etc.  until 
failure ) . 

7.2.8.  All  of  the  information  summarized  as  out¬ 
lined  in  Sections  7.2.2.  -  7.2.7.  shall  be  key¬ 
punched  in  the  proper  fields  as  indicated  in  the 
program  listings.  Since  the  same  data  cards  will 
eventually  be  used  in  conjunction  with  the  plot¬ 
ting  programs ,  special  care  shall  be  exercised 

in  establishing  the  title  card.  In  the  plotting 
programs  the  80  column  title  is  split  into  two 
40  column  titles  and  placed  one  above  the  other 
on  the  Calcomp  Plotter  output,  therefore,  the 
spacing  of  the  title  on  the  data  card  shall  be 
properly  spaced  within  columns  1  to  40,  and 
columns  41  to  80 ,  to  develop  a  neatly  spaced 
title  on  the  computer  plotted  graphs. 

7.2.9.  Upon  completion  of  all  key-punching  the 
data  cards  shall  be  arranged  in  the  following 


order : 
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7.3.  Basic  Computer  Analysis 

7.3.1.  The  required  job  control  cards,  the  "Basic 
Computer  Analysis"  program  (in  precompiled  object 
deck  form) ,  and  the  data  deck  (as  arranged  in 
Section  7.2.9.)  shall  be  arranged  as  follows: 

//  jobname  JOB  (user  chg  no  ,  2  ,  2 )  , ' BASI C  COMPUTER 

ANALYSIS ' 


//  stepname 


EXEC 


FORTHLG 
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//  LKED .  SYS IN  DD  * 

BASIC  COMPUTER  ANALYSIS  PROGRAM 
(IN  PRECOMPILED  OBJECT  DECK  FORM) 

/* 

//  GO.  SYSIN  DD  * 

DATA  DECK  (SEE  SECTION  7.2.9.) 

/* 

7.3.2.  The  computer  listing  of  the  "Basic  Comp¬ 
uter  Analysis"  program  is  referred  to  in  Section 
8.1.1.  The  output  from  the  "Basic  Computer 
Analysis"  program  shall  consist  of  one  sheet  for 
each  of  the  specimens  tested,  plus  two  sheets 
providing  the  statistical  summaries  for  this 
particular  set.  Section  8.1.2  refers  to  an  out¬ 
put  typical  of  that  which  shall  be  provided  for 
each  of  the  specimens  tested  within  this  test 
series.  When  the  last  specimen  within  a  series 
has  been  completed,  statistical  summaries  shall 
be  provided.  The  first  statistical  summary 
(Section  8.1.3.)  is  based  upon  initial  conditions 
within  the  tensile  splitting  test.  This  point 
within  each  test  corresponds  to  the  maximum 
secant  modulus  of  the  stress-strain  diagram  (i.e. 
the  maximum  stiffness  of  mix  for  each  particular 
specimen  within  this  test  series).  The  second 
statistical  summary  (Section  8.1.4.)  is  based 


■ 
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upon  failure  conditions  within  the  tensile 
splitting  test.  The  failure  conditions  within 
each  specimen  tested  shall  be  considered  as  the 
point  at  which  the  failure  strain  (Section  7.1.3.) 
has  been  attained. 

7.4  Detailed  Computer  Analysis 

7.4.1.  The  data  deck  for  the  Detailed  Computer 
Analysis  shall  consist  of  the  data  deck  provided 
in  Section  7.2.9.  plus  one  additional  card  for 
each  test  series  within  this  run.  The  one  addit¬ 
ional  card  for  each  test  series  shall  provide  the 
following  additional  information: 

1)  asphalt  content  of  the  mix  being  tested, 
in  pounds  of  asphalt  cement  per  100  pounds 
of  dry  aggregate, 

2)  bulk  specific  gravity  of  the  aggregate, 

3)  bulk  specific  gravity  of  the  asphalt  cement, 
and , 

4)  per  cent  absorption  of  the  aggregate,  in 
pounds  of  asphalt  cement  absorbed  per  100 
pounds  of  dry  aggregate. 

7.4.2.  The  information  summarized  in  Section  7.4.1. 
shall  be  key-punched  in  the  proper  fields  as 
indicated  in  the  program  listing. 

7.4.3.  The  data  card  for  each  test  series  as 
described  in  Sections  7.4.1.  and  7.4.2.  shall  be 
added  to  the  data  deck  described  in  Section  7.2.9. 


' 
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immediately  following  the  "Temperature,  Rate, 
and  Date"  card,  and  just  prior  to  the  "N"  card. 

The  data  deck  for  the  "Detailed  Computer  Analysis" 
shall  be  arranged  in  the  following  order: 

THESE  4  CARDS  ARE 


7.4.4.  The  required  job  control  cards,  the 
"Detailed  Computer  Analysis"  program  (in  pre¬ 
compiled  object  deck  form) ,  and  the  data  deck 
(as  arranged  in  Section  7.4.3.)  shall  be  arranged 


. 
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as  follows : 

//  jobname  JOB  (user  chg  no,  2,2),  ’DETAILED 

COMPUTER  ANALYSIS’ 

//  stepname  EXEC  FORTHLG 

//  LKED .  SYSIN  DD  * 

DETAILED  COMPUTER  ANALYSIS  PROGRAM 
(IN  PRECOMPILED  OBJECT  DECK  FORM) 

/  * 

//GO.  SYSIN  DD  * 

DATA  DECK  (SEE  SECTION  7.4.3.) 

7.4.5.  The  computer  listing  for  the  "Detailed 
Computer  Analysis"  program  is  referred  to  in 
Section  8.2.1.  The  output  from  the  "Detailed 
Computer  Analysis"  program  shall  consist  of  one 
sheet  for  each  of  the  specimens  tested,  plus  two 
sheets  providing  the  statistical  summaries  for 
this  particular  set.  Section  8.2.2.  refers  to 
an  output  typical  of  that  which  shall  be  provided 
for  each  of  the  specimens  tested  within  this  test 
series.  Section  8.2.3.  refers  to  an  output 
typical  of  that  which  shall  be  given  for  the 
statistical  summary  based  upon  initial  conditions, 
while  Section  8.2.4.  similarity  refers  to  the 
statistical  summary  based  upon  final  conditions. 
The  terms  initial  and  final  as  used  within  this 
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program  and  this  section  shall  be  as  defined 
in  Section  7.3.2. 

7.5.  Computer  Analysis  with  Stress-Strain  Plot 

7.5.1.  Upon  inspection  of  test  results  within 
each  test  series,  the  maximum  failure  strain,  in 
ten  thousandths  of  an  inch/inch,  attained  by  any 
specimen  within  a  series  shall  be  noted.  This 
maximum  failure  strain  value  shall  be  used  to 
establish  the  appropriate  scale  for  the  tensile 
strain,  in  inches /inch,  as  measured  along  the 
X-axis  of  the  stress-strain  plot.  The  maximum 
failure  strain  may  be  any  value,  but  since  the 
X-axis  on  the  stress-strain  plot  is  five  inches 
in  length  the  most  convenient  scales  will  be 
obtained  if  the  number  chosen  is  a  multiple  of 
five.  The  following  values  for  the  variable 
"MAXSTR"  produce  convenient  scales: 


MAXSTR 

SCALE 

10 

1"  = 

2 

ten 

thousandths 

of 

an 

in/ in 

25 

1"  = 

5 

ten 

thousandths 

of 

an 

in/in 

50 

1"  = 

10 

ten 

thousandths 

of 

an 

in/ in 

100 

1"  = 

20 

ten 

thousandths 

of 

an 

in/in 

7.5.2.  The  value  corresponding  to  the  maximum 
failure  strain  for  each  test  series  as  described 
in  Section  7.5.1.  shall  be  key-punched  in  the 
proper  field  as  indicated  in  the  program  listing. 

7.5.3.  The  data  card  for  each  test  series  as 


, 
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described  in  Section  7.5.1.  and  7.5.2.  shall  be 
added  to  the  data  deck  described  in  Section  7.2.9. 
immediately  following  the  "temperature,  rate, 
and  date"  card,  and  just  prior  to  the  "N"  card. 

The  data  deck  for  the  "Computer  Analysis  with 
Stress-Strain  Plot"  shall  be  arranged  in  the 
following  order: 


THESE  4  CARDS  ARE 


•aaSRSUKROEC 
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7.5.4.  The  plotting  operation  requested  by  this 
type  of  computer  program  is  an  off-line  operation 
that  necessitates  the  intermediate  storage  of 
data  to  be  plotted  on  a  seven  track  magnetic 
tape.  In  order  for  this  program  to  run  success¬ 
fully,  one  must  have  access  to  a  seven  track  tape 
mounted  by  the  computer  operator  prior  to  the 
execution  of  the  program.  This  can  be  accomplished 
by  the  following: 

1)  add,  CLASS  =  B  to  the  first  job  control 
card,  (See  Section  7.5.5.). 

2)  add  two  job  control  cards  to  indicate  the 
plot  parameters  and  the  code  number  of 
the  seven  track  plotting  tape.  These  two 
cards  shall  be  placed  after  the  /*  card 
following  the  object  deck,  and  prior  to 
the  //GO.  SYSIN  DD  *  card  immediately 
preceding  the  data  deck  (See  Section  7.5.5.). 

3)  complete  tape  mounting  slip  and  plot 
request  slip  as  required  by  the  Computing 
Center  (See  Section  7.5.6.  and  7.5.7.). 

These  two  slips  must  accompany  the  program 
when  being  submitted  for  processing. 

7.5.5.  The  required  job  control  cards,  the 
"Computer  Analysis  with  Stress-Strain  Plot" 
program  (in  precompiled  object  deck  form),  and 
the  data  deck  (as  arranged  in  Section  7.5.3.). 


p 

■ 
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shall  be  arranged  as  follows: 

/  /  jobname  JOB  (user  chg  no,  2,2),  'COMP  ANALYSIS-STRESS 

STRAIN  PLOT',  CLASS  =  B 

/  /  stepname  EXEC  FORTHLG 

/  /  LKED .  SYS IN  DD  * 

COMPUTER  ANALYSIS  WITH  STRESS -STRAIN  PLOT 
(IN  PRECOMPILED  OBJECT  DECK  FORM) 

/  * 

/  /  GO . PLOTTAPE  DD  UNIT  =  SYS  7T  ,  DISP  =  (NEW  KEEP), 

DSNAME  =  PLOTTAPE, 

/  /  VOLUME  =  SER  =  tape  number,  LABEL  =  (1,SL),  DCB  =  (DEN  =  1) 
/  /  GO.  SYSIN  DD  * 

DATA  DECK  (SEE  SECTION  7.5.3.) 

/  A 
/ 

7.5.6.  The  tape  mounting  slip  shall  be  as  required  by  the 
University  of  Alberta  Computing  Centre. 

7.5.7.  The  plotter  request  slip  shall  be  as  required  by  the 
University  of  Alberta  Computing  Centre. 
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7.6  Computer  Analysis  with  Stiffness-Strain  Plot 

7.6.1  Upon  inspection  of  test  results  within 
each  test  series,  the  maximum  failure  strain,  in 
ten  thousandths  of  an  inch/inch,  attained  by  any 
specimen  within  a  series  shall  be  noted.  This 
maximum  failure  strain  value  shall  be  used  to 
establish  the  value  of  the  variable  "MAXSTR" ,  as 
indicated  in  Section  7.5.1.  The  appropriate 
value  for  the  variable  "MAXSTR"  shall  be  key¬ 
punched  as  indicated  in  Section  7.5.2.,  and  added 
to  the  data  deck  as  indicated  in  Section  7.5.3. 

7.6.2  Upon  inspection  of  test  results  within 
each  test  series,  the  maximum  stiffness  of  mix, 

in  pounds  per  square  inch,  attained  by  any  specimen 
within  a  series  shall  be  noted.  This  maximum 
stiffness  of  mix  value  shall  be  used  to  establish 
the  appropriate  scale  for  the  stiffness  of  mix, 
in  pounds /square  inch,  as  measured  along  the 
Y-axis  of  the  stiffness-strain  plot.  The  maximum 
stiffness  of  mix  may  be  any  value ,  but  since  the 
Y-axis  on  the  stiffness-strain  plot  is  seven 
inches  in  length  the  most  convenient  scales  will 
be  obtained  if  the  number  chosen  is  divisible  by 
seven.  The  following  values  for  the  variable 
"MAXSTF"  produce  convenient  scales: 


- 

- 
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MAXSTF 

SCALE 

700000 

1"  =  100,000  PSI 

1400000 

1"  =  200,000  PSI 

2100000 

1"  =  300,000  PSI 

2800000 

1"  =  400,000  PSI 

3500000 

1"  =  500,000  PSI 

Etc  . 

7000000 

1"  =  1,000,000  PSI 

7.6.3  The  value  corresponding  to  the  maximum  stiffness  of 
mix  for  each  test  series  as  described  in  Section  7.6.2.  shall 
be  keypunched  in  the  proper  field  as  indicated  in  the  pro¬ 
gram  listing. 

7.6.4  The  data  card  for  each  test  series  as  described  in 
Section  7.6.2.  and  Section  7.6.3.  shall  be  added  to  the  data 
deck  described  in  Section  7.5.3.  immediately  following  the 
"maximum  strain"  card,  and  just  prior  to  the  "N"  card.  The 
data  deck  for  the  "Computer  Analysis  with  Stiffness-Strain 
Plot"  shall  be  arranged  in  the  following  order: 


»  ; 


■ 


7.6.5.  The  off-line  plotting  operation  for  the 
"Computer  Analysis  with  Stiffness-Strain  Plot" 
program  shall  he  the  same  as  that  described  in 
Section  7.5.4.  Section  7.5.5.  also  applies  to 
this  program  except  that  the  precompiled  object 
deck  shall  be  that  for  the  "Computer  Analysis 
with  Stiffness-Strain  Plot".  In  addition,  the 
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tape  mounting  slip  and  plotter  request  slip 
shall  be  the  same  as  that  described  in  Sections 
7.5.6.  and  7.5.7. 

7.7  Computer  Analysis  With  Stress-Strain  and  Stiffness- 
Strain  Plots 

7.7.1.  All  aspects  of  this  program  shall  be  the  same  as 
that  described  in  Section  7.6.  with  the  exception 
that  the  precompiled  object  deck  shall  be  that 
for  the  "Computer  Analysis  with  Stress-Strain 
and  Stiffness-Strain  Plots". 

8 .  Computer  Program  Listings  and  Sample  Output 

8.1.  Basic  Computer  Analysis 

8.1.1.  The  computer  listing  of  the  "Basic 
Computer  Analysis"  program  is  given  on  pages 
A3  3  to  A3  7 . 

8.1.2.  A  page  typical  of  that  giving  the  output 
for  each  specimen  within  a  test  series  for  the 
"Basic  Computer  Analysis"  program  is  given  on 
page  A38. 

8.1.3.  A  page  typical  of  that  giving  the  stat¬ 
istical  summary  based  upon  initial  conditions  (as 
described  in  Section  7.3.2.)  for  the  specimens  with¬ 
in  a  test  series  analyzed  with  the  "Basic  Computer 
Analysis"  program  is  given  on  page  A39 . 


A3 1 


8.1.4.  A  page  typical  of  that  giving  the  statist¬ 
ical  summary  based  upon  failure  conditions  (as 
described  in  Section  7.3.2.)  for  the  specimens 
within  a  test  series  analyzed  with  the  "Basic 
Computer  Analysis"  program  is  given  on  page  A40 . 

8.2.  Detailed  Computer  Analysis 

8.2.1.  The  computer  listing  of  the  "Detailed 
Computer  Analysis"  program  is  given  on  pages 
A41  to  A48. 

8.2.2.  A  page  typical  of  that  giving  the  out¬ 
put  for  each  specimen  within  a  test  series  for 
the  "Detailed  Computer  Analysis"  program  is 
given  on  page  A49 . 

8.2.3.  A  page  typical  of  that  giving  the  stat¬ 
istical  summary  based  upon  initial  conditions 

(as  described  in  Section  7.3.2.)  for  the  specimens 
within  a  test  series  analyzed  with  the  "Detailed 
Computer  Analysis"  program  is  given  on  page  A50. 

8.2.4.  A  page  typical  of  that  giving  the 
statistical  summary  based  upon  failure  conditions 
(as  described  to  Section  7.3.2.)  for  the  specimens 
within  a  test  series  analyzed  with  the  "Detailed 
Computer  Analysis"  program  is  given  on  page  A51. 


■ 
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8.3  Computer  Analysis  With  Stress -Strain  Plot 

8.3.1.  The  computer  listing  of  the  "Computer 
Analysis  with  Stress-Strain  Plot"  program  is 
given  on  pages  A52  to  A59. 

8.4  Computer  Analysis  with  Stiffness-Strain  Plot 

8.4.1.  The  computer  listing  of  the  "Computer 
Analysis  with  S tif fnes s -St rain  Plot"  program  is 
given  on  pages  A60  to  A6  7 . 

8 . 5  Computer  Analysis  with  Stress-Strain  and  Stiffness- 
Strain  Plots 

8.5.1.  The  computer  listing  of  the  "Computer 
Analysis  with  Stress-Strain  and  Stiffness-Strain 
Plots"  program  is  given  on  pages  A6 8  to  A7 7 . 


. 


- 


BASIC  COMPUTER  ANALYSIS 
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ISM  0002  DIMENSION  T I  II E ( 20 ) , D AT E (  5 )  ,  T I MC  I  20 ) , LO AD  I  2 0 )  ,  S TR A  I N I  2 2 ) , S TR E S S ( 2 2 

1),FLSTRS(2C)  ,FLSTRN(20)»FLTIME(20),SAMPLE(2C,3)  *  S  T  I-F  F  (  ?  0 )  ,  FLSTFF ( 2 

20)  ,  ORSTRSI  20  !  .  PRSTRN(?0  )  ,  ORSTFF  {  ?0  )  .  no.T  I  YFf  2D  .  !l\l  r  TWT(  20) 

ISM  0003  INTEGER  TI»E , SAMPLE , FA  I LTM, TEMP 

ISM  0004  REAL  LOAD 

C 

c 

C  M  =  THE  NUMBER  OF  SETS  OF  SAMPLES  TO  BE  PROCESSED  DURING 

C  THIS  R'JN  . 

C 

C 

ISN  0005  RE ADI  5 . 10)M 

I SN  0006  10  FORMAT  I  I 5) 

ISN  OC  07  DO  300  I  1=1, M 

C 

c 

C  TITLE  =  A  DESCRIPTION  I  80  COLUMNS  OR  LESS)  TO  APPEAR  AT  THE 

C  _  TOP  OF  EACH  PAGE  TO  IDENTIFY  R  A  CM  5  ARTICULAR  SET  . 

C 

C 

TSN  0008_ R  F  A  D  (  5  .  2  0  )  (TTTI  F  (  K  )  ,K=1  .20) 

ISN  0000  20  FORMAT ( 20A4) 

c 

c 

C  TEMP  =  THF  TEST  TEMPERATURE  ,  IN  DEGREES  FAHRENHEIT  ,  AT 

C  WHICH  THIS  SET  WAS  RUN  . 

c 

C  RATE  =  THE  RATE  OF  LOAD  APPLICATION  IN  INCHES  PER  MINUTE  . 

C 

C  DATE  =  THE  DATE  AT  WHICH  THE  SAMPLES  WITHIN  THIS  SET 

C  WERE  TESTED  . 

C 

c 

ISN  0010  READ! 5, 30) TEMP, RATE, (DATEIK) ,K=l, 5) 

ISN  0011  30  F0RMATI7X, 13  ,F10.3,5A4) 

r. 

c 

C  N  =  THE  NUMBER  OF  SAMPLES  WITHIN  THIS  SET  . 

£ 

c 

ISN  0012  RE  AD ( 5 , 40 ) N 

ISN  0013  40  FORMAT! 15) 

ISN  0014  DO  2C0  I  =  1 ,  M 

C 

c 

c  SAMPLE  =  THE  IDENTIFICATION  NUMBER  OF  THE  SAMPLE  . 

C 

C  V/TAIR  =  THE  WEIGHT  ,  IN  GRAMS  ,  OF  THE  SAMPLE  IN  AIR  . 

C 

C  WT WAT  =  THE  WEIGHT  ,  IN  GRAMS  ,  OF  THE  SAMPLE  IN  WATER  . 

C  ... 

C  DIAM  =  THE  AVERAGE  DIAMETER  OF  THE  SAMPLE  ,  IN  INCHES  . 

C 

c  THICK.  =  THE  AVERAGE  THICKNESS  OF  THF  SAMPLE  ,  IN  INCHES  . 

C 

C  L.  =  THF  NUMBER  OF  LOAD  AMD  STRAIN  POINTS  PICKED  OFF  THE 

C  IPAD  AND  DEFORMATION  STRIP  CHART  . 
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c 

C  FAILTM  =  THE  TIME  TO  FAILURE  ,  IN  SECONDS  . 

C 

_ c _ _ 

30  l*3  READ!  5,  50)  (  SAMPLE  (  I  ,K)  ,  K=  1  ,3  )  ,WTAIR,WTWAT,DIAM,THICK,l,FAllTM 

ISM  0016  5C  FORMAT.  (  3A4  ,  F  8.  3 , 3F  10 . 3 , 21  10) 

_ C _ , _ _ _ _ 

c 

C  SPGR  =  BULK  SR FC I F I C  GRAVITY  OF  THE  COMPACTED  MIXTURE  . 

_ C _ _ 

C 

.  ISM  0017  SP GP-WTA I R / ( WT A IR-WTWAT) 

_ C _ 

C 

c  UNITWT  =  UNIT  WEIGHT  OF  THE  COMPACTED  MIXTURE  , 

_ C _ IN  POUNDS  PER  CUBIC  FOOT  . _ _ 

C 

c 

_ ISM  0018 _ UNI  T  WT  (  I  )=S  PGR.  *  6  2 .4 _ _ _ 

ISM  0019  L M  =  M  + 1 

ISM  0020  NM=M+2 

ISM  3021 _ L'R  ITE  (  6 ,60  )  (  TITLE!  K  )  ,K=1 . 20)  ,  I  .  MM.  ISAMPI  F(I.K).K=1,3).SPGR.UN  I  T  WT  ( 

II)  ,01  AM, THICK 

ISN  0022  6  C  FORMAT!  1 H 1  ////// 20X  ,  20 A4/ /9 1 X ,  5HPAGE  ,I2,4H  OF  ,  I 2///44X, 17 H 

_ IS AM PL E  NUMBER  : ,  3A4//32X, _ 49HB ULK  SPECIFIC  GRAVITY  OF  THE  COM 

2PACTED  MIXTURE  =  ,  F5.3//38X, 1 4 HU N I T  WEIGHT  =  ,F5.l,22H  POUNDS 

2  PEP  CUBIC  FOOT  ,  / / 45 X , 1 1 HD  I  AM E TE R  =  ,  F5.3.7H  INCHES/ 

_ _ _ 3  /  4  4  X  ,  1  2  H  T  IT  I  C  K  M  E  S  S  =  ,  F5.3.7H  INCHES//) _ 

ISM  002?  WR I TE l 6 , 61 )  RAT E , FAI LTM, TEMP, ! D ATE ( K) ,K=1 , 5 ) 

ISM  0024  61  FORMAT ( 2 5X , 1 5HDETAI LS  OF  TE ST / / 2 6 X ,  1  8HR AT E  OF  LOADING  =  ,F5.3,16H 

_ _ _ 1  INCHES  /  MINUTE  .AX,  19HTI  M  E  TO  FRACTURE  =  .I3.8H  SEC0NDS//26X  ,  L9HTE 

2  ST  TEMPERATURE  :  ,I3,19H  DEGREES  FAHRENHE I T , 4X , 9HDAT E ,  :  ,5A4///) 


ISM  0025  WRITE! 6, 62) 

ISN  0026 _ 6  2  FORMAT !6CXt9HSTirFNESS/30X,6HSTPFSS.9X.6HSTRAIN.llX. 6H0F  M I  X.  I  OX. 


14HTIME/31X , 5  H ( P  S I )  ,9X, 

7H  ( 

I  N / I N)  , 10X 

,5H! PS  I )  , 9X 

, 9H( SECONDS )/// ) 

ISN 

0027 

C 

DEMCM=  3 . 14159*DIAM*THICK 

C 

C 

TIME  =  THE  TIME  , 

I  N 

SECONDS  , 

FROM  ZERO 

TIME  TO  EACH  OF  THE 

c 

POINTS  ON 

THE 

LOAD  AND 

DEFORMATION 

STR  IP  CHART  . 

ISN 

0028 

c 

c 

P.  E  A  D !  5 .70)  (  T  IME(K)  ,K=1 

ISN 

0029 

c 

c 

70 

FORMAT! 5X, 15 15) 

c 

LOAD  =  THE  LOAD  , 

IN 

POUNDS  , 

AT  EACH  OF 

THE  POINTS  ON  THE 

c 

c 

LOAD  TRACE 

• 

ISN 

0030 

c 

READ! 5 , 00) (LOAD(K) ,K=l 

» L  > 

ISM 

0031 

80 

FORMAT ( 5X, 15F5 .0) 

C 

c 

c  STRAIN  =  THE  STRAIN  .  IN  TFN  THOUSANDTHS  3=  AN  INCH  /  INCH  , 

C  ~  FOR  EACH  OF-  THE  POINTS  ON  THE  DEFORMATION  TRACE  . 

C 


. 

*  • 
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ISM  0032 
ISM  0033 
ISM  0034 

ISM  0035 
ISM  0036 


IS  M  ,  0  0  3  7 
ISM  0038 
ISM  0039 
ISM  0060 

ISM  0041 


ISN  0042 
ISM  0043 
ISM  0044 
ISM  0045 

ISM  0046 
ISM  0047 
ISN  0048 
ISM  0049 
ISM  0050 
ISN  0051 


I  S  M  0  0  5  2 
ISN  9053 
ISM  0054 
ISM  0055 

ISN  0056 
ISM  0057 
ISM  0058 
ISN  0059 
ISM  0060 


ISN  0061 
ISM'  0062 


c 

RE AQ( 5, 90 ) ( STRA IN(K )  ,K  =  1,L) 
90  FORMAT  (  5X,  15F5 .0) 

_ DP  100  .1  =  1  j  L _ 


STRESS!J)=2.*L0AD(J) /DEMON 
STRAIN(JJ=STRAIM(J)*0.0001 


C  CONSIDERING  AVERAGE  TENSILE  STRESS  OVER  A  ONE  INCH  GAUGE  LENGTH 

£_ _ AMD  ASSUMING  A  POISSONS  RATIO  OF  0.33  . 

C 

c 

_ ST  IT T( J)=0.912*STRGSS(J)/(STRAIN{J)*.5)^  1 _ 

I S  T I T  F  =  S  T ITT ( J) 

WR ITE ( 6 , 95 ) STRESS ( J ) ,STR A  INI J ) .  [STIFF, TIM E  C  J  >  ' 

_ 95  F  0  P  M  A  T ( 2  6  X  ,  F  10.2,6X,F10.6,3X.I3.9X.I5/) _ 

100  CONTINUE 
C  . 

C _ : _ _ _ . _ _  • _ 

C  SUBROUTINE  'MAX1  IS  CAl LED  TO  DETERMINE  THE  MAXIMUM  SECANT 

C  MODULUS  (  THAT  IS  THE  MAXIMUM  STIFFNESS  OF  MIX  FOR  EACH 

C _ STRESS  - STRAIN  DIAGRAM  )  . _ _ __ _ 

C 

c 

C  A  L JL _ OTA  X  CL  ,  S  T I  F  T  ,  N MNJ _ ; _ _ _ 

ORSTRS! I )  =  STRESS ( NNN ) 

ORSTRN (  !  )  =  STRA I  MI NNN ) 

_ ORSTEF (  I }  =  ST  IPF I NNN )  _ _ _ _ _ 

CRT  IMF (  I ) =  T I  ME ( NNN ) 

FLSTRSI I )= STRESS ( L ) 

_ EL  S  T P N J 1  )_=  S  T _R  A  IT;  (  L ) _ _ _ _ _ _ _ 

FLSTFFI  I  )  =  ST  I  FT ( L ) 

FLTIME ( I)=TI ME ( L ) 

200  CONTINUE _ _ _ _ _ _ _ _ 

C 

c 

_C _ SUBROUTINE  'STATS'  IS  CALL  ED  TO  DETERMINE  THE  MEAN  .STANDARD 

C  DEVIATION  ,  AND  COEFFICIENT  OF  VARIATION  OF  EACH  PARAMETER 

C  RASED  UPON  INITIAL  CONDITIONS  (  THAT  IS  THE  MAXIMUM  SECANT 

C  MODULUS  OF  EACH  STR F SS ~STR A I M  DIAGRAM  )  . _ 

C 

c 

_ C ALL  STATSIN, OPSTRS ,X RARE , S IGMA5.VAR5I _ _ _ 

CALL  STATS ( N  » ORSTRN , X0AR6 , S I GMA6 ,  VAR6 I 
CALL  STATS ( N  .ORSTFF ,XBAR7 , S IGMA7 ,VAR7 ) 

_ CAT  L  STATSIN, ORT  I  Mf  ,  XI3AR  8  ,  S  I  GMA8  ,  VAR8  ) _ 

CALI  STATS (M , UN  I  TUT , X B AR9 , S I G MA9 , VAR 9 ) 

I X  B  AP,  7-XRAR  7 

_ LS  GM  A  7  -  S  I  G_M  A  7 _ _ . _ 

WRITE! 6 ,210 )  ( T  I  TLE ( K) , K  =  1 ,20 )  ,LN, NN,N,  (  ( SAMPLE!  I  ,< )  ,K=1 ,3 ) , 1  =  1 , N) 

2  1C  FORMAT (  1H1//////20X,  20A4//91X,  5HPAGE  ,  I  2 , 4H  OF  ,I2///29X, 

_ 159HNEANS  ,  STANDARD  DEVIATIONS  .  AND  COEFFICIENTS  OF  VAR  I  AT  ION / 5 7 X 

2  »  3HFCR/44X  ,  4HTHE  ,12,2311  SAMPLES  IN  THIS  SERIES  / 4  4  X, 2  9H  BASED  UPON 
3INITIAL  COMDITIONS///15X, 17HSAMPLE  NUMBERS  :  ,3X,8!3A4)) 

WR  I  T  E  !  6 , 2 1  1  )  TEMP,  (  DATE  (  K  )  ,K=1,5) _ _ _ _ _ 

'211  FORMAT  I //26X ,19HTEST  TEMPERATURE  :  ,I3,19H  DEGREES  F AMR ENHE I T »  4X »  7 
1  HD  ATE  :  , 5A4) 


- 


. 


A36 


ISM  0063 
ISM  0064 


212 


WRITE!  6, 212) 

FCRMAT(////00X,9HSriFFNESS,21X,4HUNir/56X,SHSTRESS»6X,6HSTRAIN,8X, 
16HIT  MI  X, OX , 4HT I  ME , 9X, 6HWF IGHT/57X , 5H( PS [  ) , 6X , 7H(  l N/ IN ) , 7X , 5H( PS  I  ) 
2  t  7  X  ,  9  H  (  SECONDS)  »  5  X  ,  0  H  t  P.f.  .F.  )  ) _ 


ISN  0065 
ISM  0066 
J  SN  0067 
ISN  0068 

ISM  0069 


213 


214 


WR  I  TE  (6 ,21  3  )  XBAP.5  ,  XBARfc  ,  IXBAR7 , X  BARS , X3AX9 
FORMAT!/ /23X  ,4HMEAM,26X,F9.2,3X,F10.6,5X, I8,3X,F10.1,7X,F7.2/I 
WRIT:  ( 6, 214)  _s  IGMA5 , SIGMA6, I S0MA7 ,s IGMA3 , S I GMA9 

FORMAT ( 2  3  X , 1 8HSTANDARD  OEV I AT  ION, 12X , F 9.2 ,3X, FI  0. 6, 5X, I 8.3X.F 10 
17X,F7. 2/) 


ISN  0070 


215 


HP  I TE! 6, 215)  V AR 5 , V AR6 . V AR7 , V AR 8 , VA R9 
FORMAT (  23X ,  2  8HC0EFF I C I  ENT  OF  VARIATION 


_C_ 

c 

c 

c 


1 »3X,F10.2,7X  ,F7.2///) 


U),2X,F9.2,3X,F8.2»7X,F8.2 


SUBROUT INE 
C  EVIATIOM  i 


C 

c 

c 


'STATS'  IS  CALLED  TD  DETERMINE  THE  MEAN  ,  STANOARD 
—AND  COEFFICIENT  Q F  VARIATION  Q F  FACH  PARAMFTFR 
EASED  UPON  FAILURE  CONDITIONS  . 


I  SN 
ISN 
ISN 


0071 

0072 

0073 


CALL  STATS (N  ,FLSTRS,XBAR1 ,S IGMA1, VAR1 ) 
CALL  STATS (N  ,FL STRN , XBAR2 , SIGMA 2 , VAR2) 
CALL  ST  ATS ( N  ,  FL STFF , XBAR3 . S I GMA 3 . V AR 3 ) 


ISN 
ISN 
J  SN 
ISN 
ISN 


0074 

0075 

0n76 


0077 

0070 


2  50 


CAFL  STATS ( N  ,FLT I  ME ,XBAR4 , S IGMA4, VAR4 ) 

I  X  B  A  P  3  =  X  B  A  T  3 

I S  GM  A3  =  S  I G  M  A  3  _ _  _  _ 

WE  I  TF  (  6  ,25.0)  (  T  I  l  Ll  (  K  )  ,  K  =  1  ,70)  ,  NN,  NN  ,  N  ,  (  (SAMPLE!  I ,  K  ),  K=1  ,  3  )  ,  1  =  1,  N> 
FORMAT ( 1H1//////20X,  20A4//91X,  5HPACF  ,  I  2 , 4H  OF  .I2///29X, 

159UMFANS  ,  STANDARD  DEVIATIONS  ,  AND  COE  F  r  I  c  T  ENTS  OF  VARI  ATI.O  ’17.5  7  ' 

RIES/44X ,29HBASED  .JPOM 


I SN  0079 
ISN  ODBC 

ISN  0981 
ISN  0082 


ISN  0083 
ISN  0084 
ISN  0085 


2, 3HFCR/44X, 4HTHE  , I2,23H  SAMPLES  IN  THIS  S 
?F A  I  LURE  C0KDITI0NS///15X, 17HSAMPLE  NUMBERS 
_ WR  I  Tl-  (  6,251  )  TEMP,  (  DATE  ( K )  ,K=1  ,5  ) 


, 3X,8( 3A4)  > 


251  FORMAT ( / / 2 6 X  ,  I 9HTEST 
1 HO  ATE  :  ,5A4) 

WR I T  f ( 6,252) 


TEMPERATURE 


,  I  3 , 1 9  H  DEGREES  F AHR EN HE  I T , 4 X , 


252  FOR! 'AT!///  / OCX, 9 1!  STIFFNESS, 21X,4H  UN  IF/56X,6HSTRESS,6X, 6  H  STRAIN, 8X, 
16HPF  M I  X , 8  X , 4H  T I  ME , 9  X , 6  HW  E I GHT/ 5  7  X , 5  H ( P  S I ),6X,7H(  IN/IN)  ,7X»5H(PS.  „ 

_ 2, 7  X  ,  °H  (  SECONDS)  ,  5X ,  8H(  P  .  C  .  F  .  )  ) _ _ _ 

WR  ITS! 6,253 )  XBAR1 , XBAR2 ,  IXB AR3 , XBAR4 , XBAR9 
2  53  FORMAT ( // 2  3X , 4  H  M  E  A  N , 2  6  X , F  9 . 2 , 3  X ,F10.5,5X,I3,3X,F10.l,7X,F7.2/) 

WR  ITE! 6,2  54)  S I G  ’ 1 A l , S I G  M  A  2 , I SGMA3 , S I GMA4 , S I GMA9 _ 


ISN  0086 


ISN  0087 


2  54 


BUST  ANOARO  DEVIATIDN,12X,F9.2,3X,F10.6,5X,I8,3X,F10.1, 


ISN  0088 


ISN  0089 


ICO 


FORMAT ( 23X , 

1 7X , F  7. 2 / ) 

HR  I  T  E ! 6 , 2 5  5 )  VARI , VAR2 , VAR3 , V AR4 , VAR9 
FORMAT (23X , 2 8H COEFFICIENT  OF  VARIATION 
1  ,3X,F10.2,7X,F7.2///> 

CONTINUE  


{ %) ,2X,F9.2,3X*F8.2,7X,F8.2 


ISN  0090 
ISN  OC 9 1 
ISN  0092 


WRITE! 6  ,310) 
10  FORMAT! 1 M 1 ) 
STOT 


ISN  0093 


FMD 


******  EMD  OF  COMPILATION  ****** 


ISN  0702  SUBROUTINE  M AX ( L , ARRAY ,NNN> 

ISO  0703  01  or  05  ION  ARRAY! 50) 

ISO  0004 _ _ _ AO  AX  - APR  AY  (  1  ) _ _ _ _ 

iso  0005  nori=i 

ISO  0006  00  20  1=2, L 

ISO  0007  _ IF  (  A  P.  R  A  Y  (  I  )  .GT.  AMAX)  GO  TO  10 

ISO  0007  GO  TO  20 

ISO  0010  1C  NNN  = I 

ISO  0011 _ AM  AX  =  ARP  AY  (  I  ) _ 

ISO  0012  20  CONTINUE 

ISO  0013  RETURN 

ISN  0714 _ _ END _ _ _ 


*****  END  OF  COMPILATION  ****** 


ISN 

ISN 

I  S  N 

0 Of  2 

0003 

0°04 

SURROUT  I  NE  STATS  I  N,  X  ,  X  8  A  R  ,  SIGMA  ,  VAR  ) 

DIMENS  ION  X(  50) 

SUM=0. 

ISO 

0005 

DO  10  I  =  1 »  M 

1  SN 

0006 

SU M  =  $U M  +  X (  I  ) 

I  SO 

0707 

1  o 

CONTINUE  .  _ 

ISO 

0008 

AN  =N 

ISN 

0707 

BN=N-1 

ISO 

00  10 

XPAP=SUM/AN 

ISN 

0011 

SM  SGDF  =  0 » 

ISO 

00  1  2 

DU  20  1  =  1  *  N 

I  SN 

00  1'* 

SMSCDF=SMSQDE+( X( I J-XBAR) **2 

ISO 

0014 

20 

CONTINUE 

ISN 

00  15 

SI GMA= ( SMS  QDF/ BN ) **0 . 5 

ISO 

00  16 

V  A  R  =  S  I  GM,  A  /  X B  AR  *  1  00 . 

1  so 

0017 

RETURN 

ISN 

0018 

END 

******  end  OF  COMPILATION  ****** 


. 
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A39 


A40 


DETAILED  COMPUTER  ANALYSIS 


ahi 


ISM  0002  DIMENSION  T I TL E ( 20 ) , DAT E ( 5 ) , T I  ME ( 20 ) , LOAD ( 2 0 )  , S TRAI N { 20 } , STRESS l 2 0 

1),FLSTRS(2C),FLSTRN(20),FITIME(20),SAMPIE(20,5),STIFF(20),FLSTFF(2 
-  ....  -  2O),ORSTRS(2C),ORSTRN(20),ORSTFF(20).ORSBIT(20),UNirWT(?O).AIR(??). 

3AREA(22),W0RKI22),VDEF(22),BITST(22)»SBIT(22)»FLAREA(20),FLBTST(20 
4),FLSBIT(20)  »FLWORK(20) 

I  SM  0003  INTEGER  TIME  ,  SAMPLE , FA ILTM, TEMP 

ISN  0004  REAL  LOAD 

C 

c 

C  M  =  THE  NUMBER  OF  SETS  OF  SAMPLES  TO  BE  PROCESSED  DURING 

C  THIS  RUN  . 

C 

C 

ISN  0005  READ! 5, 10) M 

ISN  0006  10  FORMAT (15) 

ISN  0007  DO  300  I  I  =  1 ,  M 

C 

c 

C  TITLE  =  A  DESCRIPTION  (BO  COLUMNS  OR  LESS)  TO  APPEAR  AT  THE 

C  TOP  OF  EACH  PAGE  TO  IDENTIFY  EACH  PARTICULAR  SET  . 

C 

C 

ISN  0008  READ(5»20) (TITLE(K)  ,K=i,20> 

ISN  0009  20  FORMAT ( 20A4) 

C 

c 

C  TEMP  =  THE  TEST  TEMPERATURE  .  IN  DEGREES  FAHRENHEIT  .  AT 

C  WHICH  THIS  SET  WAS  RUN  . 

C 

C  RATE  =  THE  RATE  OF  LOAD  APPLICATION  IN  INCHES  PER  MINUTE  . 

.C 

C  DATE  =  THE  DATE  AT  WHICH  THE  SAMPLES  WITHIN  THIS  SET 

C  WERE  TESTED  , 

C 

C 

ISN  0010  READ(5»30)TEMP?RATEt (DATE(K) iK-1,5) 

ISN  0011  30  FORMAT! I 10, F 10.3, 5A4) 

C 

c 

C  ACMIX  =  THE  ASPHALT  CONTENT  OF  THE  MIX  ,  IN  POUNDS  OF  ASPHALT 

C  PFR  ONE  HUNDRED  POUNDS  OF  DRY  AGGREGATE  . 

C 

C  GAGG  =  THE  BULK  SPECIFIC  GRAVITY  OF  THE  AGGREGATE  , 

C 

C  GAC  =  THE  SPECIFIC  GRAVITY  OF  THF  ASPHALT  CEMFNT  . 

C 

C  PABS  =  THE  PERCENT  ABSORPTION  ,  IN  POUNDS  OF  ASPHALT  CEMENT 

C  ABSORBED  PER  ONE  HUNDRED  POUNDS  OF  DRY  AGGREGATE  . 

C 

C 

ISM  0012  RFAD(5,?5) ACMIX, GAGG, GAC, PABS 

ISN  00  13  35  FORMAT ( 40 X , 4 F 1 0 . 3 ) 

C 

r 

C  N  =  THE  NUMBER  OF  SAMPLES  WITHIN  THIS  SET  . 

C 

C  -  -  .... 

' 


' 


A42 


ISM  0Q14  READ! 5, 40 )N 

ISM  0015  40  FORMAT ( 15) 

ISM  0016  DO  200  1=1 rN 

C  _ _ _ 

C 

C  SAMPLE  =  THE  IDENTIFICATION  NUMBER  OF  THE  SAMPLF  . 

_ c _ : _ s _ _  _ 

C  WTAIR  =  THE  WEIGHT',  IN  GRAMS  ,  OF  THE  SAMPLE  IN  AIR  . 

C 

_ _ 1 _  c _ WT WAT  =  THE  WEIGHT  ,  IN  GRAMS  .  OF  THE  SA MPI  F  IN  WATER  . 

C 

C  DIAM  =  THE  AVERAGE  OIAMETER  OF  THE  SAMPLE  ,  IN  INCHES  . 

_ C  _  _ 

C  THICK  a  THE  AVERAGE  THICKNESS  OF  THE  SAMPLE  IN  INCHES  • 

C 

_ _  C _ L  =  THE  NUMBER  OF  LOAD  AMD  STRAIN  POINTS  PICKFD  OFF  THE 

C  LOAD  AND  DEFORMATION  STRIP  CHART  . 

C 

_ C _  FAN.  r-  -•  THE  TIME  TO  FAILURE  ,  IN  SECONDS  _ 

C 

c 

ISO  0  0  IT _ REAP  (  5,  SO)  (  S  AMR_L  F  (  I  ,K)  , K=1  ,3  )  ,  W  T  A  I  R  .  W  T  W  AT  ,  f)  I  AM ,  T  H  I  CK  .  L  ,  F  A  1 1  TM 

ISN  0018  50  FORMAT ( 3A4 , F8.3 , 2F1C .3 ,21  IQ) 

C 

_  c _ _ _ _ _ _ 

c  SPGR  =  BULK  SPTCIF IC  GRAVITY  OF  THE  COMPACTED  MIXTURE  . 

C 

_ _  C _ 

ISN  0019  SPGR=WTAIF /( WTA IR-WTWAT) 

ISN  0020  WTAC=100.*SPGR*ACMIX/(  100.+ACMIX) 

ISM  00  21 _ WTAGG=100.*SPCR-WTAC _ 

ISN  0022  WTABAC=PABS*WT AGO / 1 00 . 

C 

_  c _ i _ _ 

C  E  F  AC  =  EFFECTIVE  ASPHALT  CONTENT  OF  THE  MIX  IN  POUNDS  OF 

C  EFFECTIVE  ASPHALT  CEMENT  PER  ONE  HUNDRED  POUNDS  OF 

_ C _ DRY  AGGREGATE  . _ _ _ _ _ _ 

C 

C 


ISN  0°  23  E  F  A  C  =  1  W  T  AC  -  W  T  A  B  A  C  )  /  W  T  A  G  G  *  1 0  0  . 


ISN 

ISN 

0024 

0025 

C 

VAGG=WT  AGG/G AGG 

VEFAC= ( WTAC-WTABAC ) /GAC 

c 

c 

c 

AIR  =  AIR  VOIDS  IN  PER  CENT  . 

ISN 

0026 

c 

AIR!  I  )  =  100 . - (VAGG  +  VEFAC) 

i_c 

c 

c 

CSUBV  =  VOLUME  CONCENTRATION  OF 

THE 

AGGREGATE  . 

c 

' 

c 

BSUBV  =  VOLUME  CONCENTRATION  OF 

THE 

ASPHALT  CEMENT  . 

c 

c 

ISN 

0027 

C  SUBV=  VAGC-/  (  VAGG+VEFAC  ) 

ISN 

0028 

BSUBV^l .O-CSUBV 

A4  3 


ISN  0029  CP  V=CSU3V/!1  .0+0. 01  *  l  A  I R  (  I  )-3  .0  )  ) 

ISN  0030  CA=CPV/( 1.0-CPV) 


c 

c 


C  UNI T W T  =  UNIT  WEIGHT  OF  THE  COMPACTED  MIXTURE  , 

C  IN  POUNDS  PER  CUBIC  FOOT  . 

C 

c 

ISN  0031  UN  ITWT (  I  ) =SPCR*62 .4 

ISN  0032  l.N  =  M  +  l 

ISN  0033  NN=M+2 

ISN  0034  WR ITE ( 6 , 52 ) ( TITLE (K ) ,K  =  l , 20) , I ,NN 

ISN  0035  52  FORMAT ( 1 H 1 //////33X , 20A4/ /l 14X, 5HPAGF  ,I?t4H  OF  ,12//) 

ISN  00  36  WR  ITE! 6, 54 ) (  SAMPLE (  I ,K )  ,K=1 , 3)  ,  TEMP  ,  ( DATE! K) ,K=1 ,5) 

ISN  0037  54  FORMAT (22X, 16HSAMPLE  NUMBER  :  , 3A4 , 5X ,  1 9HTE ST  TEMPERATURE  :  ,13,19 

1H  DEGREES  FAHRENHE I T , 5 X , 7HDATE  :  ,5A4/) 

ISN  00  3  3  WR ITF( 6, 56 )D  I  AM, THICK, ACM  I X , GAGG , GAC , P A8S , EFAC 

ISN  0039  56  FORMAT! 45X, 1 1HDIAMETER  =  ,F5.3,7H  INCHES, 12X,12HTHICKNESS  =  ,F5.3, 

..  _  .  ....  .  17H  INCHES//22X  ,29HASPHAl.T  CONTENT  OF  THE  MIX  =  ,F4.2,2H  %,2,X,27HSP 

2.  GP.  OF  THE  AGGREGATE  =  , F 5 . 3 , 3 X , 2 5H SP .  GR.  0 F  THE  ASPHALT  =  ,F5. 

33/ /35X, 21  HASP HALT  ABSORPTION  =  ,F4.2,2H  S, 6X, 39HEFFECT I VE  ASPHALT 
4C0NTENT  OF  THE  MIX  =  , F4.2,2H  %/) 

ISN  0040  WRITE (6,58)SPGR,UNITWT!  I  )  , A  I R ( I  ) 

ISN  0041  58  FORMAT!  22.X, 40HBULK  SP .  GR .  OF  THE  COMPACTED  MIXTURE  =  ,F5.3,7X,14H 

1UMIT  WEIGHT  =  , F5.1.7H  P,C.F..7X,12HAIR  VOIDS  =  , F5.2.2H  %/) 

ISN  0042  WR ITF l 6 , 60 ) C SUBV, BSUBV 

ISN  0043  60  FORMAT! BOX, 36HV0LUME  CONCENTRATION  OF  AGGREGATE  =  , F 5 . 3 , 8 X , 3 4H VOL U 

1  ME  CONCENTRATION  OF  BITUMEN  =  .F5.3/I 

ISN  0044  WRITE(6, 62  IRATE, FAILTM 

ISN  0045  62  FORMAT! 35X, 18URATE  OF  LOADING  =  ,F5.3,16H  INCHES  /  M I NU TE , 8 X , 1 9HT I 

IMF  TO  FRACTURE  =  .13,80  SECONDS//) 

ISN  0046  WRITE (6, 64) 

ISN  0047  64  FOPM AT !42X,9HSTIFFN ESS »20X, 9 H STIFFNESS, 8X,7HRITUMEN»19X,10HW0RK  IN 

1  PUT  ) 

ISN  0048  WRITE!6,66) 

ISN  0049  66  FORMAT  (  2  1  X  ,  6  H  S  T  R  E  S  S  ,5X, 60S  TRAIN,  5X  ,6H0F  M I  X ,  8X  ,  4  H  T  I  MF-,  9X  ,  1 1  HOF  8 

1 I TUMEN, 7X. 6HSTRAIN,4X,9HT OUGHNE SS , 4X , 15HPER  UNIT  VOLUME  J 

ISN  0050  WR ITF (6,68) 

ISN  0051  68  FORMAT! 22X , 5H( PS  I ) ,4X, 7H{  IN/IN)  ,6X,5H(PSI ) ,6X,9H( SECONDS) ,4X,5H!PS* 

11)  ,1X, 10H(KG/SQ  CM)  ,4X , AH! IN/  IN) , 5X , 5H( PS  I ) . 8X.  I3H(  IN-!  B/CU  IN)  /  /) 

ISN  00  52  E>EN0M=3.1415  9*DIAM*THICK 

ISN  0053  STRESS !  1  )=C.0 

ISN  0054  STRAIN ( 1 )=0. 0 

ISN  0055  AREA! 1)=0.0 

ISN  0056  WORK! 1 ) =0 • 0 

ISN  0057  VDE  F ( 1 )=0.0 

ISM  0058  LOAD! 1 ) =0.0 

C 

c 

c  TIME  =  THE  TIME  ,  IN  SECONDS  ,  FROM  ZERO  TIME  TO  EACH  OF  THE 

C  POINTS  ON  THE  LOAD  AND  DEFORMATION  STRIP  CHART  . 

C  _ _  .  . 

c 


ISN  0059  1=1+1 

I  SM  0 060 _ P'EAD(5t70)  ! T I  ME  f  K )  , K=2 , L ) 

ISN  0061  70  FORMAT! 5X, 15  15) 

C 


■ 
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c 

C  LOAD  =  THE  LOAD  ,  IN  POUNDS  ,  AT  EACH  OF  THE  POINTS  ON  THE 

C  LOAD  TRACE  . 

_ C  _ _ _ : _ 

C 

ISN  0062  READ!  5, 80)  (LOAD(K)  ,  K=2  ,U 

ISM  0063 _ 80  FORM  A I  (  5Xj  1  5  F5 . 0  > 

C 

c 

_ C _ STRAIN  =  THE  STRAIN  ,  IN  TEN  THOUSANDTHS  OF  A N  INCH  /  INCH  . _ 

C  '  FOR  EACH  OF  THE  POINTS  ON  THE  DEFORMATION  TRACE  . 

C 

_ c _ _ _ ; _ 

ISN  0064  RE  AD ! 5  »  90 )  IS  TRA  I N  ( K  )  »i<  =  2 » L ) 

ISN  0065  90  FORMAT  (-5  X  »  I  5  F  5 . 0  ) 

ISN  00  66 _ DP  ICO  J-2.L _ _ _ _ _ 

ISN  0067  STRESS ( J ) =2 . *L0 AD ( J ) /DENOM 

ISN  0068  STRAIN! J ) =STRAI N ( J ) *0. 0001 

_ __C _ ; _ . _ . _ 

C 

C  CONSIDERING  AVERAGE  TENSILE  STRESS  OVER  A  ONE  INCH  GAUGE  LENGTH 

_ _  C _ AND  ASSUMING  A  POISSONS  RATIO  OF  0.33  . _ 

C 

C 

J  S_M  ?.  0  6  9 _ ST  I F  r  C  J  )  =0.912*STRFS  S  (_J  )  /  (  STR  A  IN_[J  )*.5J _ 

ISN  0070  1ST  IFF=ST I FF ( J ) 

ISN  0071  ATIME=TIMC(J  ) 


XL 


C 

C 

c 

AREA  =  AREA  UNDER  THE  STRE SS-S TR A I N  CURVE  ,  IN  PSI  . 

I  SN 

0072 

c 

AR  E  A  C  J  )  =AREA  (  J  -  1  )  4-0 .5*!  STRESS  (  J-l )  +  STRESS  (  J  )  )  *(  STR  A  INI  J  (-STRAIN  (  J- 
11)) 

ISN 

00  7? 

c 

C 

VDEF ( J )  =  ( AT  I  ME / 60 . ) /RATE 

c 

c 

c 

WORK  =  WORK  INPUT  PER  UNIT  VOLUME  , 

IN  INCH-POUNDS  PER  CUBIC  INCH  . 

I  SN 

0074 

c 

WORK! J ) =WORK ( J-l ) +0 . 5* ( LO AD ( J ) f LOAD! J-l ) ) * ( VDEF ( J) -VDEF ( J-l ) )/( WTA 

1 I R-WTWAT ) *0 . 061 

c 

c 

c 

8ITST  =  BITUMEN  STRAIN  ,  IN  INCHES  PER  INCH  . 

ISN 

00  7  5 

c 

c 

BITST!  J  )  =  STP.AIN(  J  J/BSUBV 

ISN 

0076 

TR  ENN=  3 . 0 

I  SI' 

0077 

c 

9  1  ENN-TR  E  NN 

c 

c 

SBIT  =  STIFFNESS  OF  THE  BITUMEN  , 

c 

IN  KILOGRAMS  PER  SQUARE  CENTIMETER  . 

o  o 

I 


. 


A45 


ISN  0078  SRIT! J)=( STIFF! J) /14.223) /( ( l .0*! 2. 5/ENN)*CA))**ENN 

ISM  0079  CC  NN  =  0 .03*0. 434* ALDG ( 400000 ./S3 IT ( J ) ) 

ISM  0080  DIFF=ABS( ENN-CENN) 

ISN  0081  IFIDIFF-O. 01 194.94.90 

ISM  0082  93  TP  FNN=CENN 

ISN  0003  GO  TO  91 

ISM  0084  94  CONTINUE 

ISN  0085  ISTFBT=SBITl J) 

ISN  0086  ISBTLB=SRIT! J ) * 1 4 . 22 

_  _I.S4_P0.87  WR  ITFI  6, 95)  STRESS!  J  )  .STRAIN!  J)  ,  I  ST  I  FF  .  T  IMF*  J)  ,  I  SBTL  R.  I  STFBT  .  8ITST  ( 

1 J )  »ARE A ( J ) , WORK  (  J) 

ISN  0088  95  FORMAT! 2GX,F7.2,Fll .5,4X* I8,4X, I7,4X, 18 , I8,4X,F10.5, Fll .4, 6X.F9.2/ 

1 ) 

I SN  0089  100  CONTINUE 

C 

c 

C  SUBROUTINE  'MAX'  IS  CALLED  TO  DETERMINE  THE  MAXIMUM  SECANT 

C  MODULUS  (  THAT  IS  THE  MAXIMUM  STIFFNESS  OF  MIX  FOR  FACH 

C  STRESS-STRAIN  DIAGRAM  )  . 

C 

C 

ISN  0090  CALL  M A X ( L , S T I F F , NN N ) 

ISN  0091  ORSTRSI I )=STRESS(NNN) 

ISN  0092  OR  STEM (  I  )  =  ST  RA I N ( NNN ) 

ISN  009?  CRSTFF (  11  =  ST  IFF { NNN  ) 

ISN  0094  ORSBIT!  I  )  =  SB  IT! NNN) 

ISN  0095  FLSTRS! I )=STRESS(L) 

ISN  0  0  96  FLSTPN!  I  )  =  STR  A  I  N  ( 1.  ) 

ISN  0097  FLSTFF!  I  )=ST  IFF(L) 

ISN  0098  F  L  T  IMF!  1  ) =  T I  ME ( L ) 

ISN  0099  FL  AREA!  I)  =AP.EA(  L  ) 

ISN  01C0  FLBTST!  I  )  =  B I TST ( L  ) 

ISN  0101  FI.  SO  IT!  I  J  =  SB  I T  (  L  ) 

ISN  0102  FL WORK ( I ) - WORK ( L ) 

ISN  0103  200  CONTINUE 

C 

C 

C  SUBROUTINE  'STATS'  IS  CALLED  TO  DETERMINE  THE  MEAN  ,  STANDARD 

C  DEVIATION  ,  AND  COEFFICIENT  OF  VARIATION  OF  EACH  PARAMETER 

C  BASED  UPON  INITIAL  CONDITIONS  (  THAT  IS  THE  MAXIMUM  SECANT 

C.  MODULUS  OF  EACH  STRESS-STRAIN  DIAGRAM  )  . 

C 

c 

ISN  0 1 C  4  CALL  STATS l N ,ORS TRS » X B AR5 , S I GMA 5 , V AR 5 ) 

ISN  0 1 C  5  CALL  ST ATS  I N  ,  ORSTRN , XB AR6 , S I GMA6 , VAR6 ) 

ISM  0106  CALL  STATS ( N  .ORSTFF ,XBAR7 ,S IGMA7.VAR7) 

ISN  0107  CALL  ST ATS ( N  ,  ORS B I T , X 0 AR 8 , S I GM A8 , V AR8 ) 

ISN  0108  CALL  ST ATS ( N »UN I TWT , XB AR9 , S I GMA9 , VAR9 ) 

ISN  0109  CALL  S T ATS ! N  ,  A I R , X B AR A , S I GM A A , V AR A ) 

ISN  0110  I XB AR7=  X BAR 7 

ISN  0111  I SGMA7=S ICMA7 

ISM  0112  I  X  B  AP,8=  X  BAR  8 

ISN  0113  I SGMA8  =  S I GMA 8 

ISN  0114  WRITE(6,210) ( T I TLE ( K J , K= I » 20 ) ,LN»NN 

ISN  0115  2 1 C  FORMAT ( IH1 //////33X ,20A4//l 14X, 5HPAGE  .I2.4H  DF  ,12//) 

ISN  0116  IF  IN  .GT,  6)  GO  TO  213 

ISN  0118  211  CONTINUE 

A46 


ISN  0119  WRITEI6, 212)M,(  ( SAMPLE l I , K) ,K= l  ,  3 I  ,I  =  1,N) 

I.SN  0120  212  FORMAT  (  5  0  X  ,  5  9  H  M  E  A  N  S  ,  STANDARD  DEVIATIONS  ,  ANO  COEFFICIENTS  OF  VA 

1R  I  AT  I0N/7CX  ,  3  H  F  OR  /  6  5  X »  4  H  T  H  E  ,  I2,23H  SAMPLES  IN  THIS  SER  I  ES  /  6  5X  ,  29  H 
_  .  .  2BASED  UPON  INITIAL  CON  D  I  T  I  ON  S  /  /  /  2  0  X  .  I  7H  S  A  MP 1  F  NUMRFRS  :  .TX.M1A4I 

3) 

I SN  0121  GO  TO  215 

I SN  01  22  213  CONT  INUE 

ISN  0123  WRI TE( 6, 214) N, (  ( SAMPLE  I  I , K)  ,K  =  1  ,3 ) , 1  =  1 ,6) , (  I  SAMPLE (  I  ,K ) ,K=l, 3) ,  1=7 

1  ,  N  ) 

ISN  0124  214  FORMAT { SOX , 59HMFANS  ,  STANDARD  DEVIATIONS  .  AND  COEFFICIENTS  OF  VA 

1 T I  AT  1 0  N  /  7  8  X »  3  H  T  0  R / 6  5  X »  4 1 1 T  H  C  ,  !2,23H  SAMPLES  IN  THIS  SCR  IES/65X , 29H 

2PASID  UPON  INITIAL  CONO  I  T  ION  S///20X  ,  17IIS  AMPLE  NUMBERS  :  ,3X,6(3A4) 

3/4CX ,6( 3A4 ) ) 

ISN  0125  215  CONTINUE 

ISN  0126  WRITE (6,216) TEMP, (DATE(K) ,  K  =  1 »  5 ) 

ISN  0127  21  6  r OR MAT ( / /35X  ,  19HTEST  TEMPERATURE  :  .I3.19H  DEGREES  F AHR FNHF I T. 9X , 7 

1  HD  ATE  :  , 5A4) 

ISN  0128  WRITEI6.217) 

ISN  0129  217  F0RMAT(////77Xt9HSTirFNESS.6Xt9HSTIFFNESS»8Xt4HUNIT»8X.3HAIR) 

ISN  0130  WRITE! 6, 218) 

ISN  0131  218  FCFMATI  53X,6HSTRESS  , 6X , 6HST R A  I N , 8X , 6H0F  MIX,6X,11H0F  BITUMEN,  6X,  6 

1HWE  IGHT , 6X , 5HV0  IDS) 

I SN  01 32  WR  ITE I  6, 219) 

ISN  01  33  219  FORMAT! 54X ,5H(PS I ) ,6X, 7HI  IN/IN) ,7X , 5HIPSI )  ,  8 X ,  1  OH ( KG/ S Q  CM )  , 5 X , 8H 

l(P.C.F.),6X,3H(!g)} 

ISN  0134  WRITE (6, 220)  X B AR5 , X BAR6 ,  I X B AR7 , I X B AR 8 , XB AR9 , XB AR A 

ISN  0135  220  FORMAT! //20X  ,4HMEAN,26X, F9.2 ,3X  ,F10 .6 ,5X, 13 ,5X,  I8,9X, F7 .2, 6X, F5.2) 

ISN  0136  WRITFI6, 221 )  S I GM A5 , S I GM A 6 . I SGM A7 , I SGM A 8 . S T GM A9 , S I GM A A 

ISN  01  37  221  FORMAT! /2CX,  18HSTANDARD  D E V I  AT  I  ON , 1 2X , F 9 . 2 , 3X , F 1 0 . 6 , 5 X ,  I  8 , 5 X  ,  I 8 , 9 X 

1 ,  F7 • 2 »  6  X  »  F  5 • 2 ) 

ISN  0138  WRITE(6t  222)  V AR 5 , V AR6 , V AR 7 , VAR  8 , VAR9  .  VAR  A 

ISN  0139  222  FORM  AT ( / 20 X ,  28 HCOEF F l C I E N T  OF  VARIATION  I  % )  ,  5X , F 6 . 2 , 5X , F 6. 2 , 5 X , F 1 0 

1.2,4X,F10.2,5X,F10.2,2X,F9.2) 

C 

C 

C  SUBROUTINE  'STATS'  IS  CALLED  TO  DETERMINE  THE  MEAN  ,  STANDARD 

C  DEVIATION  ,  AND  COEFFICIENT  OF  VARIATION  OF  EACH  PARAMETER 

C  BASED  UPON  FAILURE  CONDITIONS  . 

C 

C 

ISN  0140  CALL  S T AT S ( N  ,  FL S TRS , X B AR l , S I GMA 1 , V AR l ) 

ISN  0141  CALL  STATS (N»FLSTRM»XBAR2fSIGMA2tVAR2) 

ISN  0142  CALL  ST  ATS ( N , FL STFF , XBAR3 t S IGMA3 , VAR3 ) 

ISN  014?  CALL  STATS ( N  ,FLTIME, XBAR4, SIGMA4,VAR4) 

ISN  0144  CALL  S T ATS t N  ,  FL AR EA , XB ARB , S I GMAB , V ARB) 

ISN  0145  CALL  S T ATS ( M  .  FL BTST , XB ARC , S IGMAC t V ARC ) 

ISN  0146  CALL  ST ATS { N  ,  FL SB  I T , XB ARD , S I GMAD, VARD ) 

ISN  0147  CALL  ST ATS  I N , FLWORK , XB AR E , S I GMAE , V ARE ) 

ISM  0148  IXBAR3=XBAR3 

ISN  0149  I  SGM A 3=  S I GMA  3 

ISN  0150  I X  BARO=X  BARD 

ISM  0151  I SGM AD  =  S I GMA  D 

ISM  0152  WRITEI6.250) (TITLE(K) , K= 1 ,20) ,NN,NN 

ISN  ^  1  5 3  250  FORMAT { 1 H 1 // //// 3 3 X , 20 A4/ / 1 14X, 5HPAGE  ,I2,4H  OF  ,12//) 

isr:  n 1 5 4  irir  .gt.  o  go  to  253 

ISN  0156  251  CONT INUF 

I SN  0157  WR ITF [6,252 ) N, ( ( SAMPLE ( I ,K) ,K=1 ,3) , I = 1 , N ) 

T 


, 


. 

ISN  0158 


ISN 
ISN 
I  SN 


0159 

0160 

0161 


ISN  0162 


ISN 

ISN 

ISN 


01  6? 
0164 
0165 


ISN 
1  SN 
I  SN 
ISN 

ISN 


I  SN 

I  SN 
I  SN 


01  66 
0167 
01  6  8_ 
0169 

0170 
01  71 

017? 

0173 


ISM  01 74 


ISN  0175 


J  SN_ 
I  SN 


01  75_ 
0177 


ISN  0178 


ISN 
ISN 
ISN 
ISN 
ISN 
I  SN 


0179 

0180 

0181 

0182 

0183 

0184 


252  FORMAT (  50X , 59HMEANS  ,  STANDARD  DEVIATIONS  ,  A  NO  COEFFICIENTS  DF  VA 
IRI ATI0N/78X, 3HFOR/65X, 4HTHE  ,I2,23H  SAMPLES  IN  THIS  S ER  I  ES / 65X , 29H 
20 AS  ED  UPON  FAILURE  COND l T IONS// /20X , 1 7HS AMPLE  NUMBERS  :  ,3X,6(3A4) 
3) 


GO  TO  255 

253  CONTINUE 

_ WRITE  I  6 1 2  5  4 ) N  t ( { SAMPLE ( I  ,  K  ) ,K=1 .3) , 1  =  1 , 6 ) . I  I S AMPI  .  E (  I  .K )  ,K=1 . 3 1 . 1  =  7 

1 ,  N ) 

254  rOPMAT  (  50X ,  59IIME  ANS  ,  STANDARD  DEVIATIONS  ,  A  NO  COEFFICIENTS  OF  VA 
_ lALAIJGN/JPX  ,3HF0R/65Xt4HTHE  ,I2.23H  SAMPLES  H  THIS  SERIES/ 65X,29iL 

2RASED  UPON  FAILURE  CQNDI T IONS///20X , 17HSAMPLE  NUMBERS  :  ,3X,6(3A4) 
3/40X , 6 ( 3  A4 ) ) 

255  CONTINUE 


WP I TE ( 6,2  56) TEMP  ,  ( DATE(K) ,K=1 ,5) 
256  FORMAT! //35X  .19HTEST  TEMPERATURE 
1  HD  ATE  :  , 5A4  ) 


» I  3  »  1 9H  DEGREES  FAHR  ENHE I T»  9X »  7 


WRITE! 6,257) 

257  FORMAT (//// 7 8X,9HST  I  FF NFS S , 24 X , 4HUN I T , 7X , 3HAIR) 
_ WRITE!  6,25  8) _ 


2  58  FORMAT ( 53 X , 6HSTRESS , 7X ,6HSTRA IN , 7X , 7H0F  M I  X , l D X , 4HT I  ME  ,10X,6HWEIG 
IHT , 5X, 5HVCI0S) 

WP  ITE!  6, 2  5°) _ ' _ 


2  59  FORMAT ( 54X , 5  H( PS  I ) , 7X, 7H(  IN/  INI  , 7  X , 5  H ( P  S I  ) , 9  X , 9  H (  SECONDS)  , 6X,8H( P, 
1C  .  F  .  )  ,  6X  ,?)'(?)  ) 

WP  I  jr((,  ,2.60)  XRAR  1  ,  XBAR2  ,  IXBAR3  ,  XBAR4  ,  XBAP9  ,  XBARA _  _ 


2  60  r°P  MA  T ( / /?  0  X ,4HMFAN,28X,F7.2,4X,F10.6,5X,  I8,9X,F6.1,8X,F7.2,5X,F5. 
1?  ) 

_ WICIILI  6,261)  SIG  ■ '  A  1  ,  S I  GM  A  2  ,  I  SGMA?  i  S  IG  M  A4  ,  S  I  G  A9.SIGMAA _ 


261  F  F M  AT  !  / 2  0  >  »  10HST  AMDARD  DE V  I  A T  I  ON  ,  1  4X  ,  F7 . 2 , 4X  ,  F  1 0 . 6 ,  5X  ,  l  8 , 9X ,  F  6 .  I , 
18X,  F7.2, 5X , F  5 . 2  ) 

_ V' 0  LTE. !  6,26  2)  VAR  1  ,_VA  R  2  ,VAR3,VAR4 .  VAR9.VAR  A _ 

262  FORMAT ( /2CX, 28HC0EFFIC IENT  OF  VARIATION  ! X )  , 4X , F7 . 2 , 6X , F6 . 2 , 9X , F6 . 
12 ,9X,F6.2, 9X,F6.2,4X,F6.2 ) 

WP  ITE!  6 ,263) _ : _ ___ _ _ _ _ 


263  FORMAT ( ///78X, 9 H STIFF NESS ,6X, 10HW0RK  INPUT) 

WRITF! 6,264) 

264  FORMAT  I  66X  ,  7  HQ  I  TUM  EM ,  8  X  ,  2  HOF  ,  14X,  3HPEP. ) _ 

WRITE (6, 265) 

26  5  FORMAT  (  52X,9HT0UGIINESS,6X,6HSTRAIN,6X,7HBITUMEN,7X,11HUNIT  VOLUME) 
WR  ITE!  6  ,26  6 ) 


ISM  0 1 85 
0186 


_I_S  N_ 
I  S  M 
I  SN 
ISN 


0187 
0180 
01  89 


I  SN 
I  SN 


0190 

0191 


266  FORMAT!  54X,5H!PSI)  ,7X,7H(  IN/IN)  ,5X , 10H ( KG/SO  CM )  ,  4X  ,  1 3H < I N-L B/CU  I 

IN)  ) 

. _ W_R_I  TF  (  6,267) _ X 0  A RB  ,  XBARC  ,  I  XB  A  RD  ,  X  B ARE _ _ _ _ _ _ 


267  FORMAT ( / /20 X , 4 HM EAN , 2 8 X , F 8 . 4 , 3X,F10.6,4X,I8,9X,F7.2) 

WRITE! 6, 268)  S I GMAB , S I GMAC , ISGMAD, S I GMAE 

268  FORMAT! /2CX,  18HSTANDARD  PE V  I  A T I  ON , 1 4X , F 8 . 4 , 3 X , £j 0 , 6 , 4X ,  18 ,9X,F7.2) 


WR ITE! 6,269)  VA RB , V A RC , V ARD , V AR E 

269  FORMAT! /2CX, 28HCOEFFIC IENT  OF  VARIATION  ( Z ) , 6X , F6 . 2 , 5X , F 6 . 2 , 9X , F 6 . 


I  SN 

0192 

. . .  7Ajr  _  _  _  _ 

300  CONTINUE 

ISN 

0193 

WR ITE! 6,310) 

ISM 

0194 

310  FORMAT (  1H1 ) 

ISN 

0195 

STOP 

ISN 

0196 

END 

-  *  #  £  ft  * 

END  OF 

COMPILATION  ****** 

■ 


A48 


ISM  oro2 
ISM  0003 
J_S_M  00C4 


ISN  0005 
ISM  0006 
ISM  0007 
ISM  0009 
ISN  0010 
J SN  on l l 


SUBPOUT  I  ME  M  AXIL, ARRAY ,NNN  I 
DIO ENSIGN  ARRAY150) 

A  M  A  X  =  A  R  R  AY  I  7  ) _ 

NNN=2 

DO  20  1=3, L 

I F I  ARRAY (I)  .GT.  AMAX)  GO  TO  10 


GO  TO  20 
10  NNN= I 

AMAX=AP,R  AY  (  I  ) 


ISN  0012 
ISN  0013 
I SN  0014 


20  CONTINUE 
RETURN 
END 


*****  END  OF  COMPILATION  ****** 


ISN  0002 
ISN  0003 


SUBROUTINE  STATS!  N, X , X  BAR , S IGMA , VAR  ) 
DIMENSION  X ( 50 ) 


ISN  0°05 

ISN  0006 

ISN  0007 

1 .  .  V  1 _ 

DO  10  I  =  1  *  N 

SUM=SUM+X( I ) 

10  CONTINUE 

ISN  0008 

AN  =N 

ISN  0009 

RM=N-l 

I SN  0010 

XBAR-SUM/AN 

ISN 
I  SN 
ISN 
I  SN 
ISN 
ISN 


0011 

0012 
001  3 
0014 
0015 
00  16 


2C 


SMSQDF=0. 

DO  20  1  =  1 »  N 

S MS  Q  D  r  =  SMSQ OF  +  ( X ( I I-XBAR)* *2 
CONT INUE 

S I GMA  = ( SMSCDF/BN) **0.5 
VAR=SIGMA/XBAR*100. _ _ 


ISN 
I  SN 


0017 

0018 


RETURN 

END 


*****  END  CF  COMPILATION  ****** 


A49 


A50 


SUPPLY  NUMBER  5268  8  WATSONVILLE  AGG.  COMPACTION  :  75  BLOWS 


A51 


STANDARD  DEVIATION  0.1173  0.001051  110 


COMPUTER  ANALYSIS  WITH 


A5  2 

STRESS -STRAIN  PLOT 


I  SN 

0002 

(HOLMS  ION  T  I  TLE( 20)  ,0ATE( 5) r T  IME( 20) ,  LOAD! ?0 » , S TRA I N ( 22 ) t STRESS  (  2 2 

1 ) , FLSTRS(2C  )  ,FLSTRN( 20 ) ,FLT IMF( 20 ) , SAMPLEl 20i 3) , STI FF( 20) ,FLSTFF(2 
20 ) »  ORS  TRS ( 2  0 ) , OR  ST R  N ( 2  0 ) , ORS T  F F (20 ) t  OP T I M  F ( 20) . UN  I T  WT ( 20) • RUE ( 2048 

I  SN 

0003 

3)  ,  T  ITLE1  ( 10)  >  T  1  T  L  E  2  (  10) »  S  A  M  ( 3) 

INTEGER  TIME  , S A MPLE , F A I LT M , T EMP 

ISN 

0004 

PEAL  LOADjMAXSTR 

ISN 

0005 

CALL  PLOTS l BUF, 8192 ) 

ISN 

0006 

I  X  =  0 

_  LVL 

DOC  7 

_  JII-0 

c 

c 

C _ _  M  _=  THE  NUMBER  OF  SETS  OF  SAMPLES  TO  BE . PROCESSED  DIJR ING_ 


C  THIS  RUN  . 

C 

_ ; _ C _ _ _ _ _ 

ISN  0008  P.  E  A  D  (  5  » 1 0  )  M 

I SN  0009  10  FORMAT (15) 

ISN  0010  _ DO  300 _J  I _=  1  ,  M _ _ _ 

C 

C 

_ _ C _ TITLE  =  A  DESCRIPTION  (80  COLUMNS  OR  I  FSS  )  TO  APPFAR  AT  THF 

C  TOP  OF  EACH  PAGE  TO  IDENTIFY  EACH  PARTICULAR  SET  . 

C 

_ C _ _ 

ISN  0011  READI5.20) (TITLE(K) ,K=l,20) 

ISN  0012  20  FORMAT (20A4) 

_ c_ _ : _ _ _ 

c 

C  TEMP  =  THE  TEST  TEMPERATURE  ,  IN  DEGREES  FAHRENHEIT  ,  AT 

_ C _ WHICH  THIS  SET  WAS  RUN  . _ 

C 

r  RATE  =  THE  RATE  OF  LOAD  APPLICATION  IN  INCHES  PER  MINUTE  . 

_ C  _ _ _ _ _ _ _ 

C  DATE  =  THE  DATE  AT  WHICH  THE  SAMPLES  WITHIN  THIS  SET 

C  WERE  TESTED  . 

_ _  C _ , _ _ _ _ 

C 

ISN  0013  READ(5,30)TEMP,RATE, (DATE(K) ,K=1»5) 

ISM  0014 _ ?C  FORMAT  (7X,I3,F10.3, 5A4) _ 

C 

C 

_  _  C  '  _ HAXSTR  =  THF  MAXIMUM  STRAIN  ,  IN  TFN  THOUSANDTHS  OF  AN  INCH  . 

C  Tp>  BE  EXPECTED  WITHIN  THIS  RUN  .  THIS  IS  REQUIRED  TO 

C  ESTABLISH  THE  CORRECT  SCALE  TpR  STRAIN  ON  THE  CALCDMP 


C 

PLOT  TER 

.  THIS  MAY 

BE 

ANY  VALUF  ,  BUT 

THF  FOLLOWING 

,•  • 

c 

c 

c 

VALUES 

WILL  PRODUCE 

CONVENIENT  SCALES 

• 

C 

MAX  STR 

SCALE 

c 

10 

I"  =  2 

TEN 

THOUSANDS  Op 

AN 

INCH 

! 

I  NCH 

C 

25 

I "  =  5 

TEN 

THOUSANDS  OF 

AN 

INCH 

r 

INCH 

c 

50 

1"  =  10 

TEN 

THOUSANDS  OF 

AN 

INCH 

/ 

I  NCH 

C 

100 

1"  =  20 

TEN 

THOUSANDS  OF 

AN 

I  NCH 

/ 

I  NCH 

ISN  0015 

c 

c 

RE A0( 5 . 35) MA XSTR 

\ 


A5  3 


ISN  0016 


35  FORMAT ( F 10.0 ) 


C 

c 


c 

c 

ISN  0017  RF  ADI  5 , 40 ) N 

ISM  0018  40  FORMAT ( I 5 ) 

I SN  0019  DO  200  1=1  ,N 

c 

r 

C  SAMPLE  =  THE  IDENTIFICATION  NUMBER  OF  THE  SAMPI  E  . 

_  C 

C  WTAIR  =  THE  WEIGHT  ,  IN  GRAMS  ,  OF  THE  SAMPLE  IN  AIR  . 

C 

c  W THAT  =  THE  WEIGHT  .  IN  GRAMS  .  OF  THE  SAMP!  F  IN  WATER  . 

C 

C  DIAM  =  THE  AVERAGE  DIAMETER  OF  THE  SAMPLE  ,  IN  INCHES  . 

C 

C  THICK  =  THE  AVERAGE  THICKNESS  OF  THE  SAMPLE  ,  IN  INCHES  . 

C 

c  L  =  THE  NUMBER  OF  LOAD  AND  STRAIN  POINTS  PICKED  OFF  THE 

C  LOAD  AND  DEFORMATION  STRIP  CHART  . 

C 

C  FAILTM  =  THE  TIME  TO  FAILURE  .  IN  SECONDS  . 

C 

C 

ISN  0020  PEAD(5,50)(  SAMPLE!  1  ,K>  ,K=1 .3)  , WTAIR. W THAT  «DIAM. THICK. 1  . FA  II  TM 

ISN  0021  50  F0RMATI3A4 ,F8.3 ,3F10.3,2I  10) 

C 

C 

C  SPGR  =  BULK  SPECIFIC  GRAVITY  OF  THE  COMPACTED  MIXTURE  . 

C 

C 

ISN  0022  SPG R=WTAIR/( WTAIR -WTWAT) 

c 

c 

C  UNITWT  --  UNIT  WEIGHT  OF  THE  COMPACTED  MIXTURE  , 

C  IN  POUNDS  PER  CUBIC  FOOT  . 

C 

c 

ISN  0023  UNITWT I I)=SPGR*62.4 

ISN  0024  L N  =  N  + 1 

ISN  0025  NN=N+2 

ISN  0  3  26  WR  ITE(6,6C  ) I  TITLE! K)  ,K  =  1 , 20)  , I , NN, I  SAMPLE  I  I  ,K) , K=l, 3) , SPGR, UNITWT ( 

1 I )  , 01  AM, THICK 

ISN  0027  60  FORMAT!  1H1//////20X,20A4//91X,  5I1PAGE  ,  I  2 , 4H  OF  ,  I  2///44X  ,  1  7H 

l S  AMPLE  NUMBER  :  ,3A4//32X,  49HBULK  SPECIFIC  GRAVITY  OF  THE  COM 

2P  ACTED  MIXTURE  =  ,  F5 . 3/ / 38 X ,  I  4  HUN  I T  WEIGHT  =  .F5.1.22H  POUNDS 

2  PFF.  CUBIC  FOOT  ,  //45  X ,  1  1  HD  I  AM  ETC  R  =  ,  F5.3,7H  INCHES/ 

3/44/ ,  12HTHICKNFSS  =  ,  F5.3,7H  INCHES//) 

ISN  C028  WR I T  F ( 6 , 6 1  )  RAT E . F A  I LT M , T EMP ,  I D AT F ( K ) . K= 1 . 5  ) 

ISN  0029 


ISN  0030 
ISN  0031 


61  FORMAT! 2  5X , 15 H DETAILS  OF  TE S T / / 2 6 X ,  1 8HR AT E  OF  LOADING  =  ,F5.3,16H 
1  INCHES  /  MINUTE, 4X, 19HTIME  TO  FRACTURE  =  ,13, OH  SEC OND S //26X , 1 9H T6 
_  2ST  TEMPERAT  URE  :  _  ,J3 ,19H  DE  GREES  FAHRENHEIT  , 4X,9HDATE  :  .5  A4/ //) 


WRITE (6, 62) 

62  FORMAT  (60X,9HST  I  FEN  ESS  /  30  X  ,  6H  ST  R  E  S  S  ,  9  X ,  6H  SIR  AIN,  11X,6H0F  MI  X,  l  OX., 


A54 


14HT I  ME/ 31 X , 5H( PS  I ) ,9X, 7H{  IN/ IN) , 10X ,5H( PS  I ) ,9X  ,9H( SECONDS ) ///) 

ISN  0032  DENCM=3. 14159*01 AM«THICK 

C 

 c 

C  TIME  =  THE  TIME  ,  IN  SECONDS  ,  FROM  ZERO  TIME  TO  EACH  OF  THE 

C  POINTS  ON  THE  LOAD  AND  DEFORMATION  STRIP  CHART  . 

C 

C 

ISN  0033  L=L+ 1 

ISN  003 A  READ! 5, 70 ) ( T  IMEI K) . K  =  2  .1  ) 

ISN  0035  70  FORMAT ( 5X , 15  15  ) 

C 

c 

c  LOAD  =  THE  LOAD  ,  IN  POUNDS  ,  AT  EACH  OF  THE  POINTS  ON  THE 

C  LOAD  TRACE  . 

c 

C 

ISN  0036  READ!  5,  BOH  LOAD!  K)  ,K  =  2  ,L  ) 

ISN  0037  80  FORMAT! 5X. 1 5F5.0J 

C 

c 

C  STRAIN  =  THE  STRAIN  ,  IN  TEN  THOUSANDTHS  OF  AN  INCH  /  INCH  . 

C  FOR  EACH  OF  THE  POINTS  ON  THE  DEFORMATION  TRACE  . 

C 

c 

ISN  0038  STRESS (11=0.0 

I SN  0039  STRAIN! I )=C.O 

ISN  0° AO  READ! 5.9Q)  ( STPAIM(K) .K=7.l  ) 

ISN  0041  90  FORMAT ( 5X, 15F5.0) 

ISN  0042  DO  100  J=2,L 

ISN  0043  STRESS! J ) =2 . *LOAD! J ) /DENOM 

ISN  0044  STRAIN! J )  =  S  T  R  A  I N ( J )* 0.0001 

c 

C 

C  CONSIDERING  AVERAGE  TENSILE  STRESS  OVER  A  ONE  INCH  GAUGE  LENGTH 

C  AND  ASSUMING  A  POISSONS  RATIO  OF  0.33  . 

C 

c 

ISN  0045  STIFF! J  )  =0 . 9  12*ST RE SS t J ) / ( STR A I N ! J ) * . 5 ) 

ISN  0046  ISTIFF=STITF ( J) 

ISN  0047  WR  ITE! 6,95 ) STRESS! J ) , STRAIN! J }, ISTI FF , TIME! J) 

ISN  0048  95  FORMAT (26X  »  F  10 • 2 >  6X  T  FI C  .  6  »8X»  18  i9X,  15  /  ) 

ISN  0049  IPO  CONTINUE 

ISN  0050  FT  EMP=T  EMP 

C 

c 

C  SUBROUTINE  'GRAPH*  IS  CALLED  SO  THAT  THE  VALUES  REQUIRED  FOR 

C  THE  CALCOMP  PLOTTING  OF  STRESS  VERSUS  STRAIN  CAN  BE  STORED 

C  ON  THE  7  TRACK  PLOTTING  TAPE.  . 

C 

C 

ISN  0051  CALL  G R APH ( BUF , STRE SS , S TR A  I N , L .  I X . T IT L E . S AM PI  F .  I . F T F MP . M AX STR . N .  I I 

1,111,  XBAR  1,  XI3AR2) 

C 

c 

C  SUBROUTINE  'MAX'  IS  CALLED  TO  DETERMINE  THE  MAXIMUM  SECANT 

C  MODULUS  (  THAT  IS  THE  MAXIMUM  STIFFNESS  OF  MIX  FOR  EACH 

. 


A55 


C 

c 

STRESS-STRAIN  DIAGRAM  )  . 

ISN 

0052 

c 

CALL  MAX!!  . STIFF, NNN) 

I  SN 

0053  . 

OR  STRS ( I )  =  ST  RFSS ( NNN ) 

ISN 

0054 

CRSTRN! I )  =  S  T RAIN! NNN) 

j  55 1 : 

CRSTFFf  I  )  =  S  T  IFF! NNN  ) 

ISN 

0056 

UPTIME! I ) =T I MF ( NNN ) 

IS' 

0057 

FLSTP.S!  I  )  =  STRESS!  L) 

I  SN 

00  58 

FLSTRN!  I  )  =  S  T  R  A  I N ( L ) 

I  S  N 

0  059 

FLSTFF !  I  )  =ST  I  F  F ! L ) 

ISN 

00  60 

FLTIME!  I )  =  T I  ME! L ) 

I  S  N 

0061 

200  CONTINUE 

c 

c 

C 

SUBROUTINE  'STATS'  IS  CALLED  TO  DETERMINE  THF  MF  AN  .  STANDARD 

c 

DEVIATION  ,  AND  COEFFICIENT  OF  VARIATION  OF  EACH  PARAMETER 

C 

BASEC  UPON  INITIAL  CONDITIONS  (  THAT  IS  THE  MAXIMUM  SECANT 

C 

MODULUS  OF  EACH  STR C SS -STR A  I N  DIAGRAM  )  . 

ISN 

0062 

c 

c 

CALL  STATS ( N,  OR STRS, XBAR5 .SIGMA  5, VAR  5) 

I  SN 

0063 

CALL  STATS! N  ,  ORSTRN , X BAR6 , S IGMA6, VAR 6 ) 

ISN 

0064 

CALL  STATS (N  ,0R  STFF  ,XBAR7  ,  S I GMA 7 , V AR 7 ) 

ISN 

0065 

CALL  STATS ( N  ,ORT I  ME  ,XBAR8 ,S  I  GMA 8, VAR  8) 

ISN 

0066 

CALL  STATS ( N  ,UNI TWT , XBAR9 , SI GMA9 , VAR9 ) 

ISN 

0067 

I X  BAR 7  =  X  BAR  7 

ISN 

on6p, 

ISGHA7=SIGMA7 

ISN 

0069 

WRITE! 6 ,210)  (TITLE(K),K=1,20)  ,LN, NN,N, (  (SAMPLE!  I,K  )  ,K=l ,3 ) ,  1=1, N) 

ISN 

0070 

210  FORMAT ( 1H 1 ////// 20X , 20 A4 / /9 1 X ,  5HPAGE  ,12,43  OF  .I2///29X, 

159HMCANS  ,  STANDARD  DEVIATIONS  ,  AND  COEFFICIENTS  OF  VAR  I  AT  I  ON/ 57X 

2, 3HFGR/44X, 4HTHE  , I2,23H  SAMPLES  IN  THIS  S ER I E S /44X , 29H B A S E 0  UPON 
3INITIAL  CONDIT IONS///15X, I7HSAMPLE  NUMBERS  :  ,3X,8(3A4)) 


ISO  0071 _ WRITE!  6 , 2  1  1  )  TE  M  P  ,  (  0  A  T  E  (  K  )  ,K  =  1  ,5  ) _ _ _ _ _ _ 

ISN  0072  211  FORMAT! //26X , 19HTEST  TEMPERATURE  :  ,I3,19H  DEGREES  F AMR ENHE I T , 4X , 7 

1  HO  ATE  :  ,  5A41 

ISN  00  73 _ WR  I  IE  (6.2121 _ _ _ 

ISN  0074  212  FORMAT ( ////80X.9HST IFFNESS  »  2 1 X  » 4HUNIT/56X, 6 H STRESS, 6X»6HSTRAIN»8X* 

16H0F  MI X,8X ,4!  IT  I  ME, 9X,6HWE I GHT/ 57X , 5H5PSI ) , 6X, 7H{  IN/ IN) ,7X,5H( PSl ) 

_ 2. 7X,  9H  (  SECONDS)  ,  5X,  8H(  P.T.  .F.  )  ) _ ; _ 

ISN  0075  W  R  I  T  E  (  6 , 2 1  3  )  XBAP.5,  XBAR6  ,  I  XBAR7  ,  XflAR8  ,  XBAR9 

ISN  0076  213  FORMAT ( //23X  ,4HMEAN,26X ,F9.2, 3X , F10 .6, 5X, 18 , 3X, F10. 1 , 7X.F7.2/1 

I  S  N 0  C  77 _ WR  I  T  E  (6,214 1  S  I  G  M  A  5  I  GMA  6_,  IS  CM  A7  ,  S  I GMA8  ,  S  I  GMA9 _ 

ISN  0078  214  FORM  AT ( 23X ,  1 8HST  ANOARO  DE V I  AT  I ON , 12 X ,F9 . 2 , 3X , F 10, 6, 5X , I  8, 3 X , F 1 0 . 1 , 

17X.F7.2/1 

ISN  QP  79 _ WRITE  (6, 21  5)  V  AR  5  ,  V  AK6  ,  VA  R7  ,  VAR  8  ,  VAR9 _ 

ISN  0080  215  FORMAT !23X,  2 8HC0 EFF I C I  ENT  OF  VARIATION  ( %) , 2X.F9. 2. 3X,F8.2,7X,F8.2 

1  ,3X  ,F10. 2, 7X » F  7 . 2 / / / 1 

C 

C  SUBROUTINE  'STATS'  IS  CALL  CD  TO  DETERMINE  THE  MEAN  ,  STANDARD 

_ _ C  _ F^VIATICN  ,  AMO  COEFFICIENT  OF  VARIATION  OF  EACH  PARAMETER _ 

C  BASED  UPON  FAILURE  CONDITIONS  . 

C 

c _ _ :  _ 

ISO  0081  CALL  STATS(N,FLSTRS»XBAR1»SIG'MA1,VARI) 

IS*'  0  ^  8  ?  CALL  STATS  (NfFLSTRN,X3AR2  ,  S  I GMA2  ,  VAR2  1 


I 
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ISM  0083  CALL  STATS { N  , FLST FF , XB AR3 » S I GMA3 , V AR3 ) 

ISM  0084  CALL  5T ATS ( N  , FL T I ME » XB AR4 , S I GMA4 , VAR4 ) 

ISN  0085  IXI3AR3  =  XI3AP3 

ISM  0086  I S  Of A  3  =  S I GMA  3 

ISM  0087  111=111+1 

ISN  0088  CALL  GR APH ( BUF , STRE SS , S TR A  I N , L ,  I  X , T I T L E , S AMPL E,  I , F TE MP , MAXS TR , N , 1 1 

_  1,1 II,  XRAR1 , X8AR2 ) 

JSM  0  08  9  W  P I T  F  (  4  »  2  5  C  )  (TITLE! K) »  K  =  1 . 2  0  )  » NN.NN ,N, (  (SAMPLE!  I  ,K)  ,K=1  ,3) , t  =  l,N) 

I5M  cr  nr,  25C  FORMAT!  liu  //////20X.20  A4//91X  ,  5MPAGE  , I2,4H  OF  ,  12///29X, 

.  _  .  1  59  IP' p  AMS  j,  iTAMpAPO  DEVIATIONS.  *  AND  CnEFFICIFNJS  OF  V  AKI  AT  I  Q_N/5  7X_ 

2 1 3  H  T  0  R / 4  4  X »  4  HT  H  E  .  I2.23H  SAMPLES  IN  THIS  SFR I ES /44X , 2 9H BA SEO  UPON 
3FAILURE  COMO  IT  I0NS///15X , l 7HS AMPLE  NUMBERS  :  , 3  X , B  C  3 A  4 )  ) 

ISN  0091  WRITE! 6, 251 ) TEMP. ! DATFIK)  TK=1  f 5 J 

ISN  0092  251  FORMAT !//26X .19HTEST  TEMPERATURE  :  .I3.19H  DEGREES  FAHR  ENHE I T »  4X »  7 

1  HD ATE  :  ,5A4) 

ISM  009^  '  W R  I T E ( 6.25?) 

ISN  0094  252  F 0 RM AT ( / // / 8 OX , 9 HS T I FFNE S S , 2 1 X , 4HUN I T / 56 X , 6H STK ES S , 6X , 6H S TR A  I  N , 8X  l 

16H0F  M  I  X.  8X  »  4H 1 1  ME  »  9X.6HWE  IGHT/57X  ,  5HI  PS  I  ),6X,  7  HUN/IN), 7X,5H(PSI> 

2 . 7X.9H! SECONDS ) , 5 X . RH ( P . C  .  F . ) ) 

ISN'  0095  WRITE(6»253)  XBAR1 . XBAR2, I XBAR3 ,X BAR4, X3AR9 

ISN  0096  253  FORM  AT ( / /2 3X , 4HME AN , 26X , F 9. 2 , 3X , F 1 0 . 6 , 5X , I 8 , 3X , FI  0 . 1 , 7X , F7 . 2/ 1 

ISN  0097  W  R I T  E  (  6  »  2  5  4 )  S  I  GMA 1 , S I GMA 2 . I SGM A3 . S I GM A4 . S I GM A9 

ISN  0098  254  FORMAT ( 23X , 1 8HSTANDAR0  DE  V  l  AT  1 0 N ,  1  2  X  ,  F 9 . 2 , 3  X  ,  F  i  0 . 6 ,  5X  ,  I  8,  3 X  ,  F  1 0 . 1 , 

17X.F7.2/) 

ISN  0099  WRITE(6»255)  V AR 1 , VAR2 . VAR3 . VAR  4, VAR9 

ISN  0100  255  FORMAT!  23X.2  8HC0EFF ICIENT  OF  VARIATION  ! % ) , 2X , F 9 . 2 . 3X , F 8 . 2 , 7X , F8 . 2 

1 ,3X,F10.2,7X,F7.2///) 

ISM  0101  300  CONTI MUF 

ISM  0102  CALL  PL OT ( 0 . 0 , 0 . 0 , 999 ) 

ISN  C 1 03  WRITE! 6,310) 

ISN  0104  310  FORMAT! 1  HI  ) 

ISN  0105  STOP 

ISN  0106  END 

DOOMS  FOR  EXTERNAL  REFERENCES 

*****  Emd  of  COMPILATION  ****** 

nil 


ISM  3002 
ISM  00 G3 
ISN  0004 


ISM  0005 
ISN  0006 
.  I  S  M  0CQ7 
ISN  0009 
ISN  0010 
I SN  0011 


SUSP  NUT  I NE  M AXIL , ARRAY ,NNN) 
DIMENSION  ARRAY ( 50) 

AM  AX  =  ARRAY  (  2  I _ _ 


10 


ISN  0012 
ISN  0013 
ISN  0014 


NNN  =  2 

DO  20  1=3, L 

I  F  (  ARRAY  I  I  )  » GT  .  AMAX?  GO  TO  10 

GO  TO  20 
NNN=  I 

AMAX=ARRAY I  I  ) 


2C 


CONTINUE 

RETURN 

END 


*****  END  qp  COMPILATION  ****** 


ISN 

I  SI' 
ISN 

0002 

0003 

0004 

SU GROUT  I  ME  STATS! N,X,XBAR , SIGMA, VAR) 
DIMENSION  X ( 50 ) 

SUM=C. 

• 

ISN 

0005 

DO  10  1=1, N 

I  SN 

0006 

SUM=SUM  +  X<  I ) 

I  S.N 

0007 

10  CONTINUE 

ISN 

OOOB 

AN=N 

ISN 

00  09 

BN=N-l 

ISN 

0010 

XR AR=SUM/AN 

ISN 

0011 

SMSQDF=0. 

ISN 

0012 

DO  20  1=1, N 

ISN 

0013 

SMSQDF=SMSQDF+( XI I J-XBAR) **2 

ISN 

0014 

20  CONTINUE 

ISN 

0015 

SIGMA=( SMSQDF/BN )**0.5 

I  SI' 

0016 

VAR=SIGMA/XBAR*100. 

I  SN 

0017 

RETURN 

ISN 

0018 

END 

*  if  f/L  *  $ 

END  OF 

COMPILATION  ****** 

. 


. 
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ISN 

ISN 

0002 

0003 

SUB  ROUT  I  ME  GRAPH!  BUF  ,  STRE  SS  ,  STRA  I  N , L  ,  IX, T l TLE , S  AMPLE  ,  [  ,  FTEMP",  MAX  ST  '  ' 

IP  »N, I I  , I I I  , X  BARI , X  3 AR2 J 

- HI MC  NS  I  ON  STRESS!  2?)  •  S  T  R  A  I  N  (  2  2  )  .TIT!  F(20)  .SAM PI  Fl2n.3I.TfTI  FlflO). 

I  SN 
ISN 

0004 

0905 

1  T  I  T  L  F  2  (  1 C  )  ,  S  AM  (  3  ) 

REAL  UAXSTR 
_  .STRESS  (  L  +  l  )  =0.0 

I  SN 

0006 

S.T  R  A  I N  (  L  + 1  )  -  C .  0 

I  SN' 

0^7 

STPFSS! L+2 ) =1C0.0 

ISN 

nvo 

STRAIN! L +2  )  =  ^ A XSTR *0.0001/ 5.0 

ISN 

r  o  o n 

I r (  I  .  GT  .  1  )  GO  TO  50 

I  SN 

Nil! 

CALL  RFC  T ( L . 0,0. 0,11. 0,8. 5, 0.0, 3) 

I  ?f 

0  0  12 

CALL  PLOT { 2 . 5,2 . 5,-3) 

I  SN 

0013 

CALL  AXIsIc. 0,0.0, ’STRAIN  ,  IN/  IN' ,-14,5.0, 0.0, STRAIN!L  +  1) , STRAIN! 

1L+2)  ,  10. 0) 

ISN 

0014 

CALL  AXIStO. 0,0.0, 'TENSILE  STRESS  .  P S I  ’  . 2 0 . 7 . 0 . 90 . 0 . ST R F S S (1  + 1) . S 

1TRESS1L+2) ,20.0) 

ISN 

0015 

CALL  PLOT ( 5. 0,0. 0,3) 

I  SN 

0016 

CALL  PLOT (5. 0,7. 0^2) 

I  SN 

0017 

CALL  PLOT(0.0,7.0,2) 

I  SN 

0018 

X  =  0.C 

I  SI' 

0010 

DO  1  J  =  1 ,6 

ISN 

0020 

CALL  SYMBOL! X, 7. 0,0. 10, 13,0. 0,-1) 

ISN 

0021 

X  =  X  +  1.0 

I  SN 

0022 

1_ continue 

ISM 

0023 

Y=  0 . 0 

ISN 

00  24 

DO  2  J  =  l,  15 

ISN 

0025 

CALL  SYMBOI ( 5 . 0 , Y . 0 . 1 0 , 1 3 , 90 .0 . -I ) 

ISN' 

00  2  6 

Y= Y+0 . 5 

I  SI’ 

0027 

2  CONTINUE 

I  SI! 

0028 

DO  3  J=1  ,  10 

ISN 

0029 

T[ TLE 1 ! J)=T I TLE! J) 

I  SN 

0030 

3  CONTINUE 

I  S  N 

00  31 

K=  0 

I  SN 

00  32 

DO  5  J  -  1 1 ,  2 0 

I  SN 

0033 

K  =  K+  1 

ISN 

0934 

T I TLC2!K)=TITLF( J) 

ISN 

0935 

5  CONT  INUE 

I  SN 

0036 

CALL  SYMBOL ( 0.5 ,6. 8 ,0. 12, T I TL El ,0.0 ,40) 

I  SN 

90  3  7 

CALL  SYMBOL  (0.5. 6. 6,0.  12,  T.ITLE2 ,0.0  ,40) 

I  SN 

0038 

CALL  SYMBOL! 1 .5, 6.4,0. 10, 19HTEST  TEMPERATURE  :  ,0.0,19) 

ISN 

0039 

CALL  NUMBER! 3. 4, 6. 4,0. 10, FTEMP, 0. 0,-1) 

I  SI! 

0040 

CALL  SYMBOL ( 3.7, 6.4, 0. 10, 1HF ,0.0, 1) 

ISN 

0041 

IF  { STRESS (L)  .GT.  400.)  GO  TO  30 

ISN 

0043 

Y 1=6 . 0 

ISN 

0044 

Y  2  =  5 . 8 

I  SN 

00  45 

Y3  =  5 . 5 

ISN 

0046 

Y4=  5 . 2  5 

ISM 

0  047 

Y5  =  5 .0 

I  S  N 

0048 

GO  TO  40 

ISM 

0049 

30  CONTINUE 

I  S  N 

9050 

AN  =N 

I  SI! 

0051 

Y1  =AN<fO  .2+1 . 0 

ISN 

0052 

Y2  =  Y  1-0 .2  - 

ISM 

0053 

Y3  =  Y  2-0 . 3 

I  SN 

00  54 

Y  4=Y3- 0 .25 

* 
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ISN  0055  Y5=Y4-0.25 

ISM  0056  40  CONTINUE 

ISN  0057  CALL  SYMBOL! 2.9, Y3, 0.12, 6HLEGEND, 0.0, 6) 

ISN  00  5  8  CALL.  S  YMBOL  (  2 . 2  ,  Y4, 0 . 1  0, 2  4HSYMB0L  SAMPIF  M  J  MB  FR  .  0 . 0 . 2  4  ) 

ISM  0059  Y=  Y 5 

ISN  0060  ,  50  CONTINUE 

ISN  0061  _  _  I  Ft  1. 1  .  FQ.  Ill)  GO  TO  100 

ISN  0063  DO  7  J -  1  ,  3 

ISN  0064  SAM{  J )  =  SAMI’  L  E  (  I  ,  J  ) 

ISN  ro65  7  rrf.>TMjtiF 

ISN  0066  r>YMIU.-I+12  7 

I  r'r  C 0 6 7  CALL  L  INC!  STRAIN, STRESS, 1. ,  1  ,  1  ,  l  SYMBl) 

ISN  .0068  _  CALL  SYMBOL  (  2 . 5  .  Y  ,  0 . 1  0  ,  I  SYMRI  ,  0 . 0  ,  - 1  ) 

ISN  0069  CALL  S YMBOL ( 3 . 4 , Y , 0 . 10 , S AM , 0 . 6 , 1 2 ) 

ISN  0970  Y=Y-0. 2 

ISM  0071  RETURN 

ISN  007?  100  CONTINUE 

ISN  0073  CALL  SYMBOL ( 2 . 0 , Y l , 0 . 1 0 , 2 IHAVE  FAILURE  STRESS  =  ,0.0,21) 

ISN  00  74  CALL  NUMBER  (4.1  ,Y1 ,0.1  0,  XRAP.1  tO.Of?) 

ISN  0075  CALL  S YMBOL ( 2 . 0 , Y2 , 0 . 1 0 , 2 1  HAVE  FAILURE  STRAIN  =  ,0.0,21) 

ISN  0076  CALL  NUMBER ( 4. 1 , Y2 , 0. 1 0, XB AR2 ,0 .0 , 5) 

ISN  00  77  i  x  =  I  X  >  J 

ISN  0078  IF  (lx  .EC.  2)  GO  TO  10 

ISN  0080  CALL  PLOT ( -2 . 5 , 8 . 5 , -3) 

ISN  0081  RETURN 

ISN  r 082  10  CALL  PLOT ( 6 , 0, - 1 3 . 5 , -3 ) 

ISN  0083  IX  =  0 

I  SI'  0084  RETURN 

ISN  0085  END 

*****  END  OF  COMPILATION  ****** 

• 

COMPUTER  ANALYSIS  WITH  STIFFNESS -STRAIN  PLOT. 
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ISN  30  02  DIMENSION  TITLF(2Q)»DATE(5)rTIME(2C)  »  L  0  A  D  (  2  C  i  rSTPAIN(22)*STRESS(22 

1 )  .FLSTRS (20)  ,FLSTRN(23 ) , FLTIMl! 20  )  .SAMPLE! 2C  ,  SI  ,STI FF ( 20) .FLSTFFl 2 

- AO  )  t  OP s  TR li2_01j_nRS.T R N.L 2 0  )  ,  QRST_F_F(  2  C  )  .  OP T  I  ME  (  2  0  )  ,  0 v  I  T  W T  (  2  0  )  .  R UE7  2QV8 

3 )  , T  I  TIE  1  (  I  0  )  ,  T I TLE2 ( 10 ) , SAM! 3 ) , PLRESS ! 2  3 )  , PLRA I Nl 23 ) 

ISN  0003  INTEGER  T I M F , S AMPL E , F A  I L TM , T EMP 

_ LS N  0004 _ REAL  L PAD,  MAX STR,  M A X  SJ F _ 

ISN  0005  CALL  PL OT S ( D UF , 8  192 ) 

ISN  0306  IX=C 

ISN  OOP  7 _ 1 1  _ 

C 

c 

- £ _ H  =  the  NUMBER  OF  S  E  T  S  _  0  F_S  A  MPLES  TO  B  E  _P  ROCES  SEP  PUR  IN  G 

C  THIS  RUN  . 

C 

_ C_ _ 

ISN  0008  READ! 5 , 10 ) M 

ISN  0009  1C  FORMAT ( 15) 

_ ISN  0010 _ DC  300  1  1  =  1  . M _ _ 

C 

C 

_ _ _c _ TITLE  =  A  DESCRIPTION  (80  COLUMNS  OR  IFSS)  TO  APPFAR  AT  THF  ■ 

C  TOP  OF  EACH  PAGE  TO  IDENTIFY  EACH  PARTICULAR  SET  . 

C 

_ C _ : _ 

ISN  0011  RE  AD ( 5 , 20 ) (TITLE(K)  ,  K-l , 20 )  . * 

ISN  0012  20  FORMAT ( 20A4) 

_ C _ _ 

C 

C  TEMP  =  THE  TEST  TEMPERATURE  ,  IN  DEGREES  FAHRENHEIT  ,  AT 

_ C _ WHICH  THIS  SET  WAS  RUN  .  _ 

C 

C  RATE  =  THE  RATE  OF  LOAD  APPLICATION  IN  InC.-ES  PFR  MINUTE  . 

_  c _ ; _ _ 

C  DATE  =  THE  DATE  AT  WHICH  THE  SAMPLES  WITHIN  THIS  SET 

C  were  tested  . 

_ c _ ; _ ' _ _  _  _ _ _  __  _ 

C 

ISN  0013  RE  AD  (  5  »  3  0 )  TF  MP  ,  R  A  TE  »  ( (3  AT  E  (  K  )  ,  K  - 1 , 5 ) 

ISN  0014 _ 30  FORMAT!  7X,  I3.F1Q.3,  5A4) _ _ ; _ 

C 

C 


c 

MAX  STR  =  THE  MAXIMUM  STRAIN  ,  IN  TEN  THOUSANDTHS  OF  AN 

INCH  , 

c 

c 

c 

TD  BE  EXPECTED 
ESTABLISH  THE 
PLOTTER  .  THIS 

WITHIN  THIS  RUN  .  TmIS  IS  REQUIRED  TO 
CORRECT  SCALE  FOR  STRAIN  ON  THE  CALCOM? 

MAY  BE  ANY  VALUF  ,  -UT  THE  FOLLOWING 

c 

c 

c 

VALUFS 

WILL  PROOUCE  CONVENIENT  SCALES  . 

c 

MAX STR 

$  C  -  _  c 

c 

c 

10 

1" 

=  2  TEN  THOUSANDS  OF  AN  INCH  f 

INCH 

c 

25 

1" 

=  5  TEN  THOUSANDS  Gr  - N  INCH  / 

IN.  3 

c 

50 

1" 

=  10  TEN  THOUSANDS  OF  ;N  INCH  / 

INCH 

c 

100 

1" 

=  20  TFN  THOUSANDS  D '  \N  INCH  / 

INCH 

c 

c 

ISM  0015  PFAD( 5. 35) MAXSTR 
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ISN  0016 


35  FORM AT { F 10 . 0 ) 


C 

c 


-  - - - — - -  -  -  ■  ■  *  vv.'  1  ^  ’  *  '  •  - V"  _ 1  i:  J E_  Ji.  J  _1_U_  LU I— _^r  LU  1  LU _ 

c  WITHIN  THIS  RUN  .  THIS  IS  REQUIRED  TO  ESTABLISH  A 

C  CONVENIENT  SCALE  FOR  STIFFNESS  ON  THE  CALCOMP 

-  -  C  .  .  PLOTTER  .  THE  VALUE  CHOSEN  SHQUl D  RE  OIVISABLE  BY  7 

C  SINCE  THE  PRE-DETERM  I  NED  LENGTH  DF  THE  Y-AXIS  IS 

C  EQUAL  TO  7  INCHES  .  THE  FOLLOWING  VALUES  WILL  PRODUCE 

..  _ _ _ C  CONVENIENT  SCALES  (  NOTE:  IF  THE  MAXIM  JM  STfFFNFSS  IS 

C  NOT  KNOWN  THE  BEST  PROCEDURE  TO  FOLLOW  IS  TO  CHOSE  A 

C  VERY  LARGE  VALUE  TO  ENSURE  THAT  ALL  STIFFNESS  VALUES 

-  .  C  _  .  WILL  FALL  WITHIN  THE  LIMITS  OF  THF  PLOT  )  . 

C 

C 

F  MAX  STF  _ SCALE 

C 

C  700000  1"  =  100000  PSI 

-  _  _  . . c  .  ..  1400000  1"  =  20C000  PSI 

C  2100000  1"  =  300000  PSI 

C  2800000  1"  =  400000  PS I 

C_ 3  500000_ 1"  =  500000  PSI 

c 

c  •  ETC 

C 

C  7000000  l"  =  1000000  PSI 

c 

c 

ISN  0017  RE  AD ( 5 , 3  6 )  MAXSTF 

ISN  0018  36  FORMAT! F10 .0 ) 

C 

C 

C  N  =  THE  NUMBER  OF  'SAMPLES  WITHIN  THIS  SET  , 

C 

C 

ISM  00 ID  READ! 5, 40 )N 

ISN  0120  40  FORMAT! 15) 

ISN  00  21  DO  200  1  =  1 >  N 

C 

C 

C  SAMPLE  =  THE  IDENTIFICATION  NUMBER  OF  THE  SAMPLE  . 

C 

C  WTAIR  =  THE  WEIGHT  ,  IN  GRAMS  ,  OF  THF  SAMPLE  IN  AIR  . 

C 

c  WTWAT  =  THF  WEIGHT  ,  IN  GRAMS  ,  OF  THF  SAMPLE  IN  WATER  . 

C 

C  DIAM  =  THE  AVERAGE  DIAMETER  OF  THE  SAMPLE  ,  IN  INCHES  . 

C 

C  THICK  =  THF  AVEPAGE  THICKNESS  OF  THE  SAMPLE  ,  IN  INCHES  . 

C 

C  L  =  THF  NUMBER  OF  LOAD  AND  STRAIN  POINTS  PICKED  OFF  THE 

C  LOAD  AND  DEFORMATION  STRIP  CHART  . 

C 

C  FAILTM  =  THE  TIME  TO  FAILURE  ,  IN  SECONDS  . 

C  _  ' _  _  __ 

C 

ISN  0022  READ! 5, 50)  (SAMPLE!  I , K)  , K=  1 , 3 )  » WTA I R , WTW AT , 0 1  AM , TH I CK , L , FA  I L TM 

.  * 
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ISN  0023  50  FORMAT ( 3A4 , F 8. 3 , 3F 10. 3, 2  I  10) 

C 

c 

- - c - ST  PR  =  BULK  SPEC  I  F  I C  GRAVITY  OF  THE  COMPACTED  MIXTURF  . _ 

C 

c 

I  S  N_0  0  2  4 _ SPGR=WTAIR/ (  WT  A  I  R -WTWAT ) 

C 

C 

_ C _ UN  I  T  WT  =  UNIT  WEIGHT  OF  THF  COMPACTED  MIXTURE  . _ 

C  IN  POUNDS  PFR  CUBIC  FOOT  . 

C 

_ c _ _ 

I.SN  C 0  2  5  UN  I  T WT (  I  )  =  S P OR *62.4 

ISM  0026  LN=N+l 

ISM  0027 _ NN=N+2 _ 

ISM  0028  WR ITC(6,6C)(TITLE(K),K=1»20)»I»NN»(SAMPLE(I»K)»K=1»3)»SPGR»UNITWT( 

1  I  )  ,  0  I  A  M ,  T  H  I C  K 

ISM  ...0.0.29 _ 60  FORMAT  I  lit!  // ////20X , 70A4//91 X  r  _  .5  .HP  AGE  ,  12  ,  4H  OF  ,  l  2/ // 4  4X  ,  I  7H 

1  SAMPLE  MUMPER  :  ,3A4//32X,  49  H BULK  SPECIFIC  GRAVITY  OF  THE  CO  M 

2  PACT  ED  MIXTURE  =  ,  F5.3//38X, I4HUNIT  WEIGHT  =  ,F5.1,22H  POUNDS 

_ _ _ 2  PEP.  CURIC._FOn.I..i _ //45X,  1IHDMMETER  .  F5.3.7H  INCHE  SV_ 

3/44X , 12HTU ICKNFSS  =  ,  F5.3.7H  INCHES//) 

ISM  0030  W  R  I  T  F  (  6 , 6 1 )  RATE, FAILTM, TEMP,  (DATE(K)  , K= 1 , 5  ) 

ISN  0031 _ 61  FORMAT!  2  5X  »  1  5HDE  TAILS  OF  TE  ST//26X  ,  18HRATE  OF  L  OAO  [  MG.  .=.  ,  F5  .  3  ,  I6H 

1  INCHES  /  MINUTE, 4X, I9HTIME  TO  FRACTURE  =  ,13,33  SEC0NDS//26X , 1 9H  T  c 

2  S  T  TEMPERATURE  :  ,I3,19H  DEGREES  FA  HR  EM HE  IT ,4X,9HDATE  :  , 5A4 ///) 

ISN  00  3  7 _ WRITE  I  6, 67) _ _ 

ISM  0033  62  FORM  AT (60X, 9HST I FFMESS/30X, 6HSTRFSS ,9X, 6HSTRAI N, I IX , 6HDF  MIX.IOX, 

14HTI ME/31X ,5H( PS  I) ,9X,7H(  IN/ I M)  ,10X ,5  H( PS  I ) ,9X , 9H( S ECOMDS ) ///) 

ISN  0034 _ 0ENQM=3..14I5  9*0I  AM*THICK  _ _  _ _ 

C 

c 

_ c _ TIME  =  THE  TIME  ,  IN  SECONDS  ,  FROM  ZERO  TIME  TO  EACH  OF  THE _ 

C  POINTS  ON  THE  LUAD  AND  DEFORMATION  STRIP  CHART  . 

C 

_ c _ _ _ 

ISN  0035  READ! 5,70) ( T  IMF(K) ,K=l ,L) 

ISN  0036  70  FORMAT ( 5X , 15  15) 

_ _ C _ _ _ 

c 

c  LOAD  =  THE  LOAD  ,  IN  POUNDS  ,  AT  EACH  OF  THE  POINTS  ON  THF 

_ C _ LOAD  TRACE  . _ 

C 

c 

ISN  0037 _ READ!  5,30)  (  L  PAD  (  K  ?  ,  K  =  l  ,L) _ _ _ 

ISN  0038  80  FORMAT! 5X, 15F5.0) 

C 

c  STRAIN  =  THE  STRAIN  ,  IN  TEN  THOUSANDTHS  OF  AN  1 NCH  /  INCH  , 

C  FOR  EACH  OF  THE  POINTS  ON  THE  DEFORMATION  TRACE  . 

_ C _ 

C 

ISN  0039  R E A D ( 5 , 9 0 )  ( STPAIN(K) ,K  =  1 ,L) 

I  SN  0040 _ 90  FORMAT  (  5X  ,  15F5.0) _ _ _ _ _ 

ISN  0041  DO  100  J  -  1  ,  L 

ISN  0042  STRESS (J )=2.*L0AD(J J/OENOM 


' 


■ 


A6  3 


ISN  0043 


C 

C 

C 


STRAIN! J)=S TRAIN! J)*0.CD0l 


C 

c 

c 


CONSIDER.  LNG__A A/E  R  AGE  _TJL  N5IL£_.S  J  0  E  S  S  OVER  A  ONE  INCH  0  All  OF  I  FNG  TJi 


AND  ASSUMING  A  POISSONS  RATIO  OF  0.33 


ISN  0044 
ISN  0045 
ISN  0046 
ISN'  004  7 
ISN  0048 
ISN  0049 


ISN  0050 


C 

c 

c 


c 

c 

c_ 

c 


STIFF(J)=0.9l2’f'STRESS(J)/(STRAlN(J)*.5) 

I STI FF=ST IFF ( J) 

_ WP.m:l6  T«F331JJ_,  S  TRA  I  N(JJ  ,  I  S  TJ  FF  ,  F  I  MF  i  J  ) 

<’5  FORMAT  (  2  6  >'  ,  f  10 . 2  ,  OX  ,  FI  0 . 6 , OX  ,  18 ,9X,  15/  ) 

100  CONTINUE 
FT  EMP=T  FNP 


.SUBRPUT  I  NX ..  '.GRAPH  '  IS,  C  AjL  LLP  „S  XlLAI  _  XliL  VALUES-. REQUIRED  FUR _ 


Tilt  CALCUMP  PLOTTING  OF  STIFFNESS  VERSUS  STRAIN  CAN  BE  STORED 
ON  THE  7  TRACK  PLOTTING  TAPE  . 


CALL  GRAPH! BUF, STIFF  ,  STRA  I  N  ,  L  ,  I  X  ,  T  ITLE  ,  S.  AMPLE,  I ,  FT  E  MP ,  M  A  X  S  T  R  ,  N  ,  l  I 


C 

c 

c 

SUB  PHUT  I N  E  'MAX'  IS  CALLED  TO  DETERMINE  THE 

MAX  I  MUM 

SECANT 

C 

MODULUS  (  THAT  IS  THE  MAXIMUM  STIFFNESS  OF 

M  I  X  FOR 

EACH 

C 

c. 

STRESS-STRAIN  DIAGRAM  )  . 

ISN 

0051 

C 

CALL  MAX(L, STIFF, NNN) 

I  SN 

0052 

ORSTRS!  I  ) =  S T R F $ s ( N N N ) 

I  SN 

0053 

ORSTRN ( I  )  =  S  T  R A  I N ( NNN ) 

ISN 

0054 

ORSTFF (  I  )  =  ST  IFF ( NNN) 

ISN 

0055 

UPTIME!  I ) =  T I MF! NNN ) 

I  SN 

0056 

FLSTRS! I ) =ST  RESS (L) 

ISN 

0057 

FLSTPN!  I  )  =  S  T  R  A  I N ( L ) 

I  SN 

0058 

E L  S T F F (  I  )  - S T  IFF ( L  ) 

I  SN 

0059 

FIT  I  ME (  I )=T I  ME! L) 

ISN 

0060 

c 

200  CUNT l NUF 

c 

c 

SUBROUTINE  'STATS'  IS  CALLED  TO  DETERMINE  THE  MEAN 

,  STANDARD 

c 

DEVIATION  ,  AND  COEFFICIENT  OF  VARIATION  OF 

EACH  PARAMETER 

c 

BASED  UPON  INITIAL  CONDITIONS  !  THAT  IS  THE 

MAXIMUM 

SECANT 

c 

c 

MODULUS  OF  EACH  STRESS-STRAIN  DIAGRAM  )  . 

I  SN 

00  61 

c 

CALL  STATS (N  , OR S TRS , XB AR5 , S I GMA5 , V AR5 ) 

I  SN 

0062 

CALL  STATS (N  , ORSTRN , XBAR6, SI GMA6, VAR 6) 

I  SN 

0063 

CALL  STATS! N  .ORSTFF ,X3AR7 , S I GMA7 , VAR7) 

I  SN 

0064 

CALL  STATSIN  ,0R T I  ME , X B AR8 , S I GMA8 , VAR8 ) 

ISM 

0065 

CALL  STATS ( N  ,UNI TWT , XB AR9 , SI GMA° , VAR9 ) 

ISN 

0066 

I X  B AR7  =X  BAP  7 

ISN 

0067 

I S  GM A 7  =  S I GMA  7 

ISN 

0068 

WR I T  F ( 6 , 2 1 0 )  (TITLE! K),K  =  I, 20)  ,LN,NN,N,(  (SAMPLE 

( I,K) , K- 

I_t_3J  ,_I=_lj_NJ_. 

ISN  0069 


210  FORMAT ( 1H 1 / / // / /20X , 20 A4 / /9 1 X , 

•  1 59HME AN  S  ,  STANDARD  DEVIATIONS 


5HPAGF  i  I  2  »  4H  OF  ,  I2///29X, 

AND  COEFFICIENTS  OF  V AR  I  A T  I  ON / 57 X 


A64 


ISN 

ISN 


0070 

0071 


ISN' 

ISO 


0072 

0073 


ISO  0074 


ISO 
I  SO 
ISN 


0075 

0076 

0077 


ISN  0078 


2 , 3HF  OR/44  X , 4  HTHE  »I2,23H  SAMPLES  IN  THIS  S E R I F S / 44X , 29H BA S ED  UPON 
3INITIAL  COND  IT  I0NS///1  5X  ,  1  7IIS  AMPLE  NUMBERS  :  ,3X,8(3A4») 

WP.I  TEl  6,21 1  )  TEMP,  (  DATE  IK)  ,K=1  ,5  ) 

-211  FPFM  ATLL/.Z2  6X  t  1  9HTEST  TEMPERATURE  :  .  I3.19H  DFGRFFS  FAHRFNUF1T  .4X.7 
1HDATE  :  ,5A4) 

WRITE  I  6, 21 2) 

-21  2 _ FORMAT  (  /  ///BOX,  OUST  [FFNCSS  ,2  IX,  4HUN I T /5 6 X , 6 H ST 3 ES S  ,  6X  ,6HS.TR  A  I  N,  8X,_ 

16 HOF  MI  X ,PX , 4HT IMF, 9X  ,  6HW  E  I GH  T  /  57 X  , 5H( P  S I ) , 6X , 7H(  IN/IN)  , 7X, 5H( PSI) 

2 , 7X , 9H ( SCCONOS) ,5X,3H(P.C.F.)) 

_ WP  ITE  (6 , 21  3  )  XPAP5,  XBAR6.  IXBAR7  ,  XBAR8,  XBAR9 _ ; _ 


213  FORMAT! //23X  ,4HME AN  ,26X ,F9.2 , 3X ,F10 .6 , 5X , 18 ,3X, F10 . I, 7X ,F7.2/  ) 

WR ITE! 6 , 214)  S I  CM  A 5 , S I GM A 6 , I S CM A 7 , S  I  GM A8 , S I GM A 9 

214  FORMAT  (23X  ,  I  BUST AND ARD  OF  V I  A  T  I  0  N  ,  l  2X  ,  F  9 . 2 , 3  X ,  F  1  0 . 6 ,  5X  ,  I  8 , 3  X  ,  F  1Q_.  \  , 
17X, F7.2/) 

WP  ITE! 6, 21 5)  VAR5, VAR6 , VAR  7, VAR  8, VAR9 


c 

c. 

'  '  <-  ■>  h  1  >  LC'MI  ur  ¥  j  n  i  jj-ijN _ 1  a.  { 

1  ,  3  X , F  1 0 . 2 , 7X.F7.2///) 

c 

SUBROUTINE  ’STATS'  IS  CALLED  TO  DETERMINE  THE  MEAN  ,  STANDARD 

c 

DEVIATION  ,  AND  COEFFICIENT  OF  VARIATION  OF  EACH  PARAMETER 

c 

BASED  UPON  FAILURE  CONDITIONS  . 

ISN 

00  80 

c 

c 

CALL  S T ATS IN,FLSTRS,X BARI , S I GMA 1 , VAR l ) 

ISN 

"081 

CALL  STATS ( N ,FL STRN , X8AR2 ,S IGMA2, VAR2) 

ISN 

0082 

CALL  STATS ( N,FLSTFF , XB AR 3 , S I G MA 3 , V AR 3 ) 

I  SN 

0083 

CALL  STATS ( N  ,FLT IMF. XBAR4.S IGMA4. VAR41 

ISN 

00  84 

IXRAR3=XBAP.3 

ISN 

0085 

I S  GMA  3  =  S I GMA  3 

ISN 

0086 

I I  1  =  1 1  1  +  1 

ISN 

0087 

CALL  GRAPH! BUF, STIFF  , STRAI N , L , I X , T I TLE , SAMPLE ,  I , F TE MP , MAX S TR , N , I I 

1,1  II,  X  B  AR  3 , X  B A  R  2 , M  A  X  S  T  F  ) 

I  SN 

0088 

WRITE (6, 2 50)  (TITLE(K),K=1  ,20)  , NN.NN.N, !(SAMPLE(I,K),K=1.3),f  =  l.N) 

ISN 

0089 

250 

FORMAT !  1 H 1  //////  2  0  X ,  2  0  A  4 / / 9 1 X ,  5HPASE  ,  l 2 , 4H  OF  ,I2///29X, 

159HMEANS  ,  STANDARD  DEVIATIONS  ,  AND  COEFFICIENTS  OF  V AR I  A T I  ON/ 57X 

2 . 3HE0R/44X, 4HTHE  , I2.23H  SAMPLES  IN  THIS  S E R I ES / 44X , 2 9H B A SED  UPON 

3FAILURE  C0NDITI0NS///15X, 17H SAMPLE  NUMBERS  :  ,  3  X , 8 ( 3  A4 ) ) 

I  SN 

0090 

WRITE (.6, 251 )  TEMP,  {  DATE  C K )  ,K=1 , 5) 

ISN 

0091 

25  1 

FORMAT! //26X , 19HTEST  TEMPERATURE  :  ,I3,19H  DEGREES  F AHR ENHE I T ♦  4X •  7 

1IIDATE  :  ,  5 A 4  ) 

I  SN 

0092 

WRITE (6 ,252) 

ISN 

0093 

252 

FORMAT!  //// 8  OX,  9H  ST  IFF  NESS  ,2  1  X  ,  4HUN  I  T  /  56  X  ,  6  H  ST  P.  ESS,6X,6HSTRAIN,8X, 

16  HOF  MI  X ,eX , 4HT IME, 9X, 6HWE IGHT/57X,5H(PSI ) »  6X  ,  7H  (  IN/IN)  ,7X,5H(PSI ) 

2, 7X, 9H ( SECONDS) , 5X, 8H( P.C.F. ) ) 

ISN 

00  9  4 

WR  ITE  (  6 , 25  3  )  XBAP  1  ,  XBAR2,  IX  BARB  ,XBAR4,  XBAP. 0 

ISN 

0095 

253 

FORMAT ( //23X  , 4 HUE AN , 26 X , F9 . 2 , 3X , F 1 0 . 6 , 5X ,  18,3X,F10.1,7X,F7.2/J 

ISN 

0006 

WR  I  TE ( 6 ,25  4  )  SIGMA1 , S I GM A 2 , I SGM A3 , S l GMA4 , S I GM A9 

I  S  N 

0097 

254 

FORMAT ( 23X ,  1 8HST AND ARD  DEVIATI0N,12X,F9.2,3X,F10.6,5X,I8,3X,F10.l, 

17X , F7 . 2/ ) 

I  S  N 

0098 

WR  ITE! 6,255)  VAR1, VAR 2 , V A R3 , V AR 4 , VAR9 

ISN 

0099 

255 

FORMAT  (  23X  , 28  H  COEFFICIENT  OF  VARIATION  (  % )  ,  2X  ,  F  9 . 2 , 3  X  ,  F  ft.  2  ,  7  X  .  F  8 . 2 

ISN 

ISM 

ISN 

ISN 


0100 

0101 


0102 

0103 


1 »3X,F10.2»7X,F7.2///) 
300  CONTINUE 

_ CALL  PLOT (0.0, 0 .0 ,999) 


WRITE! 6,310) 
310  FORMAT <  1H l ) 


' 


■ 
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ISM  0104  STOP 

ISN  0105  END 

i  DC  ON  S  FOR  EXTERNAL  REFERENCES 

******  fmd  OF  COMPILATION  ****** 

ISN  0002  SUBROUTINE  M AX ( L , ARRAY  ,  NNN ) 

ISN  0003  DIMENSION  ARRAY (50) 

ISN  0004  AM AX= ARRAY ( 1 ) 

ISN  0005  NNN=1 

ISN  0006  DO  20  I=2,L 

ISN  0007  I  F I  APR  AY  (  I .)  .GT.  AMAX)  GO  TO  10 

ISN  0009  GO  TO  20 

ISN  0010  10  NNN= I 

ISN  00  11  AM  A  X  =  ARR  AY {  I  ) 

ISN  0012  20  CONTINUE 

ISN  0013  RETURN 

ISN  0014  END 

*****  END  OF  COMPILATION  ****** 

ISN  0002  SUBROUTINE  STAT S ( N , X , XB AR , S IGMA , VAR » 

ISN  0003  DIMENSION  X { 50) 

ISN  0004  S  U  M  =  0 .  . .  ...  .  _ 

ISN  0005  DC  1C  1=1, N 

ISN  0006  SUM=SUM+X( I ) 

ISN  no07  10  CONTINUE  _  _  _ _ .  ...  _  ....  _ 

ISN  0008  AN =N 

ISN  0009  BN=N-1 

ISN  0010  X 13  AR=  SUM  /  AN  . . .  .  -  -  ... 

ISN  0011  SMSQDF=0. 

ISN  0012  DO  2C  1=1, M 

ISN  0013  SMSQDF  =  SMSQDFM  XI  I  1-XBAR)  **2  . _ _ _  _ 

ISN  0014  20  CONTINUE 

ISN  0015  SIGMA=( S MS QD F / BN ) * * 0 . 5 

ISN  0016  VA  R  -  S I  GM  A  /  XB  Ap  *  1 00 .  ....  _.  .  .  - 

ISN  0° 1 7  RETURN 

ISN  0018  END 

*  *  *  *  *  end  OF  COMPILATION  ****** 

' 

■ 
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ISN  0002  SUP.  PHUT  INF  GRAPH! BUF, STIFF  ,  S  TR  A  I  N  ,  L  ,  I  X  ,  T  I  T  LE  ,  S  AM  PL  E  .  I  ♦  FT  E  MP  ,  M  AXS  T 

1R,N,  II  ,  I  II  ,  XBAR3,XBAR2,MAXSTF) 

ISR  0003  DI  PENSION  S  T I  TF  (  2  2  )  ,  S  TR  A  I  N  (  2  2 )  .  T  I  Tl  F  (  2  0  )  ,  S  AMPI  F  (  2  0 . 3  I  .  T  T  Tl  F 1  f  1  0  )  . 

1  T I T  L  G  2  (  10)  ,  S  AM ( 3  ) 

ISN  0004  .  RFAL  HA XST R , MAX S TF 

I  SN  0005  ST  IFFM  +1  I  =0  .n 

ISN  0006  STP.AINlL  +  l  )=0.0 

ISN  0007  STIFH L  +  2) =MAXSTF/7.0 

ISN  0r0?  STRAIN!!  +■?  )=MAXSTR *0.0001/5.0 

IS’  SCO?  r ' f I  .  GT .  1)  GO  TO  50 

is*.  "Oil  CALL  P.ECTtO.OiO. 0,11:0,8.5,0. 0,3) 

ISN  0012  CAM  PI  nT(2.  5,2. 5,-3) 

ISN  0013  CALL  A X I S I  0 . 0  , 0 . C , ' STR A I N  .  IN/ IN’ ,-14,5.0,0.0, STRAIN(L*1) *  STRAIN! 

1L+2) , 10.0) 

ISN  0014  CALL  AX  I  S  (  0  .  Or,  0 . 0  ,  ‘  S  T  I  F  FN F S S  OF  MIX  ,  P  S  I  •  ,  22 , 7 , 0 , 90 . 0  ,  ST  I  F  F  (  L  +  1  )  . 

1ST  IFF  C  L+2)  ,20.0) 

ISN  0015  CALL  PL OT ( 5 . 0 , 0 . 0 , 3 ) 

ISN  0016  CALL  PI  OT (5. 0,7. 0.2) 

ISN  0017  CALL  PLOTtO. 0,7. 0,2 ) 

ISN  0018  X  =  0.0 

IS1’  0  01°  00  1  J  =  1  . 6 

ISN  0020  CALL  SYMBOL! X »  7 . 0 , 0 . 10 , 13,0. 0,-1) 

ISN  0021  X=X+l.O 

ISN  0^22  1  CONTINUE 

ISV  0023  Y  =  0  . C 

ISN  0024  DO  2  J=1  ,  15 

ISN  0025  CALL  S  Y  •■'BOI  (  5 . 0  .  Y  .  0 . 1  0 , 1  3 . 90  -  0  . -1  ) 

ISN  0026  Y= Y+0 . 5 

ISN  0027  2  CONTINUE 

ISN  00  2  8  DO  3  J  =  1  , 1 0 

ISM  0029  TITLE 1(J)=TITLF(J) 

I  SN  00  30  3  CQNT  DUE 

ISN  0031  K- C 

ISN  00  32  DC  5  J= 1 1  ,  2C 

ISN  0033  K  =  K  + 1 

ISN  0034  T I TLF2 ( K ) =T I TLE ( J) 

ISN  0035  5  CONTINUE 

ISN  0036'  CALL  SYMBOL (0.5  ,£.8, C. 12, TITLE1 ,0.0,40) 

ISN  0^37  CALL  SYNBDt  10.5,6.6,0. 12, TITLE2, 0.0, 40) 

ISN  0038  CALL  SYMBOL (  1  .5 , 6.4 ,0 .  10,  19HTFST  TEMPERATURE  :  ,0.0,19) 

ISN  0039  CALL  NUMBE R ( 3 . 4 , 6 . 4 , 0 . 1 0 , FTEMP , 0 . 0 , - 1 ) 

ISN  0040  CALL  SYMBOL! 3. 7, 6. 4,0. 10, 1HF, 0.0,1) 

ISN  0041  CALL  S Y M DDL ( 2 . 9 , 5 . 5 , 0 . 1 2 , 6HL E GEND , 0 . 0 , 6 ) 

ISN  0042  CALL  SYMBDL!  2.2  ,5.25 , 0. 10, 24H SYMBOL  SAMPLE  NUMBER, 0 .0, 24 ) 

ISN  0043  Y  =  5 . 0 

ISN  0044  50  CONTINUE 

ISN  0045  IF! I I  .EQ.  Ill)  GO  TO  100 

ISN  0047  DO  7  J=l,3 

ISN  0048  SAM!  J)  ^SAMPLE!  I  ,.J) 

ISM  0049  7  CONTINUE 

ISN  0050  ISYMBL=I+127 

ISN  0051  CALL  L I NE( STRA I N , ST  IFF ,L , l » 1 , I SYMBL ) 

ISN  0052  CALL  SYMBOL ( 2 . 5 , Y , 0. 10 , I SYMBL , 0 . 0 , - l ) 

ISN  0053  CAt L  SYMBOL { 3 .4 , Y , 0 . 10 , S AM , 0 .0 , 1 2 ) 

ISN  0054  Y=  Y-0 • 2 

I 


. 


. 

; 


I  S'.!  0^55  RETURN 

I  SR'  0056  100  CONTINUE 

ISN  0057  CALL  SYMBOL!  1 .7, 6.0,0.  10.24HAVE  FAILURE  STIFFNESS  =  ,0.0,24) 

ISM  0053  CALL  NUMBFP (4.1  ,6.0,0. 10.  XBAR3,  0.0. -1  ) 

ISM  0050  CALL  SYMBOL!  2. 0,5. 8,0.  10 , 2  Hi  AVE  FAILURE  STRAIN  =  ,0.0,21) 

ISN  0060  CALL  NUMBER ( 4 . 1 , 5 .8 , 0 . 10 , XBAR2, 0.0 , 5) 

ISN  0061  IF! I 1/2*2  .EQ.  II)  GO  TO  121 

ISN  0063  CALL  PLOT ! - 1 0 . 0 , -  1 2 . 0, -3 ) 

ISN  0064  GO  TO  122 

ISM  0065  121  CALL  PL OT ! 8 . 0 , -3 . 0 , -3 ) 

ISN  0066  122  CONTINUE 

ISN  0067  RETURN 

ISN  0068  END 

******  fno  pp  COMPILATION  ****** 

COMPUTER  ANALYSIS  WITH  STRESS  -  STRAIN  AND 
STIFFNESS  -  STRAIN  PLOTS 
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isri  0  3  02 


ISN  C0C3 
ISM  0004 
ISM  0005 
ISM  000 b 
ISM  0007 


ISN  0008 
ISM  0009 
ISN  0010 


ISN  0013 
ISM  0014 


DIMENSION  TITLE!  20)  ,  0  A  T  E  (  5)  ,  T  I  M  E  (  2  0  )  ,  1 0  A  D  !  2  0  )  ,STRAIN(22  )  ,STRESS(22 
1  )  ,  FLSTP.  S  (2  0)  ,  FL  S  T  R  N  (  2  3  )  »FLT  IMF!  20)  »S  AMPLE!  20,3)  ,  S  T  I  F  F  (  20)  »FLSTFF(  2 
_Z-0-LlCJLS  TPS!  20  )  .ORSTRN!  ?C)  ,  ORSTFF  (?0)  .OPTTMF!  73  1  .UNirWTtPQI  .  RUF  I 


3)  rTITLEl! 1C )  ,  T  I  T  L  E  2  ( 10) r  S  A  M { 3 ) 
INTEGER  TIME  ,  S AMPL E , F A  I L T M , T E MP 

_ R  E_  A  L  _  L  C  A  D,  M  A  X  S  T  R  ,  M  A  X  S  TF _ 

CALL  PLOTS ( BUF, 3192 ) 

IX  =  0 
I  I  I  =C 


C 

c 

JC_ 

c 

c 

c 


M  =  THE  NUMBER  OF  SETS  OF  SAMPLES  TO 
THIS  RUN  . 


PE  PROCESSED  DURING 


READ! 5 , 10 ) M 
10  FORMAT ( I  5) 

DO  300  I  1  =  1  . M _ 


C 

C 

c 


c 

c 


TITLE  =  A  DESCRIPTION  (80  COLUMNS  OR  LESS)  TO  APPFAR  AT  THE 
TOP  OF  EACH  PAGE  TO  IDENTIFY  FACH  PARTICULAR  SET  . 


ISN 

0011 

READ! 5, 20)  ( T  ITLE(K)  ,K=1 ,20) 

ISN 

0012 

c 

20  FORMAT (20A4) 

c 

c 

TEMP  =  THE  TEST  TEMPERATURE  , 

IN  DEGREES 

FAHRENHEIT  ,  AT 

c 

WHICH  THIS  SELL  WAS  RUN 

« 

C 

C 

c 


RATE 


THE  RATE  OF  LOAD  APPLICATION  IN  INCHES  PER  MINUTE 


C 

c 

c 


DATE  =  THE  DATE  AT  WHICH  THE  SAMPLES  WIT  UN  THIS  SET 
WERE  TESTED  . 


RE  AD (5,30) TEMP » RATE, (DATE(K)  ,K=1,5) 
30  FOFMAT ( 7X ,  I  3 , F 1 0 . 3 , 5A4) 


C 

C 

c 

c 

c 

c 


MAXSTR  = 


THE  MAX  I MUM  STRAIN  ,  IN  TEN  THOUSANDTHS  OF  AN  INCH  . 
TO  BE  EXPECTED  WITHIN  THIS  RUN  .  THIS  IS  REQUIRED  TO 
ESTABLISH  THE  CORRECT  SCALE  FOR  STRAIN  ON  THE  CALCOMP 
PLOTTE R  .  THIS  MAY  ANY  VALUE  ,  BUT  THE  FOLLOWING _ 


C 

c 


VALUES  WILL  PRODUCE  CONVENIENT  SCALES 


c 

c 

c 

MAXSTR 

10 

1  "  =  2 

TEN 

SCALE 

THOUSANDS 

OF 

AN 

I  NCH 

L. 

JJiCil 

c 

25 

1 "  =  5 

TEN 

THOUSANDS 

OP 

AN 

I  NCH 

/ 

I  NCH 

c 

50 

ln  =  10 

TEN 

THOUSANDS 

c 

AN 

I  NCH 

/ 

INCH 

c 

100 

1”  =  2  0 

TEN 

THOUSANDS 

c 

AN 

I  NCH 

/ 

I  NCH 

ISM  00 1 .5 - 

c 

c 

P  P  A  0  (  5 . 3  5  ) MAXSTR 

'■  , 


' 
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I S N  0016  35  FORMAT ( F10  .C  ) 

C 

c 

_ _ C _ -  MAXSTF 

C 
C 

_ _ c _ 

c 
c 

_ r _ 

c 
c 

_ c  _ 

c 
c 

_ c_ _ _ _ MAXSTF _ SCALF 

C 

C  700000  l"  =  1CCC00  PSI 

_ .C _ 1400000 _ 1"  =_ 20000C  PS  I 

C  2100000  1"  =  300000  PSI 

C  2800000  1"  =  4COOOO  PSI 

_ C _ 3  500000 _ 1"  =  500  OCO  PSI 

c 

C  ETC 

_ c _  _  _  _ ; 


C  7000000  1"  =  1000000  PSI 

c 

c 


ISN 

0017 

READ(5,R6)  MAXSTF 

I  SN 

0018 

c 

36 

FORMAT ( F1C .0 ) 

c 

c 

N  =  THE  NUMBER  OF  SAMPLES  WITHIN  THIS  SET  . 

c 

c 

ISN 

0010 

RE  AD ( 5 »  40 ) N 

ISN 

0020 

40 

FORMAT (15) 

ISN 

0021 

c 

DO  200  1  =  1  , N 

c 

C  SAMPLE  =  THE  IDENTIFICATION  NUMBER  OF  THE  SAMPLE  . 

C 


_ c _ T.'TAIR  =  THE  WEIGHT  ,  IN  GRAMS  ,  OF  THE  SAMPLE  IN  aTR  . 

C 

c  WTWAT  =  THE  WEIGHT  ,  IN  GRAMS  ,  OF  THE  SAMPLE  IN  WATER  . 

_ C 

c 
c 

_ c 

c 
c 

_ c 

c 

C  FAILTM  =  THE  TIME  TO  FAILURE  ,  IN  SECONDS  . 

- 

ISN  00  22  RE  ADI  5, 50)  ( S  AMPLE!  I  ,K)  , K=  1 ,3)  ,WTAI R  .WTWAT  ,0  I  AM, THICK.L, FAILTM 


DIAM  =  THE  AVERAGE  DIAMETER  OF  THE  SAMPLE  ,  IN  INCHES  . 

THICK  =  TH  E  AVER  AGE  T  H  I_C.KN.E_SS  OF  THE  S  AMPLE  *,  IN  INCHES  .  . 

L  =  THE  NUMBER  OF  LOAD  AND  STRAIN  POINTS  PICKED  OFF  THE 
L  0  A  D  AND  DEFORMATION  STRI  P  CHART  . 


— — —  J.  H  E  MAXJMUJi  S  '  IF.fNFSS  .OF  M_LX  ,  IN  PSI  .  TO  BE  .EXPECTED 
WITHIN  THIS  RUN  .  THIS  IS  REQUIRED  TO  ESTABLISH  A 
CONVENIENT  SCALE  FOR  STIFFNESS  ON  THE  CALCOMP 

_ PI. OTTER  .  THE  VALUE  CHCSFN  SHOULD  BE  DIVISAI3LE  BY  7 

SINCE  THE  PR E-QE TER M I  NED  LENGTH  OF  THE  Y-AXIS  IS 
EQUAL  TO  7  INCHES  .  THE  FOLLOWING  VALUES  WILL  PRODUCE 

_ CONVENIENT  $_CALE  S  (  MOTE:  IF  TH  E_M  AXIMUM  STIFFNESS  IS 

NOT  KNOWN  THE  BEST  PROCEDURE  TO  FOLLOW  IS  TO  CHOSE  A 
VERY  LARGE  VALUE  TO  ENSURE  THAT  ALL  STIFFNESS  VALUES 
_ WILL  FALL  WITHIN  THE  L I MITS  OF  THE  PLOT  )  . 


. 

A70 


ISN  6023 


c 

c 

c 


c 

c 


5C  FORMAT  (  3  A4  ,  F  8 . 3 , 3F  1 0 . 3 , 2  I  10  ) 


-SPGR  =  BULK  SPEC  [FIC  GRAVITY  OF  THF  F.  0  *-*  P  A  C  T  E  0  MIXTURE 


ISN  0024 


S  P_G  R  =  W  T  AIR  /  (  W  TA  I  R  -  W  T  W  AT  ) 


C 

c 

c 

c 

c 


I  SN 
ISN 
ISN 


0025 

0026 

0027 


ISN 
_I  SN 


0028 

0020 


I  SN 
ISN 


0030 
00  31 


ISN  0032 


I  SN 
I  SN 


0033 

0034 


I  SN 
ISN 


0035 

0036 


UNI  TNT  =  UNIT  WEIGHT  OF  TIIF  POMP  AC  TED  _M  I  XTJJRt 
IN  POUNDS  PER  CUBIC  FOOT  . 


UNITWTI I )=SPGR*62.4 

LM=N+1 

NN=N+2 


60 


WR  I  T  C  I  6 , 6  C  )  (TITLE(K)  ,  K  -  1  »20)  *  I  , MN ,  ( S AMPLE ( I , K )  , K= l , 3 ) , S PGR , UN  I T WT ( 

1  I ) ,01AM, THICK 

FORMAT ( 1U1 //////20X ,20A4//91X ,  5MPAGE  , 12, 4H  OF  , I2///44X, 1  7  H 

1SAMPLC  MUNOER  :  ,3A4//32X,  49H BULK  SPECIFIC  GRAVITY  OF  THE  COM 

2P  ACTED  MIXTURE  =  ,  F5 . 3//3SX, 14 HUN  IT  WEIGHT  =  ,F5.1,22H  POUNDS 

2  RTF  CUBIC  E  QQJ , _ /  /  4  5  X ,  I  I  H  D  I  A  M  E  TER  =  , _ F5.3.7H  INCHES  / 


7  1 


3/44X  ,  12HTII ICKNESS  =  ,  F5.3.7H  INCHES//) 

WR  I  TE ( 6 , 61  )  RATE , FA ILTM, TEMP ,  (DATE ( K)  ,K  =  1 , 5  ) 

_  FOP.M  AT  (  2  5X  ,  1  5HDETA  1 1.  S  OF  TE  ST //26X  .  18HR  AT  E  OF  LOADING  =  ,F5.3,  16H 
1  INCHES  /  MINUTE, 4X, 19HTIME  TO  FRACTURE  =  ,I3,8H  SFCONDS//26X , 19HT E 
2 S T  TEMPERATURE  :  ,I3,I9H  DEGREES  FAHRENHEIT ,4X,9HDATE  :  ,5A4///) 

WRITE (6, 62) 


6  2  FORMAT ( 60  X , 9  HST IFFNESS/30X, 6 H STRESS ,  9  X , 6 H STRAIN , 1 1 X , 6H0F  MIX,  I  OX, 
14HT IME/31X ,5H(PSI),9X,7H(  IN/IN)  , 1  OX  ,  5H ( PS  I )  ,9X,9H( SECONDS )///) 
_ DEN0M  =  3.14I5  9*DIAM*THICK _ ; _ _ _ 


C 

C 

C 


TIME  =  THE  TIME  ,  IN  SECONDS  ,  FROM  ZFRO  TIME  TO  EACH  OF  THE 


C 

c 

c 


POINTS  ON  THE  LOAD  AND  DEFORMATION  STRIP  CHART  . 


70 


READ  (5, 70)  (T  I'-'  E  (  K  )  ,K  =  1  ,L  ) 
FORMAT ( 5X, 1515) 


C 

C 

c 


LOAD  =  THE  LOAD  ,  IN  POUNDS  ,  AT  EACH  OF  THE  POINTS  ON  THE 
LOAD,  TRACE  . 


ISN 

0037 

C 

c 

RC AD( 5, 80 ) ( L OAD ( K ) , K=  1  ,  L  ) 

ISN 

0038 

c 

r. 

80 

FORMAT ( 5X , 15F5.0) 

C 

STRAIN  =  THE  STRAIN  ,  IN 

TEN  THOUSANDTHS  OF 

AN  INCH  /  INCH  , 

C 

r. 

FOR  EACH  OF  THE 

POINTS  ON  THE  DEFOR 

MAT  I  ON  TRACE'. 

ISN 

00  39 

C 

RE  AD ( 5 , 90 )  (STRAIN iK),K  =  l»L) 

I  S  N 

00  40 

90 

FORM  AT ( 5X,15F5.0) 

ISN 

0041 

DO  100  J=1,L 

ISN 

00  42 

STRESS ( J ) =2 . *LOAD( J ) /DENOM 

I 

■ 


■ 


A71 


ISN  0043  ST  P.AIN(J)=STRAIN(J)  *0.0001 

C 

C 

.  Q...  CONS  IOTP.  IMG  AVERAGE  TENSILF  STRESS  OVER  A  ONE  [MOM  GAUGF  1  FNGTH 

C  AND  ASSUMING  A  POISSONS  RATIO  OF  0.33  . 

C 

c 

ISN  0044  STIFF!  J ) =0 . 9 12  *  ST RE SS ! J ) / { S TR A I N (  J)*.5) 

ISN  0045  ISTIFF= STIFF  (J) 

ISN  00.46  ...  .  WRITE  (f  ,95)STRESSU),STRAIN!J),ISTIFF,TIME(J) 

ISN  0047  95  FORMAT! 26X , F  10. 2 ,6X , FI  0. 6 , 8X , 18 ,9X,  15/ ) 

ISN  0048  100  CONTINUE 

ISM  0049  FT EMP  =  TFMP 

C 

c 

- .  -  c  SUBROUTINE  'GP.APH1'  IS  CALLED  SO  THAT  VALUES  REQUIRFD  FOR 

C  THE  CALCOMP  PLOTTING  OF  STRESS  VERSUS  STRAIN  CAN  BE  STORED 

C  ON  THE  7  TRACK  PLOTTING  TAPE  . 

C 

C 

ISN  0050  CALL  GRAPH1 ( BUF , STRE SS , STRA I N , L , I X , T I TL E, S AMPLE ,  I , FT EMP , MAXS TR , N ,  I 

11,111  ,  X  BAR  1  ,  XBAR2  )’ 

C 

C 

C  SUBROUTINE  'GRAPH2*  IS  CALLED  SO  THAT  VALUES  REQUIRED  FOR 

C  THE  CALCOMP  PLOTTING  OF  STIFFNESS  VERSUS  STRAIN  CAN  BE  STORED 

C  ON  THE  7  TRACK  PLOTTING  TAPE  . 

C 

C 

ISN  CO  5 1  CALL  GR APH2 ( BUF , ST  I F F , S TR A  I N , L , I X , T I TLE , S AM PLE ,  I , FTE MP , MAX S TR , N , I I 

1,111,  XBAR3, XBAR2,MAXSTF) 

C 

C 

C  SUBROUTINE  'MAX'  IS  CALLED  TO  DETERMINE  THE  MAXIMUM  SECANT 

C  MODULUS  (  THAT  IS  THE  MAXIMUM  STIFFNESS  OF  MIX  FOR  EACH 

C  STRESS-STRAIN  DIAGRAM  )  . 

C. 

c 

ISN  005?  CALL  M AX ( L , ST  I F F , NNN ) 

ISN  0053  ORSTRS ( I )=STRESS(NNN) 

ISM  0054  ORSTRN! I )=STRAIN(NNN) 

ISN  0055  ORSTFF (  I  )  =  ST  IFFINNNJ 

ISN  0056  ORTIMEl I 1 = T I  HE ( NNN ) 

ISN  0057  FL  STRS! I )  =  STRESS(L) 

ISN  0058  FLSTRN! I )=STRAIN(L). 

ISN  0059  FLSTFF!  I  )  =  ST  IFF(L) 

ISN  0060  FL  T  I  ME  (  I  )  =  TIME(U 

ISN  0061  200  CONTINUE 

C 

c 

C  SUBROUTINE  'STATS'  IS  CALLED  TO  DETERMINE  THE  MEAN  ,  STANDARD 

C  DEVIATION  ,  AND  COEFFICIENT  OF  VARIATION  OF  EACH  PARAMETER 

C  BASED  UPON  INITIAL  CONDITIONS  (  THAT  IS  THE  MAXIMUM  SECANT 

C  MODULUS  OF  EACH  STR E SS - S TR A  I N  DIAGRAM  )  . 

C  . .  . .  .  . . 

C 

ISN  006?  CAI  L  STATSfN,  OP.STRS,  X  BAR  5 , S I GMA 5 , VAR  5 > 
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I  SN 
I  SN 
I  SN' 
I  SN 
I  S  N 
I  SN 
I  SN 
I  SN 


0063 
0064 
006S 
00  66 


0067 
0068 
00  69_ 
0070 


ISN  CO  7 1 
I SN  0072 

ISN  0073 


CALL  STATS (N  ,  OR STRN , X B AR6 , S I GWA6 , V AR6 ) 

CALL  STATS ( N  ,OPSTFF , XBAR7 , S IGMA 7, VAR7) 

CALL  STATS  (  N  ,0P.  T I  ME , XBARO , S I GMA3 , VAR8 ) 

_ CAU-._iLlA.IS  C  M  ,UM  I  TWT  .  XB  AR 9_»_S  1 GMA  Q . VAR9 ) _ _ _ 

I  X  B AR7=  X BAR7 
ISGMA7=S IGMA 7 

_ WRITE ( 6  f  2 1 C )  ( T I TLE( K )  ,K=l , 20) ,LN,NN »N » <  ( SAMPLE!  I ,K)  *K  =  L  ,3 )  , I =  1 ,N) 

210  FORMAT (1  HI //////20X ,  20  A4//91X , .  5 HP  AGE  ,I?,4H  OF  , I2///29X, 

159HMLANS  ,  STANDARD  DEVIATIONS  ,  AND  COEFFICIENTS  OF  V A R  I  A T  I  ON / 5 7 X 

_ 2 ,  3 1  IF  PR  /  44  X  ,  4  IITI  IE  ,  I2.23H  SAMPLES  IN  THIS  SFRIFS  '44X  .  2  9H  BASED  UPON 

3IN1TIAL  C0NDITI0MS///15X,17IISAMPLC  NUMBERS  :  ,3X,fil3A4)) 

WP I TC( 6, 21 l ) TEMP  ,  ( DATE ( K)  ,K=l  ,5  ) 

21 1  FORMAT ( / / 2  6  X ,  19HTEST  TEMPERATURE  :  ,I3,19H  CETvEFS  F i ENHt I  T  ,  4  X  ,  7 
1  HD ATE  :  , 5  A4 } 

WRITE (6, 21 2) 


’  ^  1  '‘Ml  w  /  /  t  1  f  /  r  i  i  11  1  ',1.  )  _>  f  C  L  A  t  -t  ML7  m  1  l  1  .)  D  A  j  O  n  )  1  M)A  tDllB  1  'A'UlIl  OAi 

16HCF  M IX , 8  X , 4  H  T IMF, 9  X , 6  H  W  F I G  H  T / 5  7  X , 5  H ( P  S I ),6X,7H(IN/IN),7X,5H(PSI) 

2 ,7X,9H( SECONDS ) , 5X, 8H( P.C.F. ) ) 

ISN 

0075 

WRITE (6 ,213)  XBAR5, X3AR6,  I  XB AR7 , XBAR8 , XBAR9 

ISN 

0076 

213 

FORMAT! //23X , 4HME AN , 26X , F 9. 2, 3X, FI  0.6 ,5X, 18 ,3X, F10. 1 , 7X ,F7.2/) 

I  SN 

0077 

WRITE ( 6 ,214)  S IGMA5 , SIGMA6, I SGMA7, S IGMA 8, SIGMA9 

ISN 

0078 

214 

FORMAT! 23X, 1 8H STANDARD  DEVIATION,12X,F9.2,3X,F10.6,5X,I8,3X.F10.1. 

1  7X , F 7 . 2 /  ) 

ISN 

0079 

WR  I  TE  (  6 , 21  5  )  VARS  ,.VAR5  ,  VAR  7  ,  VAR  8  ,  VAR9 

T  SN 

0080 

215 

F  C  P  M  AT ( 2  3  X ,  2  8H  0  0  E  F  F I C I  ENT  OF  VARIATION  (%) , 2X.F9.2,  3X.F8.2. 7X,F8.2 

C 

c 

1 , 3X  ,  F 1 0 . 2 , 7X.F7.2///) 

c 

SUBROUT  IMF  'STATS'  IS  CALLED  TO  DETERMINE  THc  MEAN  ,  STANDARD 

c 

DEVIATION  ,  AMO  COEFFICIENT  OF  VARIATION  OF  EACH  PARAMETER 

c 

BASED  UPON  FAILURE  CONDITIONS  . 

I  SN 

00  81 

c 

c 

CALL  STATS! N.FLSTRS.XBARl , S I G MA 1 , V AR l ) 

I  SM 

ISN 

I  SN_ 

ISN 

ISN 

ISN 


0°  8  2 
0083 
00  8  4 
0085 
0086 
0087 


CALL  STATS  (  N  ,f 
CALL  STATS  (_N  ,_F 
I X  BAR3  =  X  BAR3 
I SGMA3=S IGMA3 
I  I  I  =  I  I  I  + 1 


,  XB  AR3  ,S IGMA3, VAR3) 
,  XBAR4 , S I GMA4.VAR4) 


ISN  0088 


ISN  0080 


ISN  0090 


CALL  GRAPH  1  (  BlJF,  STRESS,  STRAIN,  L 
11,111,  XBAP  1  ,  XBAR2  ) 


1  ,  I  I  I  ,  X (3  A P ?■  ,  XBAP2  ,  MAXSTF  ) 
K  R  I  T  E ( 6 , 2  5  C )  ( TITLE! K)  ,  K  =  1 


20) 


IX, TITLE, SAMPLE, I , F TE MP , MAX STR , N , I 

TEMP,MAXSTR,N.  I  I. 
N  N , N  N , N , (  (SAMPLE!  I , K ,  ,  X  =  1  ,  3 )  ,  I  =  1  , N ) 


1  IN 

U  '  J  V  1 

1 59HME  AN  S  ,  STANDARD  DEVIATIONS  ,  ANO  COEFFICIENTS 
2,  3IIF0P./44X  ,4HTUE  ,12,2311  SAMPLES  IN  THIS  SERIFS/4- 
3FAILIJRE  COND  IT  IGNS///15X,  17HS  AMPLE  NUMBERS  :  ,  3X,3 

•  1  1  1  C  L-J  L  :  ..  . 

OF  VARIAT I0N/5TX 

X  ,29*194660  UPON 
i  3A4) ) 

ISN 

I  SN 

0092 

0093 

WRITE (6, 251) TEMP, ( DATE (K)  ,K  =  1,5) 

251  FORMAT ( //26X , 19HTEST  TEMPERATURE  :  ,I3,19H  DEGREES 
1HD ATE  :  , 5 A 4  ) 

r  A  ,-.R  ENHE  I  T  ,  4  X ,  7 

ISN 

0094 

WRITE! 6,252) 

I  SN 

0095 

252  FORMAT (//// 80X.9HST IFF  NESS ,2 IX, 4HUN I T / 5 6 X , 5H STR ESS 

, o  X , 4  H  S  T  R  A  IN, 8  X , 

16H0F  M IX ,LX , 4HT I  ME , 9 X , 6HW  E I GHT / 57X , 5H( PS  I ) , 6X , 7 HI  r 

N/IN) *7X,5HIPSI ) 

2  *  7  X ,  9H  (  SEC  Of;  OS  )  ,5X,8H  (P.C.F.)) 

ISN 

0096 

WRITE! 6,2  53)  XBAR  1 , XBAR2 ,  I XBAR3  ,  XBAR4 , XBAR9 

ISM 

I  SN 
ISM 

0097 

0098 

0099 

253  FORMAT(//?3X,4HMEAN»26X*F9.2»3X»F10.6,5X»I8»3X»F10.1,7X,F7.2/) 

WRITE (6 ,254)  SIGNAL , SIGMA?, I S CM  A3 , S I GMA4, S I GMA9 

25  4  FORMAT  (2  3X,  1  8FI  STANDARD  DE VI  AT  I  ON, 1 2X ,F9 .2 , 3X, F10. 6 , 5X , I 8, 3X , F10 . 1 , 
lIXjJE7t2/j 

I  SN 

0100 

WRITE! 6,255)  VARl ,  VAR2 ,  VAR3 ,  VAR4.VAR9 

I  SN 

0  101 

255  FORM  AT ( 23X ,2 QHCOEFF IC I  ENT  OF  VARIATION  ( %  ) , 2X , F 9 . 2 , 3 X, F 8 . 2 , 7 X , F 8 . 2 

— 

_ 1 , 3X , FI  0.2, 7X ,F7 ,2//  /  ) 

ISM 

010? 

30C  CONTINUE 

ISM 

0103 

CALL  PLDTIO. 0,0. 0,999) 

I  SN 

0104 

WRITE (6 ,310  ) 

ISM 

0105 

310  FORMAT ( LH1 ) 

I  SN 

0106 

STOP 

I  SN 

0 107 

END 

OCONS 

FOP  EXTERNAL  REFERENCES  "  ""  ~ 

*  if  * 

END  OF 

COMPILATION  ****** 

I  SN 

000? 

SUBROUTINE  MAXIL, ARRAY, NNN) 

ISM 

0003 

DIMENSION  ARRAY ( 50 ) 

I  SN 

OnP4 

AM  AX- ARR  AY (  I  ) 

I  SN 

0005 

NNN  =  1 

I  SN 

0006 

DO  20  1=2, L 

ISN 

0007 

I F ( ARR AY (  I  )  . GT  .  AMAX)  GO  TQ  10 

I  SN 

0009 

GO  TO  20 

ISM 

0010 

10  NNN  = I 

ISM 

noil 

AM  AX  =  ARP.  AY  (  I  ) 

I  SM 

0012 

20  CONTINUE 

ISN 

00  13 

RETURN 

I  SN 

0014 

END 

>Jc  #  A  #  £ 

FND  OF 

COMPILATION  ****** 

ISN 

0002 

SUBROUTINE  S TAT S ( N , X , X  BAR , S I GMA , V AR ) 

ISM 

0003 

D I  MENS  ION  X! 50) 

ISM 

0904 

S  U  M  =  0 .  . .  _ 

ISN 

0005 

DO  10  1=1, N 

ISN 

0006 

SUM  =  SUM  +  X(  I  ) 

I  SM 

00  07 

10  CONTINUE  ..  _  -  __  -  _  _  --  ----- 

ISN 

0008 

> 

II  ' 

~T  ! 

i 

ISN 

0009 

BN=N-1 

I  SN 

0010 

X  B  A  P  =  S  U  M  /  AM  - 

ISN 

00  11 

•  SMS0DF=0. 

ISN 

0012 

DO  20  1  =  1  »  N 

1  S  N 

0013 

SMSQDF=SMSQDF+(X ( I J-XBAR) **2 

I  SN 

0014 

20  CONTINUE 

ISN 

0015 

S I CMA= ( SMSQDF/BN )**0 .5 

I  SN 

0016 

VAR=SIGMA/X3AR*100.  . - 

ISN 

0017 

RETURN 

ISN 

0018 

END 

'****  & 

END  OF 

COMPILATION  ****** 

' 

■ 


A74 


ISN  0002 


ISM  0003 


SUGRQUT I Nfc  GRAPH  1 ( BUF , STRESS , STRA I N , L , IX, T I TLE, SAMPLE ,  I , FT EMP, MAXS 
1TP  ,N  ,  I  I  ,  I  I  I  ,  XBAP.  1 ,  XBAR2 ) 


ISM 
I  SN 
I  SM 
ISM 
I  SN 


00  C  4 
0005 
0006 
0007 

oooa 


1TITIF_2(  10)  ,5  AM  (3)  ,  PL  RE  SSI  23)  , PLRAINI23) 
PFAL  MAXSTR 

_ INTEGER  U,V,H 
PL  RF  SS ( l ) -0 . 0 
PLRAIMI 1 5=0.0 
U=L»1 _ 


ISN 
I  SN 
I  SM 
ISN 
ISN 
ISM 


oooo 
0010 
0011 _ 
0012 
0013 
00  14 


W=  1 

DO  300  V  =  2  ,  U 

_ £LET;  SSI  V)  =  STRE  SSI  H) 

PLRAIMI V)=STRA INI W) 
W=W  +  l 

300  COMTIMUF _ . 


ISN 
ISN 
ISN 
ISN 
ISN 
I  SN 


0015 

0016 

0.0.17 

0018 

0019 

0021 


PLRESS (U  +  l 5=0.0 
PLRAINIU  +  1 5  =0.0 

PLPCSS(U+?)=100.0 _ 

PLRA INI U  +  2 )= MAXSTR *0.0001/5.0 

IF  I  I  .GT .  1 5  GO  TO  5  0 

IF (  I  1/2*2  .FO.  I  I  5  GO  TO  21 


I  SN 
ISN 
I  SLL 
ISM 
ISN 


0923 

0024 

0025 

0026 

0027 


CALL  RECT (0.0, 0.0,22. 0,17. 0,0. 0,3) 
CALL  PL 0110.0,9.0,-3) 

JLL_C  CNTIN  U  E _ ’ _ 

CALL 


PLOT! 3.0,3 .0,-3) 

CALL  AXISIO. 0,0.0,' STRAIN 


,  IN/  IN' ,-14,5.0 ,0.0, PLRAINIU  +  1 ), PLRA INI 


ISN 

0028 

CALL  AX  IS! 0 .0,0 .0, • TENSILE  STRESS  ,  PS  I ' , 20 , 7 . 0 , 90 . 0 , PL RESS ( U  + 1 ) , P 
1LRESSIU+2) ,20.0) 

ISN 

0029 

CALL  PLOT  15.0,0.0,3) 

ISN 

0030 

CALL  PLOT (5. 0,7. 0,2) 

ISN 

0031 

CALL  PLOT! 0. 0,7. 0,2) 

ISN 

00  32 

X  =  o  .0 

I  SN 

0033 

DO  l  J=1 ,6 

ISN 

0034 

CALL  SYMBOL! X,7.0,0. 10, 13,0. 0,-1) 

I  SN 

0035 

X=X+1 .0 

ISN 

0036 

1 

CONTINUE 

ISN 

0037 

Y  =  0.0 

ISN 

0038 

DO  2  J  =  1?  1  5 

ISN 

0039 

CALL  SYMBOL! 5.0, Y,0. 10 , 13,90.0,-1 ) 

ISN 

00  40 

Y= Y+C . 5  v 

I  SN 

0041 

2 

continue 

ISN 

004? 

DO  3  J= l , 10 

I  S  N 

0043 

TITLEH  J  )  =T  I  TLE  (  J  ) 

I  SN 

0044 

3 

CONTINUE 

ISN 

0045 

K  =  C 

ISN 

0046 

DO  5  J= 1 l , 20 

ISN 

0047 

K  =  K  +  1 

ISN 

0048 

TITLE2I K)=TITLE( J) 

I  SN 

0049 

5 

CONTINUE 

ISN 

0050 

CALL  SYMBOL  I  0.5 ,6.8,0. 12, TITLE  1 ,0 .0 ,40 ) 

ISN 

0051 

CALL  SYMBOL  I  0.5, 6. 6,0. 12, TITLE 2, 0.0, 40) 

I  SM 

0052 

CALL  SYMBOL! 1 .5, 6. 4,0. 10, 19HTEST  TEMPERATURE  :  ,0.0,19) 

ISN 

0053 

CALL  NUMBER! 3.4,6.4,0.10,FTEMP,0.0,-1) 

ISN 

0054 

CALL  SYMBOL! 3. 7, 6. 4,0. 10,  1  HE, 0.0, 1) 

' 

■ 


A75 


I  SN 

0055 

IF  {PLRESSIU)  .GT.  400.)  GO  TO  30 

ISM 

0057 

V 1  -  6 . 0 

ISM 

0058 

Y2  =  5 . 8 

ISM 

0059 

Y  3  =  5 ,5 

ISM 

0060 

Y  4  =  5 . 2  5 

ISM 

0061 

Y5=5 . 0 

ISM 

0n62 

GO  TO  40 

ISM 

0063 

30 

CONTINUE 

ISM 

0064 

AM-N 

r  sn 

0065 

Y 1 =  AM*0 .2*1  .0 

ISM 

0066 

Y2=Yl-0.2 

ISM 

P°67 

Y3  =  Y  2-0 . 3 

I  SM 

0068 

Y 4  =  Y 3 - C  .25 

ISM 

0069 

Y5=Y4-0 . 25 

ISN 

0070 

40 

CONTINUE 

r  s  n 

0071 

CALL  SYMBOL!  2 . 9  ,  Y3 , 0 . 1  2 . 6 1 11  FG  FN  0 . 0 . 0 . 6  ) 

I  SM 

00  7  2 

CALL  SYMBOL! 2.2 ,Y4»0. 1 0 »  2 41! SYMBOL  SAMPLE  NO M8E R , 0 . 0 , 24 > 

ISM 

00  73 

Y  =  Y  5 

ISM 

0  0  74 

5 CL 

CONTINUE 

ISM 

0075 

IF! II  .EQ.  Ill)  GO  TO  100 

ISM 

0077 

00  7  J  =  1  ,3 

ISM 

007  e 

SAM! J)=SAMPLE( I.J) 

ISM 

00  79 

7 

CCNT  INUE 

I  SN 

0080 

I  S  Y  M  P  L  =  I  +  1  2  7 

ISM 

00  81 

CALL  L INE ( PLRAIN.PLRESS ,U,  1, l , I SYMBl ) 

ISM 

0082 

CALL  SYMBOI  ( 2 .5 , Y,0  :i0, ISYMBL ,0.0,-l> 

ISM 

0083 

CALL  S  YMBOL ( 3 . 4 , Y , 0 . 1 0  , SAM ,0.0,12) 

ISM 

0^84 

Y  =  Y  -  0 . 2 

I  SM 

0085 

RETURN 

I  SM 

00  8  6 

100 

CONTINUE 

I  SM 

00  8  7 

CA!  L  SYMBOL!  2.0, Yl.0. 10, 21HAVF  FAILURE 

STRESS  =  .0.0.21) 

ISM 

0088 

CALL  NUMBER! 4 . 1 , Y 1 , 0 . 1 0 , X B AR 1 , 0 .0 , 2  ) 

ISM' 

0089 

CALL  SYMBOL! 2.0, Y2, 0.10, 2LHAVE  FAILURE 

STRAIN  =  ,0.0,21) 

I  SM 

0090 

CALL  NUMBER! 4. 1 , Y2 , 0. 1 0 , XBAR2 ,0 .0 . 5 ) 

I  SM 

0091 

RETURN 

ISM 

0092 

END 

END  OF 

COMPILATION  ****** 

i 

, 

■ 
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I  S  N 

I  S  N 

0002 

0003 

SUBROUTINE  GRAPH2 ( B UF , S T I  FT , STR A  I N , L , I X , 7 1  TIE , S AMPL E , I , FTEMP , MAXS T 

1 R  »  N  .  1 1  ,  I  I  I  ,  XI3AR3,  XBAR2  ,  MAXSTF) 

DI MFNS I  ON  STIFF122)  »STRAIN(22)  .  T  I  T l  F  (  2  0  )  .SAMPI  FI20.3) .TfTI  F 1  ( i 01 . 

ISN 
_ I  SN 

0C04 

non? 

IT  I T L  6 2 (  10  )  ,S AMI  3  ) 

REAL  MAXSTR, MAXSTF 

ST  I  FF  (L  +  l)  =0.0 

ISN 

ISN 

I  S  M 

0006 

0C07 

onon 

STRAIN! L+l ) =0.0 

STIFF(L+2)=MAXSTF/7.0 

STRAIN ( It? )=MAXSTR*0.0C01  /  5,0 

ISN 

0009 

CALL  PLOT (7. 0,0. 0,-3) 

ISN 

0010 

IF (  I  .GT.  1 )  GO  TO  50 

ISN 

001? 

CALL  _ A X 1 S ( 0 . 0,0.0, 'STRAIN  ,  IN/ IN', -14, 5. 0,0.0, ST  RAIN!  L+l). STRAIN! 

1L  +  2)  ,  10  .0) 

ISN 

0013 

CALL  AXI S( 0. 0,0.0, ' STI FFNESS  OF  MIX  ,  PS  I '  , 22 , 7 . 0 , 90 . 0 , ST l F F ( L +1 )  , 

1STIFFIL+2) ,20.0) 

ISN 

0014 

CALL  PLOT (5.0,0 .0,3 ) 

ISN 

0015 

CALL  PLOT  I  5 . 0,7.0 ,2 ) 

I  SO 

0016 

CALL  PL0T(0. 0.7. 0.2) 

'  ISN 

0017 

X  =  0.0 

ISN 

0018 

DO  1  J= 1 , 6 

I  SN 

0019 

CALL  SYMBOL ( X, 7. 0.0. 10.13.0.0.-1) 

ISN 

0020 

X=X+ 1 .0 

ISN 

0021 

1  CONTINUE 

ISN 

002? 

Y  =  0  .  C 

I  SN 

0023 

DO  2  J  -  1  ,  1  5 

I  SN 

0024 

CALL  SYMBOL (5. 0, Y, 0.10, 13, 90. 0,-1) 

ISN 

0025 

Y  =  Y  +0 . 5 

ISN 

0026 

2  CONTINUE 

ISN 

0027 

DP  3  J=l , 10 

I  SN 

0  0  ?  R 

TITLEl ( J)-TITLF(J) 

I  SN 

0029 

3  CONTINUE 

ISN 

9030 

K  =  0 

ISN 

0031 

DO  5  J  =  1  1  .20 

ISN 

003? 

K=  K  + 1 

ISN 

0033 

TITLE2( K)=T  I  TLF( J) 

I  SN 

0034 

5  CONTINUE 

ISN 

0035 

CALL  SYMBOL ( 0. 5,6. 8,0. 12, TITLEl ,0.0 ,40) 

I  SN 

0036 

CALL  SYMBOL (0.5,6.6,0.12, TITLE2, 0.0, 40) 

ISN 

0037 

CALL  SYMBOL!  1.5, 6. 4 ,0. 10, 19HTFST  TEMPERATURE  :  ,0.0,19) 

ISN 

0038 

CALL  NUMBER ( 3.4,6.4,0.10,FTEMP,0.0,-1) 

ISN 

0039 

CALL  SYMBOL (  3.7,6.4,0.10,  1HF, 0.0, 1) 

I  SN 

0  0  40 

CALL  SYMBOL  (  2.9,5.5,0.12, 6  H  LEGEND, 0  .Q.t6  ) 

I  SN 

0041 

CALL  SYMBOL! 2.2, 5.25, 0. 10, 24H SYMBOL  SAMPLE  NUMB ER , 0 . 0 , 2 4 ) 

ISN 

9042 

Y  =  5 .0 

ISN 

0043 

50  CONTINUE 

ISN 

0044 

IF  (  I  I  .EQ.  Ill)  GO  TO  100 

I  SN 

0046 

DO  7  J  =  1  ,3 

ISN 

0047 

SAM! J)=SAMPLE( I , J) 

ISN 

0048 

7  CONTINUE 

ISN 

0049 

ISYMBL-I+127 

I  SN 

0050 

C A 1  L  L INE! STRAIN, ST  IFF ,L , 1 , 1 , ISYMBL ) 

ISM 

9051 

CALL  SYMBOL (2.5,Y,0.10,ISYMBL,0.0,-1) 

ISN 

0052 

CALL  SYMBOL ( 3.4, Y,0. 10 , SAM, 0.0, 12 ) 

ISN 

0053 

Y  =  Y-C  .  2 

ISN 

0054 

CALL  PLOT (-7  .0,0. 0,-3) 

■ 
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ISN  0055  RETURN 

I  5 r  0056  100  CONTINUE 

ISr'  0057  CALL  SYMBOL  (  1 .7 ,  6 . 0  ,  C  .  10, 24  HAVE  FAILURE  STIFFNESS  =  ,0.0,24) 

_  ISN  0058  ...  CALL  NU  M  3  E  R (  4.  1  ,6.0,0.  10,  XBART.0.0  .-I  ) 

ISN  0959  CALL  SYMBOL  (2. 0,5. 8,0.  10, 21  hi  A  VE  FAILURE  STRAIN  =  ,0.0,21) 

ISN  0060  CALL  NUMBER ( 4.1 >5.8,0.10,XBAR2,0.0,5) 

ISM  0061  \y- ry+! 

ISN  0062  IF  (IX  .EQ.  2)  CO  TO  10 

ISN  0064  CALL  PL OT I -2  .  5 , 8 . 5 , -3 ) 

ISN  0065  RETURN 

ISN  0066  10  CALL  PLOT  I 6 . 0 , - 13 . 5 , -3 ) 

ISN  0067  IX=0 

ISN  0068  RETURN 

ISN  0069  END 

*****  EUO  OF  COMPILATION  ****** 

' 


APPENDIX  B 


COMPUTER  PLOTS  OF  TENSILE  STRESS  VS.  STRAIN 
AND  STIFFNESS  OF  MIX  VS.  STRAIN  FOR 
HIGHWAY  2-D-2/1  S  2/2  CORES 


AT  VARIOUS  TEMPERATURES  AND  AGES 


TENSILE  STRESS  ,  PSI  TENSILE  STRESS  .  PSI 

JKOO  ICO. CO  200.00  300. CO  |  400.00  |  SOO.OQ  600.00  700.00  J.00  100.00  |  200.00  ^  SCO.CO  400.00  500.00  600.00  700.00 


B1 


-+- 


-+- 


~h 


-4- 


HWT  2-0-2/ 1  4  2/2  SUPPLY' NO.  1  (200/300) 
COnES  RUGUST  1966  RT  CONSTRUCT  ION 
TEST  TEMPERATURE  t  20  F 


fiVE  FAILURE  STRESS 
FIVE  FRIEURE  STPWIN 


240. 26 
0.00133 


SYMBOL 


LEGEND 

SPITTLE  NUMBER 
C  47-B 

c  i-a 

C  23-fl 
C  27-fl 

c  n-a 

C  17-8 


0.10  0.20  0.30  0.40 

STRAIN  ,  IN/IN  (X10-?) 


0.50 


Oio 

o 


Qt 

in' 

cn 

LU 


-4- 


-+- 


-+- 


HNY  2-D-2/1  4  2/2  SUPPLY  NO.  1  (200/300) 
CORES  AUGUST  1966  RT  CONSTRUCTION 

TEST  TEMPERATURE  i  20  F 

RYE  FAILURE  STIFFNESS  =■  2S1590 

AYE  FAILURE  STRPIN  =.  0.00193 

LEGENO 


STM80L 


SAMPLE  NUH8ER 
C  47-3 
C  1-3 
C  23-A 
C  27-A 
C  11-3 
C  17-3 


Tt.OO 


0.10  0.20  0.30  0.40 

STRAIN  .  IN/IN  (X10-2) 


0.50 


FIG  .  B1 ( a) 


FIG.  Bl(b) 


-+- 


HWY  2-0-2/ 1  4  2/2  SUPPLY  NO.  1  (200/300) 
CORES  RUGUST  1966  RT  CONSTRUCTION 

TEST  TEMPERATURE  I  10  F 

AYE  FAILURE  STRESS  -  264.93 
AYE  FAILURE  STRAIN  »  0.00106 

LEGENO 

SYMBOL  SAMPLE  NUMBER 
C  47-A 
C  13-A 
C  7-3 
C  15-P 
C  49  3 


.00  0.10  0.20  0.30  0.40  0,50 

STRAIN  .  IN/IN  (X10-2) 


coo 
O-  a 


fef, 

to' 

CD 

UJ 


IT) 


-4- 


13.00 


HNY  2-D-2/1  4 
CORES  RUGUST 
TEST 


-4- 


-T- 


-t- 


2/2  SUPPLY  NO.  1  (200/300) 
1966  RT  CONSTRUCTION 

TEMPERATURE  i  10  F 


AVE  FO'LURE  STIFFNESS  = 
AYE  FAILURE  STRAIN  » 


4B0837 

0.00106 


SYMBOL 


LEGENO 

SAMPLE  NUMBER 
C  47-A 
C  13-A 
C  7-3 
C  1 5—0 
C  49-3 


0.10  0.20  0.30  0.40 

STRAIN  .  IN/ 1  NIX  1 0-2) 


o.so 


FIG.  Bl(c) 


FIG.  Bl(d) 


FIGURE  Bl:  Supply  No.  1  -  At  Construction 


. 


B  2 


85.00 


0.05  0.10  0.15  0. 20 

STRAIN  ,  IN/ IN  (XI O'2 ) 


0.25 


FIG.  B2 ( a) 


8 


cr>8 

CL_o 


x: 

szr 


CO 

CO 

UJ 

fc8 

t— »o 

CO 


-4- 


-4- 


-4- 


8).  oo 


HNY  2-0-2/ 1  4  2/2  SUPPLY  NO.  1  (200/300) 
CORES  AUGUST  1966  AT  CONSTRUCTION 
TEST  TEMPEflfiTURE  i  0  F 


PVE  FAILURE  STIFFNESS  « 
five  FRiLune  strain  - 
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FIGURE  B2 :  Supply  No. 
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HHT  2-D-2/1  <,  2/2  SUPPLY  NO.  2  (200/300) 
CORES  AUGUST  1966  RT  CONSTRUCTION 

TEST  TEMPERATURE  i  20  F 
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FIGURE  B 3 :  Supply  No.  2  -  At  Construction 
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HUT  2-0-2/1  4  2/2  SUPPLY  MO.  2  (200/3001 
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TEST  TEMPERATURE  i  0  f 
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FIG.  B4 ( a) 
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FIG.  B  4  (  d ) 


FIGURE  B4 :  Supply  Mo.  2  -  At  Construction 
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FIG.  B5 ( a) 


FIG.  B  5 ( b ) 
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KWT  2-0-2/1  «,  2/2  SUPPLf  NO.  3  (200/300) 
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TEST  TEMPERATURE  I  10  F 
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HHT  2-0-2/ 1  6  2/2  SUPPLY  NO.  3  (200/300) 
CORES  AUGUST  1966  AT  CONSTRUCTION 
TEST  TEMPERATURE  I  10  F 
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FIG.  B5(c) 


FIG.  B5(d) 


FIGURE  B5 :  Supply  No.  3 


At  Construction 
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STRAIN  .  IN/IN  (XI O’2) 
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FIGURE  B6 :  Supply  No.  3  -  At  Construction 
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TEST  TEMPERATURE  i  20  F 
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HWY  2-0-2/1  4  2/2  SUPPLY  NO.  1  (200/300) 
CORES  FEBRUARY  1960  10  MONTHS  SERVICE 

TEST  TEMPERATURE  i  20  F 
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FIG.  B7  (a) 


FIG.  B  7 ( b ) 
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FIG.’  B  7  (  c  ) 
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FIG.  B  7 ( d ) 


FIGURE  B 7 
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HNY  2-0-2/1  4  2/2  SUPPLY  NO.  1  (200/300) 
CORES  FEBRUARY  1963  18  MONTHS  SERVICE 
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RHY  2-0-2/ 1  4  2/2  SUPPLY  NCI.  1  (200/300) 
CORES  FEBRUARY  1968  13  RCNTHS  SERVICE 

TEST  TDiPEftHTUfE  |  -10  F 
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FIG.  B9(a) 
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HWY  2-0-2/ 1  4  2/2  SUPPLY  NO.  2  (200/300) 
CORES  FEBRUARY  1968  10  MONTHS  SERVICE 

TEST  TEMPERATURE  i  10  F 
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HUY  2-0-2/ 1  4  2/2  SUPPLY  NO.  2  (200/300) 
CORES  FEBRURRY  I960  10  MONTHS  SERVICE 

TEST  TEMPERATURE  i  10  F 

RYE  FRILURE  STIFFNESS  -  683S7S 

RYE  FRILURE  STRRIN  -  0.00128 
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FIG.  B9(c) 
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FIGURE  B9 :  Supply  No.  2 
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FIG.  B10(b) 
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COPES  FEBRUARY  1963  10  MONTHS  SERVICE 

TEST  TEMPERATURE  i  20  F 
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HMT  2-0-2/ 1  «,  2/2  SUPPLY  NO.  3  (200/300) 
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HUY  2-D-2/1  L  2/2  SUPPLY  NO.  3  (200/300) 
CORES  FEBRUARY  1960  10  MONTHS  SERVICE 

TEST  TEMPERATURE  I  10  F 
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FIG,  Bll(c) 


FIG.  Bll(d) 


FIGURE  Bll :  Supply  No.  3 
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HUT  2-0-2/ 1  2/2  SUPPLY' NO.  3  (200/300) 

cones  FEBRUARY  1968  10  MONTHS  SERVICE 

TEST  TEMPERBTUfT£  i  0  F 


BYE  FBI LURE  STRESS 
RYE  FB 1 LURE  3TF5BIN 


303.58 

0.00031 


LEGENO 

STM80L  STSPLE  NUH6EB 
235 35  L.C. 
Z3S39  U.C. 

23531  L.C. 

23532  U.C. 
23530  L.C. 


0.02  0.04  0.05  0.03 

STRAIN  .  IN/IN  (X10-2) 


o.io 


X 
* — 1 


o«.. 

<M 

CO 

CO 

UJ 

2: 

l*-Q 

>—*  o 
i-p: 
07 


-+- 


‘Tj.oo 


HHY  2-0-2/ 1  <  2/2  SUPPLY  NO.  3  (200/300) 
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FIG.  B1 2 ( a) 


FIG.  B1 2 (b ) 
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FIG.  B1 2 ( c ) 
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BUY  2-0-2/ 1  2/2  SUPPLY  NO.  3  (200/300) 

CORES  FEBRUARY  1963  18  MONTHS  SERVICE 

TEST  TEMPERATURE  1  -10  F 
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FIGURE  B12 :  Supply  No.  3 
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FIGURE  B13:  Supply  No.  1  -  32  Months  Service 
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HWT  2-0-2/ 1  4  2/2  SUPPLY  NO.  2  (200/300) 
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FIG.  B14 (a) 


FIG,  B 1 4 ( b ) 
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FIG,  B14 (d) 


FIGURE  B14 
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HMT  2-0-2/ 1  s,  2/2  SUPPLY  NO.  2  (200/300) 
CODES  APRIL  1969  32  MONTHS  SERVICE 

TEST  TEMPERA TUfiE  i  20  F 


RYE  FAILURE  ST1FFNES3  ■ 
RYE  FAILURE  STFWN 


128464 

0.00372 


SYMBOL 


LEGEND 

SAMPLE  NUMBER 
23  U.C. 
23  U.C. 
29  L.C. 
22  U.C. 

22  L.C. 

23  L.C. 
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STRAIN  .  IN/ IN  (XI O'2 ) 
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FIG.  B1 5 ( a ) 


FIG.  B15 ( b ) 
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HWT  2-0-2/1  4  2/2  SUPPLY  NO.  2  (200/300) 
CORES  APRIL  1969  32  MONTHS  SERVICE 

TEST  TEMPERATURE  i  10  F 


RYE  FAILURE  STRESS 
AYE  FAILURE  STRAIN 


253.54 

0.00257 


SYMBOL 


IECEN0 

SAMPLE  NUMBER 

26  L.C. 

24  U.C. 

25  U.C. 
24  L.C. 
30  L.C. 

26  U.C. 


t.OO 


- 

0.10  0.20  0.30  0.40 

STRAIN  .  IN/IN  (X10-2) 


0.50 


FIG.  B1 5 ( c ) 


FIGURE  B15 :  Supply  No.  2 


STRAIN  ,  IN/IN  (X10-2 ) 

FIG.  B15 ( d ) 


32  Months  Service 


TENSILE  STRESS  .  PSI  TENSILE  STRESS  .  PSI 
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HHT  2-0-2/ 1  «,  2/2  SUPPLY  NO.  2  (200/300! 
CORES  APRIL  1969  32  MONTHS  SERVICE 

TEST  TEHPERATURE  i  0  F 

RVE  FBI LURE  3TRES3  •  301.71 

BVE  FBI LURE  STRAIN  »  0.00102 

LEGENO 

SYMBOL  SAMPLE  NUMBER 

17  L.C. 


O.OS  0.10  0.1S  0.20 

STRAIN  .  IN/IN  (X10-2) 

FIG.  B 1 6 ( a ) 


0.23 


FIG.  B1 6 ( b ) 
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HWY  2-0-2/1  «,  2/2  SUPPLY  NO.  2  (200/300) 
CORES  RPRIL  1969  32  MONTHS  SERVICE 

TE3T  TEMPERATURE  i  -10 F 


RYE  FAILURE  3TRES3 
BVE  FAILURE  STRAIN 
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4- 


4- 


366.33 

0.00031 


LEGENO 

SAMPLE  NUMBER 

27  L.C. 
30  U.C. 
16  U.C. 
16  L.C. 


0.03  0.10  0. 15  0.20  0.2S 

STRAIN  .  IN/IN  (X10-2) 

FIG.  B16 ( C ) 

FIGURE  B16 :  Supply  No.  : 


FIG.  B 1 6 ( d ) 

32  Months  Service 


TENSILE  STRESS  .  PSI  TENSILE  STRESS  .  PSI 

cp. 00  100.00  200.00  500.00  400.00  |  500.00  _  SOO.CO  700.00  Jl.CO  100.00  200.00  300.00  400.00  500.00 
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HHT  2-0-2/ 1  i,  2/2  SUPPLY' HO.  3  (200/300) 
CORES  RPR] l  1969  32  MONTHS  SERVICE 

TEST  TEMPERATURE  ,  20  f 

RVE  FAILURE  STRESS  •  IBS. HI 

RVE  FAILURE  STRAIN  .  0.00238 

LEGEN9 

symbol  sample  NU-stn 

34  U.C. 

33  U.C. 

31  L.C. 

3!  U.C. 

HS  U.C. 


0.10  0.20  0.30  0.40 

STRAIN  .  IN/IN  (X10‘2 ) 

FIG.  B17 ( a) 
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FIG.  B1 7 ( b ) 


HWY  2- 0-2/1  K  2/2  SUPPLY  NO.  3  (200/300) 
CORES  APRIL  1969  32  MONTHS  SERVICE 

TEST  TEMPERATURE  i  10  F 


me  FAILURE  ST1FFNES3  ■ 
me  FAILURE  STRAIN 


183207 

0.00243 


SYMBOL 


LECENO 

SAMPLE  NUMBER 

37  L.C. 
36  L.C. 
40  L.C. 
40  U.C. 
32  U.C. 
35  L.C. 


cb.00  0.10  0.20  0.30  0.40  0.50 

STRAIN  .  IN/ IN  (XI O'2 ) 


FIG.  B17 ( C ) 


FIG.  B17 ( d) 


FIGURE  B17 :  Supply  No.  3-32  Months  Service 
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TENSILE  STRESS  .  PSI  TENSILE  STRESS  ,  PSI 
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Ht-JT  2-0-2/1  4  2/2  SUPPLY  S3.  3  (200/300) 
CORES  APRIL  1969  32  MOUTHS  SERVICE 

TEST  TEMPERATURE  i  0  F 


RVE  FAILURE  STRESS 
RVE  FAILURE  STRAIN 


270.33 

O.COObS 


SYMBOL 


LEGEND 

SAMPLE  NUMBEft 
S6  U.C. 
37  U.C. 
44  U.C. 

41  U.C. 

42  L.C. 

42  U.C. 

43  U.C. 
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STRRIN  .  IH/INIX1Q-2) 

FIG.  B1 8 ( a ) 
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FIG.  B18 ( b ) 


-t- 


-+- 


HUT  2-D-2/1  4  2/2  SUPPLY  HO.  3  (200/300) 
CORES  APRIL  1969  32  MOUTHS  SERVICE 

TEST  TEMPERATURE  i  -10 F 


RVE  FAILURE  STRESS 
BVE  FAILURE  STRAIN 


270.24 

0.00073 


SYMBOL 


LEGEND 

SAMPLE  NUMBER 

38  U.C. 
33  L.C. 
35  U.C. 
38  L.C. 
32  L.C. 
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STRRIN  .  IN/INIXIO-2) 
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FIG.  B 1 8 ( c ) 

FIGURE  B18:  Supply  No 
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HUY  2-  0-2/1  4  2/2  SUPPLY  NO.  3  (200/300) 
CORES  APRIL  1969  32  MONTHS  SERVICE 

TEST  TEMPERATURE  i  -10 F 


flVE  FAILURE  STIFFNESS  ■ 
AVE  FAILURE  STRAIN 


673112 

0.00073 


SYMBOL 


LEGEND 

SAMPLE  NUMBER 

38  U.C. 
33  L.C. 
35  U.C. 
S3  L.C. 
32  L.C. 


0.05  0.10  0.15  0.20 

STRRIN  .  IN/ IN  (XI 0"2 ) 


0.25 


FIG.  B18 ( d) 


32  Months  Service 


•  ■ 

TENSILE  STRESS  ,  PSI  TENSILE  STRESS  ,  PSI 
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HNY  2-0-2/ 1  L  2/2  SUPP.  NO.  1  (200/300) 
CORES  (U.C.)  JUNE  1969  34  MONTHS  SERVICE 
TEST  TEMPERATURE  i  0  F 


fiVE  FAILURE  STRES3 
flVE  FAILURE  STRAIN 


522.32 

0.00042 


STHSOL 


LEGEND 

SflNPU  NUH8EA 
K1BC  7  U.C. 
NBL  5  U.C. 
NBC  4  U.C. 
HOC  2  U.C. 
NBC  1  U.C. 
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HWT  2-0-2/ 1  i.  2/2  SUPP.  NO.  1  (200/300) 
CORES  (U.C.)  JUNE  1969  34  MONTHS  SERVICE 
TEST  TEMPERATURE  i  0  F 


nv£  fit i lure  stiffness  - 

RVE  FAILURE  STRAIN 


2674948 

0.00042 


STMuOL 


LEGEND 

SAMPLE  NUMBER 
NBL  7  U.C. 
N8L  5  U.C. 
NBL  4  U.C. 
NBL  2  U.C. 
NBL  1  U.C. 


hi.oo 


O.OS  0.10  '  0.15  0.20 

STRAIN  .  IN/IN  (XI 0-2) 
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FIG.  B19 (a) 


FIG.  B19 ( b ) 


FIG.  B19 ( c ) 


FIG.  B19 ( d ) 


FIGURE  B.19 


'uoply  No .  1  -  34  Months  Service 
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B  2  0 
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STRAIN  .  IN/IN (X10-2  ) 


FIG.  B20 (a) 


FIG.  B  2  0 ( b ) 


FIG.  B 2 Q ( c  ) 


FIG.  B  2  0 ( d ) 


FIGURE  B20 :  Supply  No.  2  -  34  Months  Service 


TENSILE  STRESS  .  PSI  TENSILE  STRESS  .  PSI 

103.00  200.00  300.00  400. CO  500.00  600.00  700.00  41.00  100.00  200.00  300.00  400.00  500.00  600.00  700.00 


B21 


FIG.  B  2 1 ( a ) 


FIG.  B21 (b) 


FIG.  B21(c) 
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i  0  F 
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c  3TP.0IN  *  0.00104 

LEGF.N0 

STMBOL 

Sflin  E  MJMBER 

a 

SOL  10  L.C. 
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SBL  7  L.C. 

A 
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FIG.  B  2 1 ( d ) 


FIGURE  B21 :  Supply  No.  3  -  34  Months  Service 


APPENDIX  C 


COMPUTER  PLOTS  OF  TENSILE  STRESS  VS.  STRAIN 
AND  STIFFNESS  OF  MIX  VS.  STRAIN  FOR 
WATSONVILLE  AGGREGATE 

CALIFORNIA  ASPHALT  AND  KNEADING  COMPACTION  METHOD 

(ANDERSON ,  1968) 
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FIG.  Cl (a) 
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1/2  INCH  MAX  HEO  GRAD.  WATSONVILLE  FIGG. 
SUPPLY  NO.  '1176  (85/100  PENETRATION) 

TE3T  TEMRERATLFIE  .  10  F 

RVE  FAILURE  STIFFNESS  -  S5737 

RYE  FAILURE  STRAIN  -  0.00062 

LEGEND 

STROa  SAMPLE  NUMBER 
G3  01  CAL 
G8  02  CAL 
GQ  OS  CAL 
GO  ON  CAL 
G3  OS  CAL 


0.00 


0.02  0.CI  0.06  0.09 

STRAIN  ,  IN/IN  (X10-1  ) 
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FIG.  Cl ( b ) 
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1/2  INCH  MAX  MED  GRAD.  WATSONVILLE  AGG. 
SUPPLY  NO.  '1176  (85/100  PENETRATION) 

TEST  TEMPERATURE  i  30  F 

AVE  FAILURE  STIFFNESS  -  162833 

AVE  FAILURE  STRAIN  -  0.00472 

LEGEND 


0.02  0.09  0.06  0.03 

STRAIN  .  IN/IN  (Xt0-») 


o.to 


FIG.  Cl(c) 


FIG.  Cl(d) 


FIGURE  Cl 


Golden  Bear 


(85/100  Penetration) 
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FIG.  C2 ( a) 


FIG.  C  2 ( b ) 
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1/2  INCH  MAX  MED  GRfiD.  WATSONVILLE  RGG. 
SUPPLT  NO.  11176  (85/100  PEHETfflT ION) 
?ERTTUr£  «  0  F 


RVE  FAILURE  3Tr£33 
fiVE  FAILURE  3TPPIN 


593.73 

0.00110 


37 M80L 
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SfihPLE  IJJTOEfl 
08  16  CAL 
03  17  ca. 
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1/2  INCH  MAX  HEO  CRAO.  WATSONVILLE  RGG. 
SUPPLT  NO.  >1176  (85/100  PENETRflT  I  ONI 
TEST  TETPEPfiTlfiE  t  0  F 


0.05  0.10  O.tS  0.20 

STRAIN  .  IN/IN  (X1Q-*) 
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FIG.  C2(c) 


FIG.  C  2 ( d ) 


FIGURE  C2 


Golden  Bear  (85/100  Penetration) 
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1/2  INCH  MAX  HEO  GRAD,  WATSONVILLE  RGG. 
SUPPLY  NO.  5268  (85/100  PENETRATION) 

TEST  TEMPERATURE  i  40  F 

RVE  FAILURE  STIFFNES3  -  42908 

flVE  FAILURE  STRAIN  -  0.00776 

LEGEND 

STRB0L  SAMPLE  HUMBER 

.  sm  oi  ca 

.  SM  02  CPL 

.  sm  os  ca 

♦  sm  cn  ca 

,  sm  os  ca 


Oloo 


0.02  0.04  0.03  0.03 

STRAIN  .  IN/ IN  (XI 0'1  ) 
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FIG.  C 3 ( b )  • 


FIG.  C 3 ( c )  ' 


FIG.  C3(d) 


FIGURE  C 3 


Santa  Maria  (85/100  Penetration) 
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1/2  INCH  MAX  MED  GRAO.  WATSONVILLE  PGG. 
SUPPLY  NO.  5268  (05/100  PENETRATION) 
test  temperature  ,  ?o  r 

RVE  FAILURE  STRESS  *  357.41 

RVE  FAILURE  STRAIN  *  0.00470 


SYMBOL 


LEGENO 

SAMPLE  NUMBED 
SH  11  CAL 
SH  12  CAL 
SH  13  CAL 
SH  14  CAL 
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1/2  INCH  RAX  HLD  GRAO.  WATSONVILLE  RGG. 
SUPPLY  NO.  5268  (05/100  PENETRATION] 
TEST  TEHPEPATLRE  I  20  F 


AVE  FAILURE  STIFFNESS  « 
RVE  FAILURE  STRAIN  - 


1110273 

0.00470 


LEGENO 

STHEML  SAMPLE  NUMBER 
SH  It  CAL 
SH  12  CAL 
SH  13  CAL 
SH  l'l  CAL 
SH  IS  CAL 


l.co 


0.02  0.04  0.05  0.03 

STRAIN  .  IN/ IN  1X10-’) 


0.13 


FIG.  C4 ( a ) 


FIG.  C4(b) 
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1/2  INCH  MAX  RED  GRAO.  WATSONVILLE  AGO. 
SUPPLY  NO.  5268  (85/100  PENETRATION) 
TEST  TEMPERATURE  i  0  F 


RVE  FAILURE  STRESS 
RVE  FAILURE  STRAIN 


543.11 
O’.  00230 


LEGENO 

SYMBOL  SAMPLE  NUMBED 
SH  16  CAL 
SH  17  CAL 
SH  13  CAL 
SH  13  CAL 
SH  20  CAL 

- ( - 1 - 
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FIG.  C4(c) 
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1/2  INCH  HA<  HE0  GRAO.  WATSONVILLE  AGG. 
SUPPLY  NO.  5268  (05/ 100  PENETRATION) 
TEST  TEMPERATURE  i  0  F 


AVE  FAILURE  STIFFNESS  ■ 
AVE  FAILURE  STRAIN 


435515 

0.00230 


SYMBOL 


LEGENO 

SAMPLE  NUMBED 
SH  16  Crt. 
SH  17  CAL 
SH  18  CAL 
SH  19  CAL 
SH  20  CAL 


1.00 


0.02  0.04  0.06  0.03 

STRAIN  .  IN/ IN  (XI 0'1  ) 


o.io 


FIG.  C4(d) 


FIGURE  C4 :  Santa  Maria  (85/100  Penetration) 
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STRAIN  ,  IN/IN  IX 10-2) 
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1/2  INCH  MAX  MED  GRA0.  WATSONVILLE' AGG. 
SUPPLY  NO.  6571  (60/70  PENETRATION) 

TEST  TEMPERATURE  i  40  F 

AVE  FAILURE  3TIFFNES3  *  133400 

flVE  FAILURE  STRAIN  •  0.00437 

LEGENO 


STHBOL 


SAMPLE  NUMBER 
CH  613  CHL 
ch  614  ca 

CH  615  CHL 


0.10  0.20  0.30  0.40 

STRAIN  .  IN/IN  (XI 0-2) 


0.50 


FIG.  C5 ( a ) 


FIG.  C5(b) 
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FIG.  C  5 ( c )  FIG.  C5(d) 


FIGURE  C5 
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(60/70  Penetration) 


TENSILE  STRESS  .  PSI  TENSILE  STRESS  .  PSI 
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1/2  INCH  MAX  MEG  GRA0.  WATSONVILLE  AGG. 
SUPPLY  NO.  6571  (60/70  PENETRATION) 
TEST  TEilPERflTUnE  i  20  F 


OVE  FBI  LUBE  STIFFNESS  > 
BVE  FBILUlE  STHOIH 


560726 

0.00177 


STHBOL 


LEGEND 

SAMPLE  HJHSEft 
CH  buS  COL 
CH  606  CAL 
CH  607  COL 
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STRAIN  .  IN/IN  (X10-2  ) 


0.2S 


FIG.  C  6 ( b ) 
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STRAIN  .  IN/IN  (X10-2) 


FIG.  C6(c) 


FIG.  C6(d) 


FIGURE  C6 
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FIG.  C7 ( a ) 
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1/2  INCH  MRX  MED  GRRD.  NRTS0NVILLE  RGG. 
SUPPLY  NO.  6584  (05/100  PENETRRTION) 
TEST  TEhFEftTTUnE  *  40  F 


PVE  FPILURE  STIFFNESS  = 
RYE  FRILURE  STRAIN 


103057 

0.00537 


STHB0l 


LEGEND 

SAMPLE  NUMBER 
CM  813  CPL 
CM  814  CPL 
CM  815  CPL 
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STRAIN  .  IN/IN  (X10->  1 
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FIG.  C7(b) 
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HYE  FOILIJRE  STIFFNESS  = 
HYE  FPILURE  STRP1N 


24MGS8 
0. 00354 


STF.BCT. 


LEGEND 

SRMFLE  NUMBER 
CM  804  CPL 
CH  816  CPL 
CH  817  CPL 
CH  810  CPL 
CH  819  CPL 


t-DO 


0.02  0.04  0.06  0.08 

STRAIN  .  IN/IN  (X10-* ) 


o.io 


FIG.  C7(d) 


FIGURE  C7 :  Chevron  (85/100  Penetration) 
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1/2  INCH  MAX  NED  GRAD.  WATSONVILLE  AGG. 
SUPPLT  NO.  6584  (85/100  PENETFflTlONl 
TEST  TEMEEfinTURE  t  0  F 


RYE  FAILURE  STIFFNESS  -  1445773 

RYE  FAILURE  STRAIN  -  0.0C037 
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FIGURE  C8 


Chevron 


(85/100'  Penetration) 


• 

TENSILE  STRESS  ,  PSI  TENSILE  STRESS  .  PSI 

0.50  100.00  200. CO  300. OC  400. CO  500.00  600.00  700.00  J3.00  100.00  200.00  300.00  400.00  500.00  600.00  700.00 

^ a _ i _ 1 _ 1 _ -J - 1 _ I _ ■  -  ■  -l - 1  - A  - ■+.  .  - 1--  .  — i - i - - - 4 - ►  ■  -t  -4 - 1 - 1  t  4 - r 


C9 


-4- 


-4- 


-4- 


-4- 


1/2  INCH  MAX  MED  GRAD.  WATSONVILLE'  AGG. 
SUPPLY  NO.  6572  (120/150  PENETRATION) 

TEST  TEMPERATURE  ,  no  F 
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fig.  C9 ( a) 
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1/2  INCH  MAX  ME0  GRAD.  WATSONVILLE  AGG. 
SUPPLY  NO.  6572  (120/150  PENETRATION) 
TEST  TEMPERATURE  I  40  F 
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FIG.  C9(b) 
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FIGURE  C9 :  Chevron 
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1/2  INCH  MAX  MEO  C-RAO.  WATSONVILLE  AGG. 
SUPPLY  NO.  6572  (120/150  PENETRATION) 
TEST  TEMPERATURE  i  30  F 
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FIG.  C9(d) 
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FIG  .  C10(a) 
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1/2  INCH  MAX  MED  GRftO.  WATSONVILLE  AGG. 
SUPPLY  NO.  6572  (120/150  PENETRATION) 
TEST  TEMf’EraTUne  I  0  F 
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FIGURE  CIO:  Chevron 
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APPENDIX 
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COMPUTER  PLOTS  OF  TENSILE  STRESS  VS.  STRAIN 
AND  STIFFNESS  OF  MIX  VS.  STRAIN  FOR 
WATSONVILLE  AGGREGATE 
CALIFORNIA  AND  ALBERTA  ASPHALTS 
VARIOUS  MECHANICAL  IMPACT  COMPACTIVE  EFFORTS 
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SUPPLY  NUMBER  5268  4  WATSONVILLE  RGG. 
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FIG.  Dl(a) 

FIGURE  Dl:  Santa  Maria 
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(85/100  Penetration)  75  Blows 
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FIGURE  D2 
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FIG.  D3 ( a) 
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SUPPLY  NUMBER  5268  <  WOTSONVJU  E  AGG. 
COMPACTION  :  10  BLOWS 

TEST  TENPemruie  I  0  F 
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FIG.  D3(b) 


FIGURE  D3 
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FIG.  D4(b) 


FIGURE  D4 :  Alberta  Supply  No.  1  -'-75  Blows 


D3 


FIG.  D5(a) 

FIGURE  D5 :  Alberta  Supply  No. 
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WATSONVIllE  AGGREGATE 
ASPHALT  SUPP.  NO.  3  (200/300  PEN) 

TEST  TO'PcftVrUfie  t  0  F 
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FIG.  D6(b) 


FIGURE  D6 :  Alberta  Supply  No.  3 


75  Blows 


